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Simplified methods are presented for calculating the number of arrays of close-packed 
spheres having a repeat distance of n layers; the distributions over the 8 symmetry 


groups are also given. 


Belov [1] has given a complete list of close 
packings with repeat over n layers for values of n 
up to n‘= 12; there are 43 such packings for n = 12; 
for 13 layers and above, the problem appeared too 
difficult and consideration of it was deferred. The 
methods previously used have now been simplified 
somewhat and made more systematic; the solutions 
for n = 13 and 14 are obtained almost mechanically. 

The essential minimum symmetry of anyclose- 
packed array is re) = P3m (I); the packings for this 
purpose may be denoted by h and c, of whichh di- 
vides the formula of h and c into two symmetrical 
halves and denotes the passage of a symmetry 
plane through the corresponding layer of spheres 
perpendicular to the packing axis. The centers of 
symmetry lie either in the c layers, in which case 
the formula is divided into two symmetrical halves, 
or between two cc or hh layers, if the correspond- 
ing pair divides the packing formula into two sym- 
metrical halves. 

The minimum repeat distance in a packing 
formula consisting of h and c corresponds either 
to a complete repeat distance in the array of 
spheres or to a half or third of the latter. That 
minimum repeat may contain an odd number of h, 
in which case it is always half the full repeat dis- 
tance, and the array is characterized by 
a 63 axis (not a triad one), and sc belongs to group 
oe = P63mce (II) if there are no additional symmetry 
elements. Alternatively, it may contain an even 
number of h , in which case the full repeat is de- 
termined by the number of c in the minimum re- 
peat, which is found as follows: we count the ¢ be- 
tween each pair of #, the corresponding numbers 


ait 


Detailed calculations are presented for n = 13 and 14. 


being formed into a sum with alternating plus and 
minus signs (not excluding zeros, which correspond 
to runs of h). The algebraic sum may be a mul- 
tiple of 3, in which case the repeat distance of the 
packing coincides with the minimum repeat dis- 
tance of the formula, the space group being ch = 
P3m in the absence of additional symmetry ele- 
ments. Alternatively, the sum may not be a mul- 
tiple of 3, in which case the full repeat is three 
times the minimum, and the space group is C3y = 
R38m (I) in the absence of additional symmetry 
elements. 

Additional centers of symmetry are possible in 
all three cases: Ch, = P3m becomes Diy = Pém = 
P3m (IV), and fs = R38m becomes Dig = Rém = 
R3m (V), but those in group Ch. = P63mc also gen- 
erate horizontal symmetry planes, so the space 
group becomes Dah = P63/ mmc (VI). 

Horizontal symmetry planes may be added in- 
dependently only to C$, = P3m, which gives Dip = 
= P3/mm2 = Pém2 (VII). 

The packing formula may not contain hf at all, 
in which case we getcubic close packing ... cec... 
with symmetry O; = Fm3m (VIII). 

These features give us the following system- 
atic way of deducing all packings of repeat distance 
n. We seek all groups with an even number p of h 
and with sums of the c that are multiples of 3. We 
first draw up all combinations of even numbers of 
h possible for the givenn. For instance, for n= 
1S 

Dae 
is Bill Pear ii beeag Be 


6, ot, 42,°4117 733,321, 3111, 222) 
2214 e2i2dee ddd, 114411), 
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8, 71, 62, 614, 53, 524, 5114, 
431, 422, 4244, 4121, 
3311, 3131, 3221, 3242, 32411, 
31211, 2222, 22241, 21212, 

10, 91; 82, 814, 73, 721, 64, 634, 
622, 55, 544, 532, 442 


44, 


41111, 


Sony 


Some c must come between the h, and here 


there cannot be less than 3, so p= 10 


is the limit. 
For each combination of the h we list the dif- 


ferent numbers of c, and check that the algebraic 


sum is divisible by 3. 
have for the 321 combination that 


te tey ieee SOP eed, ot 
Sieg PG prec eae ee ee ee 
See Ba See Vc eB = es a 
Ce ey ea ee pe ee erin) 

eto ee eee yee aed Seal ony as ie 


Bail = See ced 
= eee 
yee eee beeen 
Tie aye) 

So pe ee SESE 


Five formulas for the 321 combination have 
multiples of 3, which in letter form are 


de | AAhchhcheecce|.,., 2) ...|hhhchheccchee | -- ., 
3) . | Ahhechhecchcc | ees 4) sieve Wennecenhechee | eeey 
5) | hhheccchhchee | Shep t/= 


The combinations of the h written as numbers 


may have symmetry planes, in which 


case inser- 


tion of the c gives us two formulas that are, re- 
spectively, mirror images, and such arrays are to 


be taken as identical. For instance, 
for the following in combination 141: 


. | Achhhhcecchec |... and . 


All the combinations of the h are examined in 


this occurs 


. | echcccchhhhch | Ae 


this way to give us the total number of n-layer ar- 
rays, but some of these are only multiples of ar- 


rays of shorter repeat distance. 


out by specifving the minimum repeat distances, 


which may [1] only be n/2 or n/3, subject to cer- 


For instance, for n = 13,we 


These are struck 


tain conditions. Of course, no minimal periods are 
possible at all for n = 138, while for n = 14 they each 


have 7 members and correspond to independent 7- 


layer arrays provided that the number of h is even 


in each half. Finally, the arrays are 


assigned to 


their space groups by identification of centers of 
symmetry and mirror planes, in accordance with 


the above rules. 
This gives the number of 13- ae 


63, of which 21 have symmetry Dg = = P6m = P3m:! 


1) | eccechchecece i 2) | echecccccchee |, 3) 
4) \ccechhchheece |, 5) | chhccccccchhe |, 6) 
7) | echhececchhec | é 8) | heechccchccch ee 9) 
10) | hhchccccchchh | ’ 11) | chhheeccchhhe | ’ 12) 
13) | neehchchchech |, 14) | chcchhchheche ies 15) 
16) | fhhechchcchhh| , 17) | Acchhhchhhcch| , 18) 
19) | hhhhhccchhhhh , , 20) | hhhchhchhchhb| - 24) 


r arrays as 


hhececccccchh 
ecchhccchhcce 
checchchccchce 
hehchccchchch 
chchhccchhchc 


hchhhhchhhhch| 


| hchhchchchhch| , 


and 42 have symmetry Oe = 


22) hhhecchcceccc | , 
20) hhechccchcccc : | 4 
28) |hhhhhccheccce | , 
31) |hhhchhcheecce | , 
34) | hhhechhecchce ies 
37) | hhhchcheccche | , 
40) | hhchchchchcce | ; 
43) hhchecchheche | ; 
46) |nhhhhhhchecce | 3 
49) |nhhhhchhecche | , 
52) |hhhhchhhchccc | , 
55) |hhhhchhehekce | , 
58) |ahhchhhchhece | , 
61) |Ahhechhchhchc | , 
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P3m: 
23) | hhehcccheecce te 
26) |achchechcecce | 
£9) |nhhhchechcece | , 
32) | Ahhchhcccehce | , 
35) |hhhecchhcchce | , 
38) hhhehcchehccc | , 
41) | hhchchechchee li; 
44) |hhchhchcchcce | , 
47) |hhhhhhchehccc | , 
50) | Ahhhhcechhohe | ; 
53) |Ahhhcecchhhchc | , 
56) |hhhhchhcchchc | , 
59) |hAhchhchhcchc | , 
62)_ |hhhhhhhhchhcc! , 


24) hhchececcchce | - 
, 27) |ehechcchcece | , 
30) | hhhhechechcce ee 
33) |hhheccchhchce | , 
36) |hhhchchehcece | , 
39) |nhhcheechchce | , 
42) hhchchhchccecc | ' 
45) | nhechhchechcc | 
48) |hhhhhchhhecce | , 
51) |Ahhhhechhcche | , 
54) |hhhhechhhchcc | , 
57) |nhhhechhchche | , 
60) hhhchhechhche H - 
63) |hhhhhhchhhhcc! . 


The 14-layer arrays number 126, of which 7 


have symmetry Dép = P6;/mme: 


L) 
|chhhhhechhhhhc |! , 
6) 


3) 


| ecchcceccchcce 
4) 
| hhchchhhhchchh | 


2 


is 2) 
| chehchechchche | ; 
7) 


| echhhcccchhhcc | , 


!hechechhcchcch | , 


5) 


|hchhhchhchhhch | . 


One has symmetry Céy = P63/ mc: 


8) 


| hchheechchhece he 


Twenty-four have symmetry Din = P6m2:° 


a) Rheceecheccecchh ‘ 
2) hhhecccheccchhh b 
5) Ahecchchchccchh , 
8) Rhhhecchecchhhk, 
21) Rchhhechechhhch, 
24) hchhchchchchhch, 
7) Achchhhhhhhchch , 
SO) Ahhhhhchchhhhhh , 


i 
1 
4 
2 


Twenty-four have 


33) | ceceeechhececce 
36) | cheecchheecche 
39) échchceecchchcé , 
éhhecchchecchhe , 
écchhchchchhcceé , 
hhhecchhecchhh ' 


\ |echhchhhhchhcc | , 


; 


4) | hchhechhcchhch | , 


de) 
49) 
yi 


el 


10) Aehececheccchch , 


11) hechccehecchcch’ ‘ 


) Ahchccchecchchi, se hae Sakweaes 
16) Rehheccheechhch, 17) Rchchcchechchch , 
19) Ahheechhhecchhi, (0) Ahhechchchcchhh, 
22) Rhechhehchhechh, 23) Hhchchchchchchh, 
25) Ahhhhcchechhhhk, 26) fhhhechhhcchhhi, 
28) Ahhchchhhchchhh, pe Bhchhhchchhhch, 
34) HhhAhchhhchhhhh, 32) Ahhhchhhhhchhhh. 


symmetry D3q = P3m:4 


‘ i échccecceccchcé , 


37) éhchecccecchché , 


40) | hhechcccchcchh | , 


43) ) | chhehcecchchhe | 
46) | nchechechecheh | ; 
49) écchhhhchhhhecé, 

52) | chhechhhhccAhe | ; 
09) |Ahhchchhchchha| , 


39) | eccchchhchccccec | , 
38) | chchccccechche | ; 
41) écechhhchhheccé , 
44) échhechchechhcé , 
47) \hhhhececcchhhh | , 
90) | cchhhhcchhhhce | , 
3) chhchhecchhchhe , 


5 
06) | chhhchhhhchhhe! ; 


’ 


‘The dots before and after the formulas are omitted. The centers 
of symmetry are found each time from the middle letter and the 


end strokes [1]. 


The rules show that group Déh has a center of symmetry at each 
end stroke and between the two central letters, while symmetry 
planes pass through letters 4 and 11 (see [1]). 


3The first letter (accented) and the ttanslationally identical 15th, 
and also the 8th, have symmetry planes [1]. 


“A center of symmetry lies-at the start (and end) of the period,and 


another in the middle. 
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and 70 have symmetry C4, = P3m:° 


105) | hhhhchhchecche | , 106) | hahhchchhechcc ; , 107) hhhhechechhchce| , 
108) | hhahchchchceche |; , |()9)| hhheechhhchhcc| , 110)| hhhchhhchcchec | , 
07) | Rhhecheececece le 58)! hhchechccecccc | , 59) hhchcceechccce } 111) | hhhechhhechche | , L12) | hhhchechhheche | Fecal G Ws) | hhhchhchhchece | ’ 
G0) hhechecheccece | , 61)! hhechcecechece | , 62) | Ahecechechecce | , 114) | ahchhechchhce | , 115) | hhhchhchchchcec |, 116) !hhhchhehcchehe | , 
63) | Acheehecechece | , (4)|hAhAhcheccecce | , 69) |hhhhchchececee |,  117)| Ahhchchhchechc | , 18) |hhhechchhchche | , 119) | hhchchhechchhe | , 
(6) Ahhhchceechece | | 67)|hAhhechchcccce | , 65) | hhAheechchecee | , (20) | phhhhhhhcheche; , 121) |hhhhhAhchhchec| , 122) |hhhhhhchhhechc | , 
9) | hhhechehecheee | , 70)|hhhchhecehcece |, 71)! hhhchhecceeche | , 123) | Ahhhhchhhhchee| , 124) | hhhhhechhhchhc| , 125) ‘hhhhhchhhchchc| , 
1) Ahhcchhechecce } , 73) | hhhcchhecccche | , 74) | hhhecchhchccce | F 126) ‘Rhhhchhhchhche | ‘ 
75) |Ahheechheccche | , 76) | hAheecchhecche | , 77) | Rhhchchcechece | , 
78) | Ahhececchheche | , 79)| Ahhehcchccchcc | , 80) | hhhececechhche | , 
81) ‘Rhchhehcheccce | , 82) | hhchhcheccchcc | , 73) \\Ahchhechchcccc | , LITE RA TURE CITED 
HEY) Bhochhekeheece: |, 85)| hheahhechchese | ,. 96) Miccehhehehcee.|-— 4. NW Belov. Structures.of Jonic Crystals and 
87) | hhecechhechche | , 88)| hhchechhechcce | , 89) | hhchchchecchcc | , : . 
90) | Ahchchcchceche | f 91) | hhchechchechcc | x 92) | hhchecchchchcc | , Metallic Phases [in Russian] (1947). 
3) | hhhhhhchheeeec| , 94)| hhhhhhechheecec! | 5) | hhhhhhcheechee i 
6) | hhhhhchheehece |, 97) | Ahhhhechhcheec| , )8) | hhAhhehchechce | , 
99) | Ahhhhchccheche | , 100) | hhhhchhchhecce ; , 101)| hhhhchhechhece | , There are no additional symmetry elements; in each group x= 


102) | hhhhchheechhee | , 103) | hhhhchhheechee | , 10 4)| hhhhcchhhecche | , 0 (mod, 3). 
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A method previously described is applied to give working formulas for the variance 
in the potential as deduced via Fourier series from electron-diffraction data, to- 
gether with the errors in the peak heights arising from termination and the vari- 


ances of the atomic coordinates. 


Fourier series derived from experimental am- 
plitudes do not give the true distribution of the lat- 
tice potential on account of experimental error, and 
also of series cutoff. The accuracy of the poten- 
tial and of the atomic coordinates may be deduced 
by a method similar to that given previously for the 
electron density [1], which is as follows. The 
atomic f, curves [2] are used to compute the inte- 
gral characteristics G of the potential ~ within the 
atoms for several B (thermal-vibration parameter): 
the mean-square potential, the potential at the cen- 
ter of the atom, and the mean square of the first 
derivative of the potential. These can be repre- 
sented via a universal relation to the atomic num- 
ber: 


G= kl" 4, (1) 


in which k and a are governed by the Debye factor 

B characteristic of the crystal, q being determined 
by the cutoff; q = 1 when there is no cutoff, decreas- 
ing as the cutoff increases. The factor q reflects 
the loss of Fourier terms by cutoff. Sometimes it 
is necessary to take account of the contribution to 

a characteristic from terms in the cutoff region, in 
which case the q of (1) is replaced by (1 —q). 

The accuracy of a Fourier series for the po- 
tential may readily be evaluated via the character- 
istics of the atoms in the unit cell, which are easily 
found from (1); here we need to know the extent of 
the experimental error, for which purpose we may 
use the factor 
— ZI®o|—] 9, | 


b x| Do] 


(2) 


Only an approximate value for B is needed in or- 
der to evaluate the accuracy of the series; B ~ 1 for 
most silicates and other materials with strong 
bonds and high melting points, 1-2 for ionic crys- 
tals and metals, and around 4 for organic crystals. 
Heavy atoms may reduce B to 0.5-1. 


Error in the Potential 


The Fourier series here is constructed from 
®hk] tAPpK] , With A®pK] = b&p_K] as the errors 
in deducing the amplitudes and b the mean relative 
error given by (2). The standard deviation in the 
potential is [1] dependent on b and is given by 


aes (3) 


in which G, is the appropriate characteristic for the 
atom in question, the sum being taken over all the 
atoms in the cell, whose volume is 2. The cutoff 


error is 
Vad =V A —9)36,,/2. (4) 


The over-all standard deviation for a three- 
dimensional series on the basis of the G, of (1) is 


then 
+A =V (09 +1—@)kBz2/0. (5) 


Vad =Viek 

The Q of (5) is replaced by S (cell area) for 
the two-dimensional series used in projections. 
Table 1 gives k, @, and q as functions of B for the 
two cases. 


XG; ES, 
vu 
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TABLE 1. Values of k, a, and q for (5) 


q for (sin@/A)o + 10°8 


Series B h 
dimensions 
Three 4 0.0557 
2 0.0490 
3 0.0443 | 
4 0.0399 
Two 4 0.0860 
2. 0.0840 
3 0.0815 
4 0.0790. | 


0,5 0,6 0.7 0.8 

39 | 0,780 0,861 0.913 | 0.948 
30 | 0,920 0.960 0.982 | 0.995 
32 1) 0,980 0,380 |°0,995 | 4 

vw | 0,965 0.987 OPS97>) 4 

Zip Oe Oui: 0.901 0.978 | 0.989 
24 | 0.960 0.983 ().994 | 4 

2! ().980 0.993 0.998 | 4 

AWW AGE Silent! 0,997 t 1 


TABLE 2. Values of k, a, and q for Use in (1) for g(0) and w*(0) 


q for (sin@ /A)g + 10° 
B hk r a 


a 


(: clewerssarlt Oc8t 

2» | 0.495 | 0.79 
? (0) OU ag MNO 77 
A h0. 3494 ad” Oe8 

t | 0,480 | 0.75 

; » | 0.403 | 0°72 
9’ (0) 3 | 0.362 |» 6 70 
4d 6.395 | 0.68 


The G are calculated from the fe given in [2], 
expressed in p-units | so the error given by (5) 
will be in p/ A (three dimensions) or p/ A? (two); 
conversion to volts (or volt-A for projections) 
simply involves multiplication by 114.5 [3]. The 
peak heights are of the order of a few hundred 
volts, while Vag is around ten volts or so in or- 
dinary structural studies, with b © 0.2-0.3 and 
(sin 0/A), © 0.7. The G given by (1) and Table 1 
may be used to normalize Fourier series for the 
potential, i.e., to reduce them to the absolute 
seale [1,4, 5]. 

The V Ag? of (5) gives an estimate of the error 
in @ as an average over the whole volume (or area) 
of the cell; but sometimes we need to know how the 
cutoff affects the peak heights (0), as its effect is 
greatest for these. This ~(0) = G, is the second in- 
tegral characteristic that can be put in the form of 
(1), as can the analogous quantity ~'(0) for projec- 
tions. The q in (1) indicates the loss in peak height; 
Table 2 gives the parameters of 1 for the cases of 
(0) and ¢'(0). 

This and (1) together enable us to find the ex- 
pectation value for the potential at the center for 
any B. The peak heights in electron diffraction are 
proportional to Z%, with @ < 1, which explains 
why light atoms are more readily detected in 
the presence of heavy ones (2 > 1 for x rays). 


0.5 0.6 0.7 8 
33 0.45 0,90 0.64 
50 0.63 0.74 0.81 
60 0.73 0,83 0.89 
69 0.81 0.90 0,94 
i) 0,69 0.77 0.83 
72 0,81 ().88 0.92 
80 ().38 0.93 0.96 
86 0.92 0.96 0.98 


The peaks for the light atoms are readily 
seen, together with those from the heavy 
atoms, via Fourier series for the potential, on ac- 
count of the weak dependence on atomic number [6]. 
The q of Table 2 show that cutoff reduces the heights 
substantially; as for the V Ag? of (5), the k of Table 2 
calculated from (1) give the 9 (0) inp/A®, and the 

¢'(0) in p/A*, while multiplication by 114.5 con- 
verts these respectively to volts and volt-A. 


Error in the Atomic Coordinates 


The general formula for the error Axj in a co- 
ordinate of atom j is 
/ (wei); 
ue (6) 


Ax; is directly proportional to the perturbation 
produced by the error and cutoff waves (i.e., to the 
first derivative of these) and is inversely propor- 
tional to the curvature of the peak for atom j (Second 
derivative at the peak). The AgY have a random dis- 
tribution, so the dAy/ 8x do also; the numerator is 
therefore to be evaluated as its mean square, as in 


j dA 
SUES 


‘A p-unit is the electron scattering for a proton for (sin 6/2) = 0,1- 
10° cm™ (2, 3]. 
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the evaluation of Ay from (3)-(5). 
ror introduced via JA¢/ 0x from errors in the am- 
plitudes may [1] be found from 


eS en 


) = |? 2G. / Q 4 


(7) 


while that introduced via the cutoff wave is, by an- 
alogy with (4), given by 


i c elena) a. (Q. 


Ox (8) 

The over-all standard deviation V (OA g/Ox) 
may be calculated from the sums of (7) and (8) by 
analogy with (5); but here we need a certain refine- 
ment, because the cutoff wave has a regular struc- 
ture, nota random one. The atom in question 
must lie outside spheres of radius about 2 A drawn 
around all the other atoms; but the cutoff wave and 
the first derivative of this are maximal within such 
spheres, so the quantity we need to find, Ax}; is 
only a fraction of that given by (8). It is found that 
the significant part of the mean square of the first 
derivative outside the sphere (as regards Ax;) may 
be found by multiplying the result from (8) by wee so 
it accounts for V4 ~ 4% of the total. Then, 


VC -V (eo Saag 


The second derivative may also be put as in 


(1): 


a (10) 
By substituting (9) and (10) into (6), we find, 
when B = 1-4, Eqs. (11a)-(11d) for particular 
values of the parameters Vk/k', @ and @' with the 
introduction of an effective error beff given by Eq. 
(12); because of cutoff, the magnitude of this error 
is always larger than the error b given by Eq.(2). 


As in (8), theer- 


For three-dimensional series 


for B=1 Ax, = 0.025 (bett/ VO) V ZF" zs", (11a) 
% 

for B=2 Az; = 0.045 (bese) VQ) Vir Wagar t e's) 

for B=3 Az; = 0.065 (betef VQ) V BZ} /Z}™, (110) 


for B=4 Aa; = 0,093 (bett/ VQ) V 22)" (23, (11d) 


be 1—q)/8 
bett= sea _ (12) 


where 

The q and q' are given by Table 3 and Fig. 1. 
The sum under the root in (1la)-(11d) is taken over 
all atoms in the unit cell; the atomic numbers a 
in the denominator being for the atom appropriate 
to Axj; Q is the cell volume. 

The low q' show that cutoff markedly reduces 
the second derivative, especially for small B, 
which increases Ax;. 

The analogous formulas for two-dimensional 
series are 


10 
(sin G/A)o 


Fig. 1. Relation of q (solid lines) 
and q' (broken lines) to (sin@/ A) 
and to B for three-dimensional 
series. 


TABLE 3. Values of q and q' for (12) for Three-Dimensional Series 


sin@/A, : 10° 

B 
0.4 0.5 0.6 0.7 0.8 0.9 1.0 
{ 22 36 49 62 72 80 87 
ate 2 39 57 72 84 91 95 98 
, 3 52 72 85 93 97. 99 100 
4 62 81 92 97 99 1CO 100 
n 3 6 10 16 24 33 49 
hoe 2 6 13 24 35 AT 7 67 
ay e40 3 43 2A 37 52 66 THI 86 
4 19 35 ih 1 He 79 89 94 
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TABLE 4, Values of q and q' for (12) for Two-Dimensional Series 


sin 0/Xg* Mone 


0.4 0.5 

L 4 57 

, 2 59 75 
oy 3 68 84 
4 77 90 

{ 8 14 

i 2 17 28 
1 %o 3 25 A() 
4 33 ot 


for B=1 Ax; = 0.057 (bets/ VS) / BZ}? /Z2", (18a) 
for B= 2 Ax; = 0.095 (lete/ VS) V 22} /22*, (13b) 
for B=3 Az, = 0.130 (lese/ VS) V 2Z5" jz", (180) 
for B=4 Ax; = 0.165 (beg VS) V 32)" 29” (134) 


Here, beff is still given by (12), but q and q' 
are different (Table 4 and Fig. 2); Sis the cellarea. 
It is to be appreciated that (13a)-(13d) are not ap- 
plicable unless atom j is represented by a separate 
peak. 

Equations (11) and (13) apply to Ax, Ay, and 
Az; Ax = Ay = Az if the cutoff is the same in all di- 
rections, but q and q' will vary with direction if it 
is not, as will the bef of (12), so the errors will be 
different (largest in the direction of greatest cutoff). 


Equations (11) and (13) give the inherent er- 
rors of the series for the potential. The errors in 
the atomic coordinates in addition include contribu- 
tions from errors in the cell constants and in the 
interpolations involved in locating the maxima. 
These additional sources of error are usually rela- 
tively trivial and increase the Ax; of (11) and (138) 
only slightly; reasonable intensity estimates (b © 
0.2) and optimal cutoff at (sin 0/A) - 107® = 0.6-0.7 
have been found by experience and by calculation by 
these methods or others [7,8] to give errors not in 
excess of 0.01 A (from three-dimensional series) 
for the coordinates of the heaviest atoms in a struc- 
ture, or 0.02-0.038 A for projections. The struc- 
ture may contain heavy, medium, and light atoms, 
in which case the errors for the medium atoms are 
much as for x-ray methods, those for the heavy 
ones are rather larger, and those for the lightones 
are reduced by factors of two or three. This again 
shows that light atoms are readily detected in the 
presence of heavy ones by electron diffraction. 
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Fig. 2. Relation of q (solid lines) 
and q' (broken lines) to (sin@/A)¢ 
and to B for two-dimensional 
series. 


Conclusions 


A universal relation kZ%q (in which Z is the 
atomic number) may be used to compute the accu- 
racy for two- and three-dimensional Fourier series 
in electron diffraction via integral characteristics 
for the distribution of potential in the atom (mean 
squares of the potential and of the first derivative, 
potential and second derivative of the potential at 
the center of the atom). Here, k, @, and q are pa- 
rameters, the last relating to cutoff; values are 
given for these for each of the characteristics listed 
above for several values of B, the temperature fac- 
LOSS 

These characteristics are used in working for- 
mulas for the accuracy of series for the potential 
and atomic coordinates, the reliability factor be- 
ing used as an estimate of the experimental error. 
The parameters are presented as tables and graphs. 
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We have already [1] deduced the possible struc- 
ture types with the closest packing of A and B 
atoms for the case in which each A atom lies ina 
cube-octahedron of B atoms, while each B atom 
has four A atoms in its coordination sphere. 

The same principle is used here: the closest 
packing allows only structure types in which each 
A atom is surrounded only by B atoms, while each 
B atom has at least one A atom in its coordination 
sphere. Detailed consideration is given to the case 
in which each B atom has three A atoms adjacent 
to it; then each B has '4 of an A atom, while each 
A atom has ,-12=4 Batoms. The composition 
is AB,. 

Consider a close-packed layer; if a B atom is 
surrounded only by B atoms in such a layer, it can 
have not more than two A atoms adjacent to it inits 
coordination sphere in the three-dimensional array. 

It may be shown that, if a B atom in a layer 
has three A atoms adjacent to it, there are Batoms 
each having four A atoms in the coordination sphere 
in space. Figure la shows that atoms of types 1 
and 2 must be A in cubic close packing, because 
any deviation from this would mean that atoms of 
types 3 and 4 would have no A atom adjacent in the 
layer. But then the B atoms of type 5 each have 
four A atoms in the coordination sphere in space. 
The same occurs for hexagonal close packing (Fig. 
1b). Atoms of type 1 in the layers above and below 
must be A, but then the B atoms of type 5 have four 
A neighbors each. 

Thus, a B atom has either one or two adjacent 
A atoms within a layer; in the latter case, the hex- 
agons of B atoms around the A atoms are joined at 
vertices or along sides. 
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Joints at vertices give rise to B atoms each 
having four A atoms in the coordination sphere. 
Figure 2 shows that the type 1 atoms cannot be B, 
because they cannot have more than one adjacentA 
within a layer, and hence not more than two in 
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Fig. 1. (a) Cubic close packing; (b) hex- 

agonal close packing. The centers of the 

hexagons in the middle layer are A atoms, 
while the vertices are B atoms. 
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Fig. 2. (a) Cubic close packing; (b) hex- 
agonal close packing. The centers of the 
hexagons are A atoms, the vertices being 
B, 


space. If type 1 is represented by A, each B atom 
of type 2 will touch four A atoms. 

Further, two junctions along sides that are not 
parallel inevitably imply junction at a vertex;atoms 
of type 1 or 2 (Fig.3) must be A, because any other 
results would mean that a B atom of type 3 would be 
surrounded only by B atoms. But then point 4 or 5 
will represent junction at a vertex, which we have 
shown above to be impermissible. 

There remain three types of junctions (Fig. 4); 
in the first (Fig. 4a) we have four B atoms per A 
atom; in the second (Fig. 4b), five; and in the third 
(Fig. 4c), six. A layer on average has at least 
four B atoms per A atom, and this occurs when 
there are only junctions of the first type. These 
may occur together with the second and third types, 


N. L. SMIRNOVA 


Fig. 3. The centers of the hexagons are A 
atoms, the vertices being B; the broken 
lines indicate hexagons of B atoms thatoc- 
cur if atoms 1 and 2 are of type A. 
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Fig. 4. Possible types of junction in a close- 
packed layer: junction along sides. 


superposition is possible in cubic close packing and 
in hexagonal (Figs. 6 and 7); both figures show the 
unit cells and the hexagons in the lowest layer, the 
numbers 0,1,2, etc., representing A atoms in 
layers 0,1,2, etc. The repeat distance perpendicu- 
lar to a layer is 15 layers for the cubic case and 10 
for hexagonal. The structure type in the first case 
is described as follows: a = aceV%; c = ac; space 
group 14/m; n = 2; A at 2(a); B at 8(h): x= 0.4, 

y = 0.2. This is known under the name Ni, [2] and 
has several representatives. The unit cell is mono- 
clinic in the second case (Fig. 8): a = V 4c} as 
b = ayV3, c= V 4a}, +c}; B = 67°49'; space group 


in which case the A :B ratio in a layer would be less 
than 1:4; to have four B atoms per A atom inspace, 
we must have either that each layer has four B 
atoms per A atom or, if there are layers withA:B > 
1:4, there must also be ones with A:B < 1:4. The 
first is impossible, so the second is as well; junc- 
tions of the second and third types are therefore 
eliminated, to leave the first type as the only one 
possible within a layer (Fig.5). The next layer of 
this type is superimposed in such a way as to give 
each B atom three A neighbors. Only one mode of 


Fig. 5, Array of A and B atoms in a layer for junc- 
tions of the type of Fig. 4a; the centers of the hex- 
agons are A atoms, the vertices being B atoms, 
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W012 Fig. 8. Unit cell of the possible struc- 

o§ xg ture type for AB, with hexagonal close 
kd of 08 packing, the numbers being the heights 
as fractions of the cell parameter: O— 
Nk = Vey (EY 18 re 


OTF 97 of SIM oF 1 
Fig. 6. Possible structure type for AB, with cubic 
close packing: x — layers 0,3,6,9,15; @— layers 
1,4,7,10,13; O— layers 2,5,8,11,14, 
C2/c; n= 4; A at 4(e):0y 4, y= Ves By at 8(/):xyz, 
x= 0515 Viewsat Olde Brat 80): x= 0.3, y= 
*/,, z= 0.3. The literature appears to carry no 
evidence for representatives of this type. 

There are thus two structure types possible for 
AB, subject to the above conditions. 

We are indebted to Academician N. V. Belov 
and Professor G.B. Bokii for valuable advice and 


comments. 
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Three compounds, CoSb, CoSb,, and CoSb3 are established in the cobalt—antimony 
system. The compound CoSh; crystallizes in the cubic system with a = 9.016 kx. 


According to data of Podkopaev [1], Lewkonja 
[2], and Losev [3], two compounds, CoSb and CoSh,, 
can crystallize in the Co—Sb system. In the range 
70 to 95 wt.% Sb, between the peritectic and eutec- 
tic horizontals (894 and 618°C, respectively, Fig. 1), 
Losev [3] noticed a thermal effect at varying tem- 
peratures (~770°C); he ascribed this to a polymorph- 
ic transformation of the compound CoSh,. Since 
two other workers [1,2] had failed to observe this 
phenomenon, it was inserted into the phase diagram 
as "possible" only [4]. 

Our own microscopic and x-ray structural 
studies of several cobalt—antimony alloys offered 
no confirmation of the polymorphism of compound 
CoSb,. In our view, the thermal effect at ~770°C 
should be attributed not to a polymorphic trans - 
formation of the compound CoSh,, but to a peritec- 
tic reaction involving formation of crystals of a 
third compound corresponding to composition CoSbs. 


Preparation of the Alloys 


Cobalt and antimony were fused together in 
various ratios from small concentrations of Co to 
50 at.% in quartz ampoules with an argon atmos- 
phere. Both slowly and quickly cooled alloys were 
Studied. 


Structural Investigation 


Alloys containing less than 13.89 wt.% Co con- 
sist of CoSb3 crystals and the eutectic CoSb3—Sb. 
With increasing Co content, the quantity of precipi- 
tating CoSbs crystals rises. The peritectic reac- 
tion for the formation of CoSbs; crystals proceeds 
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very slowly, and in order to pass to completion re- 
quires annealing for some hours at a temperature 
close to 770°C. 

Annealed alloys of composition CoSbg3 (13. 89 
wt.% Co) are homogeneous and have a polyhedral 
structure. Figure 2 shows a micrograph taken at 
magnification 140 with light polyhedral crystals of 
CoSb3; the dark regions are pores. 
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STRUCTURE OF COBALT — ANTIMONY ALLOYS 


Fig. 3 


An x-ray photograph with more than 40 clear 
lines was obtained from an annealed sample of 
CoSb3 by the powder method in a 114.4-mm diam-~- 
eter camera, using CuKg radiation. On calculating 
and indexing the diffraction picture, it proved that 
the crystals belonged to the cubic system. The unit 
cell was body-centered with spacing a = 9.016 + 
0.003 kX and contained eight CoSb3 units, i.e. , 32 
atoms of Co and Sb. Indexing the powder photo- 
graphs and noting the extinctions showed that the 
structure belonged to one of the two x-ray groups 
mI — (N102) or m8mF...(N112) [5]. Rapidly cooled 
alloys of composition CoSh3 were heterogeneous. 
Figure 3 shows a micrograph of such an alloy taken 
at magnification 300; in the etched field of the eu- 
tectic, gray crystals of CoSb, can be seen, sur- 
rounded by a thin layer of CoSbs crystals. Annealed 
alloys of composition CoSb, (19.49 wt.% Co) are 
homogeneous and constitute the compound CoSh,. In 
Bokii's handbook [6], there is an indication that the 
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TABLE 1 
e a Mean micro-| Spread of 
Compound hyd; x”, | hardness, | microhardness, 
g/cm* | g/cm 2 2 
kg/mm kg/mm 
CoSb3 T.2 7.68 459 420—190 
CoSby Tey. = 565 520 —610 
CoSb 8.6 8.93 580 520—660 


compound CoSb, crystallizes in the rhombic system 
and has the structural type of marcasite, but 

no values are given for the unit-cell parameters. On 
rapidly cooling an alloy of composition CoSh,, 
heterogeneous structure occurred including crystals 
of all three compounds (CoSb, CoSb,, and CoSbs). 
Annealed alloys of composition CoSb (32.62 wt.% 
Co) are homogeneous and represent the compound 
CoSb, which has the structural type of nickel arsen- 
ide [6]. 


Density and Microhardness of the 
Compounds 


The density of the compound Opyq was deter- 
mined for bars weighing 0.5 to 2 g by hydrostatic 
suspension [7] in carbon tetrachloride. The table 
shows the densities ©hyq and the densities ox cal- 
culated from the x-ray diffraction data at room tem- 
perature. The x-ray density of CoSb was calculated 
from data given in [6]. 

The microhardness of the alloys was measured 
on the PMT-3 with a constant indentor load of 20 g. 
The table shows the mean microhardness together 
with its scatter. We see from the table that with 
increasing cobalt content the density and microhard- 
ness of the compounds increase. 
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It is shown that the compounds CoSb3, RhSb3, and IrSb3 crystallize in the cubic system 


and are isostructural with CoAs3. 


It was shown in [1] that, in addition to the com- 
pounds CoSb and CoSh,, the cobalt—antimony sys- 
tem contained a compound CoSb3 with a = 9.016 kX. 
The chemical analogs of cobalt, rhodium, and iridi- 
um should also form compounds of composition 
MeSbs with antimony. 

In order to prove this, we fused rhodium and 
iridum together with antimony in the ratio MeSb3 in 
quartz ampoules in an argon atmosphere. The 
alloys of composition RhSb3 and IrSb3 were annealed 
in an argon atmosphere at 700-800°C for 5-8 h and 
cooled in the furnace. These annealed alloys ap- 
peared practically homogeneous in the microscope; 
only in individual small regions did the peritectic 
reaction fail to pass to completion. Powder x-ray 
diffraction photographs were taken of these alloys 
in a 57.3-mm diameter camera in CuKy radiation; 
calculation and indexing of these showed that the 
three cubic compounds CoSb3, RhSb3, and IrSb3 
were isostructural. 

For further x-ray study of the crystal struc- 
ture of CoSb3, RhSb3, and IrSb3, it was essential to 
have fairly pure powders of each phase. In order to 
obtain these we prepared alloys richer in antimony 
than MeSb3. From microscopical data, these al- 
loys consisted of MeSbs and antimony crystals. The 
MeSb; crystals were separated from the alloy by 
dissolving the antimony in aqua regia. 

Using the pure crystalline phases of CoSbs, 
RhSb3, and IrSb3 obtained in this way, x-ray pow- 
der photographs were obtained in a 114.4-mm 
diameter camera in CuKg radiation; these contained 
more than 40 sharp lines. The x-ray lines were in- 
dexed on a bec cell. The accurate lattice spacings 
are shown in Table 1. 


Theinterplanar spacings calculated from these 
data agreed closely with the measured values 
(Table 2). 

The density of the CoSbs3, RhSb3, and IrSbg 
crystals was determined by hydrostatic weighing in 
carbon tetrachloride, and the results are given in 
Table 1. The unit cells of CoSb3, RhSb3, and IrSb3 
contain eight formula units. Inspection of known 
structures of composition AB; and comparison of 
their unit cells with those of CoSb3, RhSb3, and 
IrSbs showed that these latter compounds are iso- 
structural with CoAsg [2,3]. 

The positions of the atoms in the space group 
Im3(T?) for the CoAsg3 type of structure are as fol- 
lows [8] (Fig.1): 

Co 8 (c) 4%, %, 4, ete. (without parameters) 

As 24: (g) Oyz, etc. 

The structure amplitude depends on parameters 
y and z, determining the positions of the Sb atoms. 

The line intensities were calculated from the 
formula 


I ~ Ly: p-F?, 


where Ly is the Lorentz factor and p the recurrence 
factor for parameter values y = 0.35 and z = 0.15. 

Comparison between the calculated and meas- 
ured intensities (Table 2) shows satisfactory gen- 
eral agreement. Occasional discrepancies are 


TABLE 1 
Compound a in kX ° hyd», oes 
g/cm g/cm 
CoSb; 9.016-++0.002 st 7.69 
RhSbs 9.211+0.002 7.9 7.96 
IrSb3 9.230-+0.002 8.9 9.41 
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TABLE 2 


hkl 


110 
200 
211 
220 
310 
222 
321 
400 
411; 430 
420 
332 
422 
510; 481 
521 
440 
530; 433 
600; 442 
611; 532 
620 
541 
622 
631 
444 
710; 550; 543 
640 
121; 633; 552 
642 
730 
732; 651 
800 
554; 741; 811 
820; 644 
653 
822; 660 
831; 750; 743 
662 
152 
840 
833; 910 
842 
921; 761; 655 
664 
930; 851; 754 
932; 763 
844 
941; 853; 770 
00 10; 068 
11 10; 277 
02 10; 268 
059; 349 
Z2 10; 666 
13 10; 259; 567 
178; 558; 477 
04 10; 468 
169; 33 10 
24 10 


01 11; 378; 459 
1211; 15 10; 369 


088 
079; 0311 
288; 44 10 


279; 2311; 35 10; 677 


06 10; 668 
14 11; 578 
26 10 
569 
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6.37 
4,51 
3.68 
3.19 
2.85 
2.60 
2.41 
2,29 
2.12 
2.02 
1.92 
1.84 
AT 
1.65 
1.59 
1.54 
1.50 
1.46 
1.42 
1.39 
1.36 
1.33 
1.30 
1.27 
1.25 
1.227 
1.205 
1.184 
1.145 
1,127 
1.110 
1.093 
1,078 
1.062 
1.048 
1.034 
1,021 
1,008 
0,995 
0.984 
0,972 
0.961 
0.950 
0.930 
0.920 
0.911 
0.902 
0.893 
0.884 
0.876 
0.867 
0.859 
0,844 
0.837 
0.830 
0.823 
0.816 
0.8032 
0.7969 
0.7907 
0.7847 
0.7789 
0.7731 


CoSb, 


3.67 


2,84 


2.40 
2.12 
2.01 
2.92 
1,84 
1.76 


0.876 
0.867 
0.859 
0.845 
0.837 
0.830 
0.833 
0.816 
0.8032 
0.7907 
0.7847 
0.7789 
0.7331 


418 
14,5 
233 
0.4 
2682 
10.6 
647 
2.3 
241 
706 
264 
309 
952 
12 
21 
256 
217 


191 
407 


- CO 


342 

81 
264 
177 
425 
351 
160 


V.W. 


6.51 
4,60 
3.76 
3,25 
2.91 
2,66 
2.46 
2,30 
2.17 
2.06 
1.96 
1.88 
1.80 
1.68 
1.63 
1,58 
1.53 
1.49 
1.45 
1,42 
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0.993 
0,982 
0.971 
0.950 
0.940 
0,981 
0,921 
0.912 
0.903 
0.894 
0.886 
0.878 
0.862 
0.855 
0,848 
0,840 
0.834 
0,820 
0,814 
0,808 
0.8017 
0.7957 
0.7898 
0.7841 
0.7784 
0.7730 


RhSb, 


6.48 
4.59 
3.73 
3.25 
2.90 
2.65 
2.46 
2.16 
2.06 
1.96 
1,87 
1.80 


1.018 
1.003 
0.993 


0.971 
0.949 
0.939 
0.931 
0.920 
0.912 
0.903 
0.894 
0.878 
0,862 
0,855 
0,848 
0.840 
0.834 
0.820 
0.814 
0.808 
0.8015 
0.7954 
0.7897 
0.7840 
0.7784 
0.7730 
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IrSb, 
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probably connected with slight deviations of the an- 
timony parameters from the values taken above. 
The hemihedral space group of skutterudite and its 
analogs is responsible for the appearance of lines 
with the same absolute indices but different signs 
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(and hence different structure factors) on the Debye 
photographs. Hence, the parameters can only be 
refined by using the intensities of reflections from 
single crystals. 

This investigation shows that a whole series of 
isomorphic compounds crystallize in the skutterud- 
ite structure with square arrangements of the type 
(Asal a It is extremely probable that compounds 
RhAsz; and IrAs3 with analogous structures also 
exist. 

We are grateful to Engineer R.N. Kuz'min for 
help in carrying out the experiments. 
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The unit cell and space group of lead diethyldithiocarbamate were determined. A 


general solution is given for the structure. 


position. 


Lead diethyldithiocarbamate 


C,H5 S 
< / 
. N — CK Pb 
2 


CoH; Ss) 5 


SH 


in which the hydrogen atoms of an amino group are 
replaced by organic radicals, C,H; groups. Free 
and thiosubstituted carbamic acids are unstable 


compounds and cannot be separated in the free state; 


their salts, however, are stable compounds. These 
dissolve easily in organic solvents (some of them 
simply in water), melt easily, and enjoy wide prac- 
tical use in accelerating the vulcanization of rubber, 
as insectofungicides, medicinal substances, etc. 
This group of compounds has been little studied as 
regards physical and chemical properties and atom- 
ic structure. 

Cambi and Malatesta [1] studied a number of 
dithiocarbamate compounds. With the help of x-ray 
analysis they established isomorphism between the 
salts of tervalent iron and cobalt. A study of mag- 
netic moments led to certain conclusions regarding 
the nature of the chemical bond in the molecules. 


(RaNCSp)3Co, 
(RgNCSo)sFe, 


=) 
u == 5.5 (Bohr magnetons) 
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The Pb—S bonds have a pyramidal dis- 


In the compound with tervalent cobalt, the co- 
valent type of bond predominates, phyd = 0. Inthe 
second compound the bond was assigned an ionic 
nature, since this agrees with the magnetic mo- 
ment of the Fe’’’ ion in aqueous solution, uw = 5.9 
[2]. 

Later investigations, however, associated the 
magnetic moment in iron complexes with the exist- 
ence of covalent bonds [3]. The isomorphism was 
explained by the fact that the octahedral bonds of 
the metal atoms had similar directions. 

The lead atom is characterized by the following 
occupation of. valence orbits: 


leaabal a 


this allows a different configuration of the bonds in 
the molecules and crystalline compounds, namely, 
a triangular configuration with p* electron orbits, 

as observed in the structure of PbCl,. As a result 


om!) 


® Pb 


Fig, 1. Structure of PbO, 
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TABLE 1 
Spacings, A Density 
Crystal P 6 = hs ARS) 
Formula | Oe SOROS ETON) |S cE 
system | 4 | i | ° | Be SO5| oy 
Ze Aso 
i . | 
{(CoH;)oNCSs]gPb |Mono- | 9.55 |11.75|14,72) 96°] C3,—P&/e | 4 | 2.02) 2.04 
clinic 


of the unpairing of the s* electrons, tetrahedrally 
directed sp® bonds are formed, as established by 
Zhdanov and Ismailzade [4] in the structure of tetra- 
phenyl lead Pb(C,H;), and also in Pb(CHs), [5] and 
Pb» (CH3)¢ [6]. 

Pyramidally directed bonds were also estab- 
lished for Pb, for example in the oxide PbO. Moore 
and Pauling [7] showed that in the PbO structure 
the lead atoms lay at the vertex of a square pyra- 
mid, the base of whichcomprised four oxygen 
atoms, the angles O—Me-—O being 118 and 75° (Fig. 
iy 

The metal atom forms four electron orbits di- 
rected toward the base of the pyramid, the vertex 
of which is occupied by an unshared pair of elec- 
trons. The oxide SnO has the same structure. 

In another form of lead oxide, the lead atomis 
connected with three oxygen atoms in the rhombic 
system. 

Cox [8] considered the dimensions of the unit 
cell of the complex lead oxalate K[Pb(C,O,).], and 
concluded that the four Pb—O bonds had plane con- 
figuration; the atomic structure of this compound 
was not solved, however, and this bond configura- 
tion could not be proved. 

Bystrom and Evers [9] studied the structures 
of Ag,PbO, and Ag;Pb,O0,. The first compound be- 
longs to the monoclinic system, and the position of 
the oxygen atoms relative to the lead is the same 
as in tetragonal PbO. The oxygen atoms form se- 
verely distorted tetrahedra; the lead atoms, how- 
ever, lie outside these. The structureof Ag; Ph)Og 
is tetragonal, and the lead atoms in this are sur- 
rounded octahedrally by oxygen atoms. 

Shirane [10] and colleagues studied the struc- 
ture of PbTiO; by x-ray and neutron diffraction.The 
configuration of the lead and oxygen atoms was the 
same as in tetragonal PbO. The shortest Pb—O 
distance was 2.53 A. 

Compounds of lead with oxygen analogs (sulfur, 
selenium, tellurium) belong to the NaCl-type struc- 
mules (LIL IVA 

Specimens of lead diethyldithiocarbamate ob- 
tained from E.N. Gur'yanova (Molecular-Structure 
Laboratory) were recrystallized in benzene. The 


crystals had the form of hexahedral prisms of di- 
mensions 0.6 x 0.8 x 1 mm. Table 1 gives the x- 
ray data. 

In addition to this, crystals with unit-cell di- 
mensions a = 9.55, b=11.75, c= 44.1 A, contain- 
ing 12 molecules, were found. 

In x-ray rotation photographs taken along the 
[100], the odd layer lines were very severely 
weakened, owing to the disposition of the lead atoms 
in positions x = 44 and 4. The Laue photographs 
reveal pseudo-symmetry planes due to the same 
heavy lead atoms. 

The structure was solved by means of F’ series. 
The initial data included the measured intensities of 
82 Oki and 86 hOl reflections. 

Figure 2 shows the projection of the interatom- 
ic—vector function on the (100) plane, from which 
the configuration of the Pb—S bonds was determined 
The Pb Atom lay at the vertex of a four-sided pyr- 
amid having a base made up of four S atoms. 

The maxima in the (010) projection were poor- 
ly resolved owing to the large superposition num- 
ber of interatomic vectors for a small number of 
terms in the series. The maxima corresponding to 
the S atoms in the (100) projection were unresolved 
as a result of the superposition of two pairs of atoms. 
Sj, Syy, and Sy, Spy. 


Tetrahedral disposition of the S atoms around 
the lead atom was excluded owing to the positioning 


@ 


Fig. 2. Projection of F series on the (100) plane, 
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TABLE 2 
Projection (100) | Projection (010) 
Atoms 
y z x 
Pb Otte 0 0.250 0) 
$5 0.16 0.13 | 0.47 0.43 
4 0.29 | 0.94 | 0.07 0.94 
S45 0.16 0.13 | 0.03 0.12 
ore 0.29 0.94 | 0.42 0.9 


Yuc 
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Fig. 3, Electron-density projection on the (100) plane, 


of the Pb atoms indicated above. Plane configura- 
tion of the Pb—S bonds was not supported by the pro- 
jection on the (010) plane. 

The electron-density projection on the (100) 
plane, constructed with structure-amplitude signs 
determined from the coordinates of the heavy atom 
(lead), excluded triangular configuration of the Pb—S 
bonds and supported pyramidal configuration.A re- 
fined electron-density series was then constructed, 
the signs of the structure amplitudes being cor- 
rected in accordance with the S coordinates obtained 


Fig. 4, Electron-density projection on the (010) plane. 


DIETHYLDITHIOCARBAMATE 409 
TABLE 3 

Atoms a 
Pbh— Pb 4 25 
Phesetie Dee 
Pb—Sry ry ea) 
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from the pyramidal model. A (010) electron-den- 
sity projection was also constructed. 

From these projections (Figs. 3 and 4) the co- 
ordinates of the Pb, S;, Syz, Spyz, and Syy atoms 
were obtained (See Table 2). 

The interatomic distances obtained are given in 
Table 3. 

The (100) electron-density projection contained 
maxima corresponding to some of the light atoms 
and also superfluous maxima due to truncation 
waves. 

The dependence of these superfluous maxima on 
the truncation wave is made clear by constructing 
the theoretical electron-density series, using the 
coordinates of the heavy lead atom. 

The resolving power of the series is reduced 
owing to the absence of reflections 0k/ and h0l with 
k+l =2n+1 and h= 2n+1 (required for space 
group C3, — P2;/c) as a result of the large differ- 
ence between the scattering power of the Pb atoms 
and the other atoms in the structure. 

The maxima on the (010) projection were more 
poorly resolved on account of the superposition of 
atoms of neighboring molecules. 

Geometrical analysis of the structure con- 
firmed the model obtained from the electron-density 
projection. Figure 5 shows the disposition of the 
molecules in the unit cell. 

In this way we found that the structure of lead 
dithiocarbamate contains a pyramidal configuration 
of the lead—sulfur bonds analogous to that earlier 
established for lead oxide. 

A review of crystal chemistry [13] indicated 
the desirability of studying whether this configura- 
tion could occur in finite complexes, since it had 
formerly only been established for layered structures 
of metal oxides with infinite lattices of oxygen atoms. 
This we have done in the present paper, which is 
the first to indicate tetrahedral—pyramid configura- 
tion in a lead complex. 

The interatomic distance Pb— Pb = 4.95 A in- 
dicates the absence of any chemical bond between 
the lead atoms, such as occurs in the compound 
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Fig. 5. Disposition of the molecules in the 
unit cell [(100) plane]. 


Pb )(CH3)g, where we have 2. 88 A [6]. Here this 
distance exceeds the atomic diameter of lead still 
more on account of the way in which the molecules 
are packed in the cell. 

The Pb—S bond has a predominantly covalent 
character, as confirmed by the directional proper- 
ties of the valences and the interatomic distances. 
For the ionic type of bond with coordination number 
6, the Pb—S distance equals 3.08 A eae 

Analysis of the interatomic distances in lead 
dithiocarbamate shows that the Pb—S bonds are re- 
alized with a deficit in the number of electrons; 
hence, on hybridization of the bonds,an unshared 
pair of electrons may be retained in the Pb, since 
on formation of the tetrahedral bonds all the val- 


ence electrons of the Pb atom take partin these bonds. 


Table 4 [14] gives the values of the dipole mo- 
ments for a series of compounds of dithiocarbam- 
ate salts. 

The Pb—S bond has an intermediate sheers 
between covalent and ionic; hence,the Pb 
bond is polar. Thus, the increase in dipole mo- 
ment on passing from zinc to lead is explained by 
the different directional properties of the valences 


TABLE 4 
(CsHyoNSo)oZn | w- 10-28 4.89 
(C5Hi oN Ss)oNi 1.94 
(CsHi9N So)oPb Dace 
( 


) 
C5Hy,NS2)3Bi 4.05 
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in the molecules. As the structural studies show, 
the pyramidally directed Pb—S bonds explain the in- 
crease in the dipole moment. 


Conclusions 


1. Single crystals of lead diethyldithiocarb- 
amate were obtained in the form of colorless hexa- 
hedral prisms. 

2. X-ray diffraction revealed a monoclinic 
iy Bis Ov Tike (isp S148 72 A:B =96° 
The space group 


unit cell: 
with Z = 4 molecules per cell. 
was C8, — P2,/c. 

3. A pyramidal configuration of the Pb—S 
bonds was established by means of an F’ series. 
From electron-density projections on the (100) and 
(010) planes the coordinates of atoms Pb and Sj_ly 
were obtained, together with the interatomic dis- 
tances in the pyramidal complex. 

4, On the basis of these structural data, the 
dipole moments of a number of dithiocarbamates 
are considered. 


LITERATURE CITED 


1. L. Cambi and L. Malatesta, "Magnetismus 
und Polymorphie der inneren Komplexsalze. 
Eisensalze der Ditiocarbamidsauren," Ber. , 
70; 2067 (1937): 

2. L. Pauling, Nature of the Chemical Bond 
[Russian translation] (1947). 

3. R.S. Newholm, Magnetism and Inorganic 
Chemistry [Russian translation] (Advances in 
Chemistry, 25, 4, 517-544, 1956). 

4, G.S. Zhdanov and I.G. Ismailzade, "Crystal 

structure of metalloorganic compounds. I. 

X-ray study of crystals of tetraaryl compounds 

of silicon, tin, and lead," Zh. Fiz. Khim. ,24, 

12, 1495-1501 (1950). am, 

A.F. Wells, Constitution of Inorganic Sub- 

stances [Russian translation] (1949). 

6. H.A. Skinner and L.E. Sutton, "The electron 
diffraction investigation of the molecular struc- 
ture of dilead-hexamethyl," Trans. Faraday 
Soc. ,36, 1209 (1940). 

. W.J. Moore and L. Pauling, "The crystal 
structures of the tetragonal monoxides of lead, 
tin, palladium, and platinum," J.Am.Chem. 
Soc. ,63, 1392 (1941). 

8. E.G. Cox, A.J. Shorter, and W. Wardlaw, 
"Stereochemistry of bivalent tin and lead," 
Nature, 139, 71 (1937). 

9. A. Bystrom and L. Evers, "The crystal struc- 
tures of Ag, PbO, and Ag; Ph.O,," Acta Chem. 
Seand. ,4, 4, 613 (1950). 


on 


- 


10. 


ike 


12. 


13. 


STRUCTURE OF LEAD DIETHYLDITHIOCARBAMATE 


G. Shirane, R. Pepinski, and B.C. Franzer, 
"X-ray and neutron diffraction study of ferro- 
electric PbTiO;," Acta Cryst. .9, 137 (1956). 
G.B. Bokii, Introduction to Crystal Chemistry 
[in Russian] (1954). 

B.F. Ormont, Constitution of Inorganic Sub- 
stances [Russian translation] (1950). 

G.M. Powell, "Inorganic structures," Collec- 


14, 


411 


tion: New Studies in Crystallography and 
Crystal Chemistry [Russian translation](1951), 
Vols 3, ). cue 

L.S. Kuzina. Dissertation: Isotopic Exchange 
of Sulfur and Certain Derivatives of Dithio- 
carbamic Acids and Their Structure [in Russian] 
(L. Ya. Karpov Physicochemical Institute). 


SOVIET PHYSICS — CRYSTALLOGRAPHY 


VOLY TNO? S SEPT .-OCT., 1956 


ISOMORPHISM OF THE FERRATES AND TITANATES 
OF STRONTIUM, BARIUM, AND LEAD 
E. G. Fesenko and O. I. Prokopalo 


Scientific-Research Physicomathematical Institute at the 


V.M. Molotov State University, Rostov 


Translated from Kristallografiya, Vol. 1, No. 5, 


pp. 520-523, September-October, 1956 
Original article submitted July 14, 1956 


Lead ferrate was obtained and its electrical properties studied, together with those of 


barium and strontium ferrates. 


The remarkable dielectric properties of com- 
pounds with the perovskite-type structure stimulate 
a search for new compounds with this type of struc- 
ture. A recent communication [1] described the 
synthesis of barium ferrate in two forms, tetragon- 
al and hexagonal, and also the production of stronti- 
um ferrate [2]. The x-ray structural work de- 
scribed in these papers showed that the parameters 
of the compounds in question were extremely close 
to those of the corresponding titanates. 

The present investigation was undertaken in 
order to synthesize lead ferrate (this has not hither- 
to been done, so far as we know) and study its elec- 
trical properties together with those of barium and 
strontium ferrates. 

As original components for the synthesis we 
took carbonates of Sr, Ba, and Pb, and also Fe,O3 
(all of analytical grade). The latter was subjected 
to chemical analysis for iron content and the data of 
the analysis were taken into account in calculating 
the components. Roasting took place in an oxygen 
atmosphere. A sample of stoichiometric composi- 
tion SrFeO3 was obtained on roasting for 1 h at 
1200°C; BaFeOs; was roasted for 4 h at 950°C, and 
PbFeOs for 1 h at 1000°C. 

The x-ray photographs were taken by the pow- 
der method in an 86-mm diameter x-ray camera us- 
ing Cu and Fe Kq radiation. The photographs of all 
three substances contained very weak lines of Fe,O3 
and the oxides of Sr, Ba, and Pb, respectively. 
This indicated that the conditions specified had 
taken the reaction to a fair conclusion, although the 
samples obtained were not analyzed chemically. 

In all our attempts at preparing barium ferrate 
the cubic form was always obtained. We were un- 
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able to produce the hexagonal and tetragonal modi- 
fications, in spite of observing the conditions specified 
in [1]. We did establish a variation of the lattice 
spacing of cubic barium ferrate with roasting tem- 
perature. The spacing a varied linearly from 

4.007 A (950°C) to 4.010 A (1150°C). The cubic 
form of BaTiO; showed a variation of the opposite 
kind. As we know [8], the cubic form of barium 
titanate can be obtained from Ba(OH), - 8H,O and 
TiO, at 400°C. We found that a sample obtained at 
this temperature had a = 4.014 A and one obtained 
at 800°C had a = 4.010 A. The different effects of 
the roasting temperature may be explained by the 
fact that, in barium ferrate, there is a partial trans- 
formation of the quadrivalent iron into tervalent on 
heating above the formation temperature, corre- 
sponding to a weakening of the interatomic forces. 

In the case of BaTiO;, heating above 400°C leads to 
ordering of the structure. 

The parameters obtained for our substances are 
given in Table 1, which,for comparison ,also con- 
tains data for the forms of barium ferrate and titan- 
ate given in [1] and [4]. 

The data in the table illustrate the complete 
isomorphism of the ferrates and titanates and show 
that the cell volumes of the ferrates are smaller 
than those of the corresponding titanates. At the 
same time we notice the extremely small change in 
spontaneous deformations on replacement of the 
octahedral ion. 

This similarity is based on the common nature 
of the chemical bond in titanates and ferrates; in the 
BO3 groups, where B= Fe or Ti, the bond in both 
cases is effected by 4s and 3d electrons. The dif- 
ference observed is due to the fact that the state of 
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TABLE 1 
, Stoichiometric composition 
Modification Ae 
(A SrTiO, | SrFeO, | BaTiO, | BaFeO, | PbTiO, | Pbreo, 
Cubic al 3.905 3,869 4.010 4 O07 — a 
Tetragonal a ~ 3.994-] 3.98 | 3.900 | 3.887 
c — — 4.034 4.01 4.151 4.123 
c/a 1.010 1.01 1.064 1.060 
Hexagonal a _- _ 5.74 5.68 = = 
c — 14.05 13.86 = os 


the Ti ion is energetically more favorable in estab- 
lishing the bond. 

Comparison of the lattice parameters of the 
titanates and ferrates enables us to estimate the 
radius of the Fe** ion, regarding which data ob- 
tained by interpolation with respect to the radii of 
Mn and V ions exist in [2]. The radius of the Fe** 
ion may also be estimated from the volume of the 
octahedral cavity which it occupies. By comparing 
the lattice spacings of BaFeO3 and BaTiO; [1], it 
was earlier suggested that the radius of Fe** was 
smaller than that of Ti‘. This selection of materi- 
al, however, was unfortunate, since, for barium 
titanate Gol'dshmidt's geometrical criterion t ex- 
ceeds unity, and hence the dimensions of the cen- 
tral ion can only indirectly influence the cell size. 
It is better to estimate the radius of Fe!” from 
strontium titanate and ferrate, for which t < 4. 

The lattice spacings of the cubic forms found 
experimentally and calculated from the ionic radii 
are, respectively, 4.010 and 4.10 A for Ba TiO3 
and 3.905 and 4.00 A for SrTiO; [5]. The theoreti- 
cal overestimates indicate the necessity of correct- 
ing the tabulated ionic radii used. This should be 
done primarily for oxygen and titanium ions. In 
fact, since the nature of the bond in the BOs groups 
is partly covalent [6], the oxygen —ion radius should 
be taken rather smaller and the titanium —ion radi- 
us rather larger than those indicated above. The 
same conclusion may easily be drawn from the fol- 
lowing. It is known that the octahedral cavity in 
BaTiO; and the cubo-octahedral cavity in SrTiO; 
are somewhat free, respectively, for Ti and Sr 
ions. In the first of these the parameter is deter- 
mined by a framework of Ba and O ions and in the 
second by Ti and O. Since the experimental pa- 
rameters in both BaTiO; and SrTiO; are smaller, 
the radius of the oxygen ion must be taken rather 
smaller than that indicated. 

From the experimental value of a for cubic 
BaTiO,, regarding the bond between the Ba and 
other ions as covalent, we can estimate the aver- 
age radius of oxygen as approximately 1.28 A. The 


octahedral cavity occupied by titanium in SrTiO 
can contain an inscribed sphere of radius 0.67 A. 

The isomorphism of SrTiO3 and SrFeO3 en- 
ables us (by taking the oxygen radius found) to de- 
termine the radius of the sphere corresponding to 
the octahedral cavity in SrFeO3. This proves to be 
0.65 A. If we consider the radius of the Fe** ion 
to be as much below the radius of the sphere in the 
octahedral cavity of SrFeO3 as the radius of the 
Ti** ion is in SrTiO;, we may estimate the radius 
of the Fe** ion as about 0.62 A. The Fe" radius 
may possibly be rather smaller, since we took the 
averaged parameter value for SrFeO3, the lattice 
of which is probably defective. 

In view of the isomorphism of the titanates and 
ferrates, we may expect that tetragonal BaFeOs and 
PbFeO3 have phase transformations at temperatures 
higher than the corresponding titanates. In fact, 
bearing in mind the instability of quadrivalent iron, 
we may assume the existence of defects in the 
BaFeO3 lattice. A result of this will be the weaken- 
ing of the internal field. With increasing concentra- 
tion of defects the spontaneous deformation will fall 
to zero. The cubic modification corresponds to 
this limiting case. Considering the similarity of the 
spontaneous deformations in the titanates and fer- 
rates, and also the existence of defects in the lat- 
ter, we come to the conclusion that the ferrates 
have a high temperature of the phase transformation. 

The ferrate samples obtained had a high elec- 
trical conductivity. Figures 1-3 show Ino = f(1/T) 
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curves for our ferrates of strontium, lead, and 
barium, respectively. The activation energies or 
(allowing for the large electrical conductivity) ener- 
gy gaps between the normal band or local levels and 
the conduction band are calculated from the formula 
w = [nog —Ino,)/ 1/ T,; —1/T,)]k. For strontium 
ferrate, w= 0.10 eV. In the temperature range 20 
to 290°C, there are no changes in the mechanism of 
electrical conduction. The In c= f(1/T) curves for 
lead and barium ferrate have a break at Ty equal to 
180 and 130°C, respectively. The value of w at the 
break point changes from 0.50 to 0.63 eV for 
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PbFeO; and from 0.34 to 0.60 eV for BaFeO3. The 
existence of the break may be explained by a change- 
over from the normal band to the conduction band, 


while the electrical conduction at temperatures be- 


low Ty is mainly determined by the distance from 
the conduction band to the local levels arising as a 
result of impurities. It is possible, however, that 
the break in the Ino = f(1/T) curves is connected 
with structural changes. 

At the present time, the conditions for synthe- 
sizing ferrates are being optimized with a view to 
studying their properties further, with the closest 
possible approach to the composition ABO3. 
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Diphenyliodonium chloride and iodide are isomorphic crystals: (CgHs),ICl, a = 20.81, 
bisa /e2 1620.26 A, Pf > 102.5°) m= 8. (CpHa)] <1, a =22208, b1600 7 ,uer= 

20.42 1G B =101°, n= 8. The space group is C2/c. These were studied by x-ray 
diffraction using Mo Ka radiation. In view of the isostructural connection between the 
two substances, only (CgH;),IC1 was studied completely. The coordinates of the atoms 
were determined by calculating a three-dimensional electron-density series containing 


more than 2000 independent Fy] terms. 


The accuracy in determining the positions of 


the light atoms (C) was +0.03 to 0. 04 A; that of determining the heavy atoms (I, Cl) was 


SAO VE 


Introduction 


The chemical nature of iodonium compounds 
manifested in exchange and decomposition reactions 
is of considerable interest. Thus, for diphenyl- 
iodonium chloride and iodide, as well as for some 
of their derivatives substituted in the benzene rings, 
there are some typical reactions [1,2] involving a 
stage in which free radicals are formed (I) and 
others involving heterolytic decomposition of the 
molecule (II): 


I C,H4I#C! + Hg 
GyHi 
cl 


I [ (CeHs)P + Cots! oT: eq gt) se ((CgH3),P]* CI+CgHs! (a) 
t — 


CH, CH; CH, CH, CH, 
+ 
+ 
I is I sd I 
—— i +1> —2 (b) 


Certain physical properties and physicochemi- 
cal data [3] indicate an ionic character for the 
halide —halide bond in diphenyliodonium halides. 
Apart from studying the character of this bond, it 
is also interesting to determine the configuration 
of the molecule (valence configuration of tervalent 
iodine, valence angles C—I-C and C—I-X). 


CoH ,HgC1+C,Hl 
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The problem of the present x-ray diffraction 
study lay in establishing the configuration of the 
molecule in the crystal and elucidating the connec- 
tion between its structure and chemical properties. 


Experimental Part 


Crystals of diphenyliodonium chloride and iod- 
ide ,(CgH;),.I—Cl and (C,Hs).I—I,were obtained by 
slow evaporation of the saturated solutions in ethyl 
alcohol. The crystals had the form of platelets 
[colorless for (C.gHs5).I Cl and yellowish for (C,Hs)2I— 
I] strongly drawn out in the direction of the b axis 
and reaching several mm in length (transverse di- 
mensions of the order of 0.2 x 0.4 mm). The 
simple forms were a {100} and c{001}. The unit- 
cell parameters were determined from oscillation 
diffraction photographs in the reciprocal-lattice 
camera. 

The systematic extinctions of hkl reflections 
with h+k # 2n and h0l with h ~ 2n andl # 2n lead 
to two space groups,C2/c and C/c. The fact that 
the first of these relates to close packing, and that 
the crystals have no piezoelectric effect, however, 
leads us to choose C2/c. 


The great similarity between the ac electron- 
density projections of the two compounds and the 
similarity between their unit-cell parameters indi- 
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(Cells) T—Cl | (Go¥s)21 —I 
“ 20.81-+-0.06 A | 22.08--0.10 A 
b 5.82+0.02 A | 6.27-40,03 A 
@ 20. 26-+0.06 A | 20.42--0.08 A 
6 Or era LOT => 
Pmeas| 1-67 g/cm aa 
ex 1,76 g/em® | 4.95 g/cm 
n 5 s 


cates that they are entirely isomorphic. Hence, 

only the diphenyliodonium chloride was subjected to 
a full structural study. The structure of the iodide 
was analyzed on the principle of complete isomorph- 
ism with the chloride. This is justified in view of 
the fact that the x and z coordinates of the heavy 
atoms practically coincide in the two crystals. 

The crystals were photographed in Mo Ka@ ra- 
diation, using the reciprocal-lattice camera. The 
whole experimentally attainable set of reflections 
was recorded; by rotation around the b axis, the 
hol, hil, h2/, h3l, h4l, hdl developments were 
obtained. The total number of reflections recorded 
was around 2000. In order to reduce the measured 
intensities for various layer lines to a single rela- 
tive scale, a series of oscillation photographs 
around the b axis was used. 

Forconverting intensities Inkl into structure 
factors Fed , the Lorentz and polarization factors 
were taken into account. 


Solutions of the Structure 


The x and z coordinates of the heavy atoms (I 
and Cl) in the two substances were found from two- 
dimensional series of interatomic vectors P(x0z) 
and used in determining the signs of Fpo] for sub- 
sequent construction of electron-density series p 
(x0z). The electron-density projection on the ac 
face for diphenyliodonium chloride is shown in Fig1l. 


Fig. 1, Diphenyliodonium chloride, Electron-density 
projection p(x0z) for the independent part of the unit 
cell. 
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This gives a clear picture of the benzene-ring 
maxima as well as those corresponding to I and Cl. 

The y coordinates of the iodine and chlorine 
atoms (in diphenyliodonium chloride only) were de- 
termined by means of a three-dimensional series of 
interatomic vectors (Fhkj series). Attempts to de- 
termine the y coordinates by means of two-dimen- 
sional series of interatomic vectors P (0yz) and P 
(xy0) and to calculate the one-dimensional section 
A (02y'4) of the three-dimensional series by the 
Harker method gave no positive results. 

The three-dimensional series of interatomic 
vectors was calculated by the formula 


wae 
y a > i (Fini + Fini) cos 2xhz cos 2alz 


h=o0 k=0 |[=—oo 


VAGHIS) = 


— (Fixit — Frxi) sin 2xhz sin 2xlz] cos 2ky. 


The symmetry of the series was C2/m and the 
independent part a/4, b/2, c. 

In solving the structure we were only interested 
in aaxima corresponding to vectors between heavy 
atums, a list of which is given in Table 1. 

For the C3, = C2/c group, the FAky series does 
not enable us to find the positions of the glide planes 
c and n directly. The correctness of the positions 
chosen (y = 0 for c and y = y/, for n) was confirmed 
by calculating the I-I, I-Cl, and Cl—Cl distances 
and ata later stage the C —C distances in the 
"broken" benzene ring (Fig.1). 

Starting from the x and z coordinates of the I 
and Cl atoms determined from the electron-density 
projection p (x0z) (these values were close to those 
found from the F},7 series) and the y coordinates 
taken from the Fed series, we calculated the signs 
of all the structure amplitudes, making the justified 
assumption that in the overwhelming majority of 
cases the sign of the structure amplitude was de- 
termined by the heavy atoms. The signs obtained 
were used to construct a three-dimensional elec- 
tron-density series: 


4 ee 
0 (yz) = + >» s > A; cos 2ahy 4+- B, sin 2eky, 
h=0 k=0 1=—oo 
where 


Ay = (Prt + Frat) cos 2xhvr cos 2rlz 


— (Faxi —F xi) sin 2ehex sin 2lz = 0 for 1 + Qn, 


By = (Pag + Faz) sin 2rhx cos 2nlz 


+ (Part — Prri) cos 2ahx sin 2ilz = 0 forl = 2n. 
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TABLE 1. Maxima of Three-Dimensional 
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Interatomic-Vector Series 


Coordinates of | - |e, 
c : Kw» /G 
maximum in |@ 2 |x ¢ 
= relative units 5 5 Os 
E ise ae hopes ane Theoretical coordinates of maxima 
Pera r a 76 
Bo seh Ty Na task 
35 ob ol zea | 8s 
= oO o 2 Su 
a> GadjHo 
: 
| 
1-10 10.0940, 5001 745 | 422 ) es 
I ( | ) 140 (37 MES is 
eel Ey 1 4 | 
= 46|0.500)0,625} 674 | 422 a =a 
I --I |0.146)0.500)0,625} 674 4ay Sih ee a ite 
1--1|0.44610.406|0.125] 316 | 22} ee a ST 
T-—Cl}0.148)0 $0,023) 256 | 42, 20 (2, — %e1) — GT = Yo — GE — Ze) 
ye ie j ay opel , | 
T= CHO Wei 001013921 246) ser | dor Ol” oy YT Ser hs ec 
=e 4 : Pago > 4 ; - 7 ! | 
I --CHO0.14810 0.544) 161 N22 Goi 27s VT + Yep Zien 2) pee 
| | = 
| 


The independent part of the series is a/4, c/2, b. 

Since the a and c spacings of the cell are large, 
the corresponding axes were divided into 120 parts 
for calculating p (xyz), while the b edge was divided 
into 48 parts. The electron-density values were 
not calculated for the whole independent region, but 
only in sections ("prickings") parallel to the b axis 
passing through the regions of electron-density 
maxima which appeared on the ac projection. These 
sections contained all 14 maxima corresponding to 
the carbon, iodine, and chlorine atoms. Theexact 
coordinates of the maxima were determined by in- 
terpolation, using the method described earlier [4]. 
The coordinates of the maxima (in relative units 
and A) for the atoms of one molecule are given in 
Table 2, together with the heights of the maxima. 
The heights of the carbon-atom maxima vary be- 
tween 2.1 and 3.4 (in arbitrary units); an exception 
is that belonging to CJ, whichis 1.3. This is a 
consequence of an incorrect choice of signs for 
some of the structure amplitudes and also the exist- 
ence of a "dead zone," i.e., the omission of cer- 
tain reflections which the method of photography 
used failed to record. 

Table 2 also shows the "ideal" coordinates of 
the carbon and hydrogen atoms' calculated on the 
assumption that the benzene rings were regular 
plane hexagons with side 1.40 A and C—H bonds di- 
rected aiong the diagonal of the benzene ring 1.08 
A in length. The positions of the hydrogen atoms 
(calculated on the assumptions indicated) are re- 
quired in order to study steric problems in the 
molecule and analyze the stacking. All the dis- 
tances between atoms not connected by valences 


were calculated on the basis of the "ideal" posi- 
tions of the carbon and hydrogen atoms. 

The benzene rings observed experimentally are 
practically plane hexagons with the following bond 
lengths: 


I I 


1.39 1,39 1.43 1.40 
1,37 1.43 1,38 1.40 
1,40 1.38 1,36 1.45 


The accuracy of the heavy atom coordinates 
found is 0.01 A. The accuracy for the carbonatoms 
(determined by the deviation of the benzene rings 
from the coplanar state) was 0.04 A; hence, the 
differences in the C—C distances are not subject to 
analysis. 

"Tdeal" 


Calculation of the Carbon 


Atoms 


The scatter in the C—C bond lengths of the 
benzene rings was due to experimental error (main- 
ly incorrect choice of the signs of certain structure 
amplitudes) since the symmetry of the benzene 
rings with respect to the lines C;—C, and C/—C{ 
was slightly disturbed. In order to determine the 
signs of the structure amplitudes for constructing 
the next approximation to the electron-density 
series, slight corrections should be introduced in- 
to the experimental values of the carbon-atom co- 
ordinates so as to transform the benzene rings into 
regular plane hexagons with sides of 1.40 A. “This 


1See below for the method of calculating "ideal" C and H coordi- 
nates. 
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TABLE 2. Coordinates of the Atoms of One Molecule in the Structure of Diphenyliodonium Chloride 


Actual molecule Height of Ideal molecule 
A Relati i aa | A | maxima fn Relative units A 
tom elative units : = 3) arbitrary 
x y Zz x y Zz units x | y Zz x ve g 
Te + + ~— + 
Cy 0.183 |—-0,187 | 0.356 8,80 | —1.09 | 7.21 3.4 0.183 | —0,187 | 0.356 3,80 | —1.90 |: 7,21 
G 0.231 |—0,357 | 0.370 4,80 | —2.08 | 7,50 3,3 0.2380 | —0.359 | 0.369 4,80 | =2,09 | 7,49 
C3 0,232 |—0,521 | 0,318 4,83 | —3,03 | 6.44 2.1 0,238 | —0,521 0,319 4,86 | —3.03 | 6.47 
Cy 0,187 |—0.517 | 0,257 3,89.) —3,01 5,20 2.5 0,187 | —0,516 | 0.258 8.89 | —3,00 | 5,22 
Cs 0.144 |—0,330 | 0,243 2.99 | =-1,92 4.93 2.5 0,144 | —0.328 | 0,243 2.99 | —1.91 | 4,93 
Cs, 0.189 }—-0,176 | 0,294 2589) |\=1 02 | 5,95 2.8 0,188 | —0,168 | 0.294 2.88 | —0.98 | 5.95 
Center of gravity 
of I ring 3.87 | —8,80 | 6,21 
Hp 0,265 | —0.368 | 0,419 5,52 | —2,14 | 8.48 
Hg 0.270 | —0.655 | 0.329 | 5.62 | —3.81 | 6.67 
Hy 0.188 | —0,646 | 0.220 | 3,91 | —3.76 | 4,46 
Hs 0.111 | -0,314 | 0.195 | 2,31 | —1.83 | 3,94 
He 0,102 | —0.031 0,284 2.12 | ~0,18 | 5.75 
ex 0.077 | 0,058 |0.420 | 1.61 0.34 | 8,50 2.1 0.079 | 0,067 | 0,421 1,64 | 0,389 | 8,54 
C4 0.039 | 0.238 |0.386 | 0.82 | 1.39] 7,81 1,3 0,037 | 0,242 | 0.389 | 0.78 | 1.41 | 7.88 
C's —0,029 0,214 0.370 —0.61 jibe) 7,49 Be —0,030 0,207 0,369 —0,64 1.20 7.48 
eh —0,058 |—0,003 | 0.388 |—1.20 | —0.02 | 7.87 2.3 —0,058 | —0.003 | 0.887 |—1,20 | —0,02 | 7.84 
(Gh —0.018 |—0.174 | 0.420 |—0.37 | —1.01 | 8.50 2.1 —0,016 | —0.178 | 0,419 | —0.34 | —1,04 | 8.49 
C'>5 0,050 |—0.151 | 0.437 1,04 | —0.88 | 8.86 2.8 0,052 | —0.143 | 0,436 1.08 | —0.83 | 8,84 
Center of gravity 
of I’ ring 0,22 0.19 | 8.18 
Hy 0,058 0.407 | 0.378 4.21 2.37 | 7.65 
Hl, — 0,063 0,342 | 0.343 | —1.31 1.99 | 6.94 
H', =O Ti |e ed 0,374 | —2.30 | —0,18 | 7.58 
H's —0,037 | —0,340 | 0.431 | —-0,77 | —1.98 | 8.74 
H'g 0.084 | —0,278 | 0.462 1,75 |) —1,62 | 9.35 
I 0.1802 0.0430 | 0.4344 3,75 0.25 | 8,80 53 
Cl 0.3316 0.0636 | 0.4585 6,90 Ose87 | 9.29 13.3 
is achieved by analyzing the experimental coordi- 1.0.6452 + 0.584y — 0.493z — 1.994 = U, 
nates as described below. I’,.0.087z — 0.359y — 0.9292 — 7.614 = 0. 
1. First of all, the method of least squares 
was used to find the equations of the "best" planes The deviations of the carbon atoms from the 
for the two benzene rings corresponding planes do not exceed 0.04 A; this 
quantity characterizes the accuracy of determining 
6 the carbon-atom coordinates. 
1 5 te 2 ‘ : A 
>) (At + By + Cz + DP min, 2. The benzene rings were subjected to an af- 
td fine transformation into the coordinate plane x'0z'. 
whence The transformation requires transfer of the origin 
of the coordinate system to the center of one of the 
6 6 6 6 : A 
2 < a i 
y y afeLaR y Pacer y RACGeNs y sea: toms (in our case, C; and C/), after which the 
ra = ea = formulas for transforming from the "old" system 
6 6 6 6 of coordinates to the "new" and conversely, take 
Gl es 2 "al aml 
A >» vy + BY vir CS Yizi + DY yi =9, the form 
i=1 i==1 t=] i= 
6 6 6 6 1) vw =ayx + ayy + a3 12 
, Oe ] a. N) 2? a. sae feds 
A > mitt B 2 Hits > ai D py % = 0, 2) Y! = yk + yoy + Agoz (a) 
jel i=1 i=1 i=] 3) 2! = Aygt + Aggy + A352 
fy 6 6 6 
al Sy jis tle Ie >) yi tC » Zit yes =), L = Oy;2' + Ayoy’ += A192" 
i=t i=1 i=1 i=1 Y = Ag, L' + Aggy’ + Ayg2' ¢ (bd) 
2 j Z == llyi@ tb Ago’ l= Aone 
where xjyjzj are the experimental coordinates of . ar % 
the carbon atoms. Here, x', y', z' are the coordinates in the new sys- 
Setting D equal to an arbitrary constant, say tem and x, y, z those in the old. 
D=1, we obtain four equations with three unknowns A, Let us denote hj — hy = y', where h;, and hjare 
B, andC.On determining these we find the equations the distances of the first and i-th atoms from the 
for the "best" planes for the Iand I' rings. In the "best" plane, and let us write equations (a) for 


orthogonal coordinate system a*be: atoms C3, Cy, and C; (the most distant from Cj): 
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Ys = Qo%3 + AooY3 + 43223, 
Ya = Qy2%q TT AgeYq + Ag024, 


Ys = Qyo%s + AnoY5 + A305. 


Hence,we find aj, a, and agp. 

Equations (b) 1,2,3 for atom Cy, (xf = 0, y{,z{ = 
Zo, wherez =v [ dy, — (hy — hy)*], and dy is the dis- 
tance between atoms C, and C,) will be 


Ly = AyYa 4+- 320, 
Ya = A244 F A320, 


Zq = AzoYa + A332. 


These give coefficients a3, a3, and a33. 

The coefficient a; is determined from expres- 
sion a,= +V (1— qd, — a’.), the sign being selected 
arbitrarily. Coefficients a, and a3; are foundfrom 
the relations 


1491 + A229 + A303 = O, 
1431 + Ay232 + Ay3A33 = 0. 


3. The next stage was to find the coordinates 
of the carbon atoms in the ideal benzene ring, this 
being so arranged in the plane as to minimize the 
mean-square deviation of its atoms from the ex- 
perimental positions. For this purpose the center 
of a regular hexagon of side 1.40 A was made toco- 
incide with the center of gravity of the experiment- 
al benzene ring, after which the origin of coordi- 
nates was moved to the center of gravity of the lat- 
ter in such a way that the 0Z" axis should coincide 
with the diagonal C;—Cy. 

In the new system of coordinates, x"y"z",the 
regular hexagon was rotated around the 0Y" axis so 
as to make the deviations of its carbon atoms from 
the experimental positions in this system a mini- 
mum. This kind of rotation, as we know, corre- 
sponds to the formula 


expr LiqcOSe + Zjqsing, 

ex — Ziq Sing + Zjqcos¢, 

where sie i Zap ore the experimental coordinates 

and x;q, Zjq the ideal coordinates of the vertices of 
the regular hexagon. 
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The required rotation may be found by minimi- 
zing the square error in the form 


S (a exp “id ©O8 ? — Ziq Sin wo) é 
— min 


ya (Zexpr Tig Sing — Zig Cos 9)? 

Having found the position of the regular hexa- 
gon corresponding as closely as possible to experi- 
ment, we must now retransform to the original co- 
ordinate system. This will give the coordinates of 
the vertices of a regular hexagon, i.e., the ideal 
positions of the carbon atoms in the original system. 

The fully solved structures will be compared 
with the structures of other compounds of tervalent 
iodine in a following article. 

The present paper is the first of a series on 
the structure of onium compounds originating in 
the X-Ray Analysis Laboratory of the Institute of 
Heteroorganic Compounds, Academy of Sciences of 
the USSR. In addition to diphenyliodonium boro- 
fluoride, [(CgHs),IJ° [BF4]~ (which has an ionic 
structure), the structures of certain bromonium 
and chloronium compounds were studied. 

We consider it our duty to thank Yu. T. Struch- 
kov for direct participation in the work and to 
acknowledge the constant interest and attention of 
A.I. Kitaigorodskii, the Director of the Labora- 
tory. 
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The lattice parameters of leucosapphire were determined by x-ray diffraction in the 
temperature range of 20-1000°C; the lattice parameters of ruby were determined as 


a function of Cr,O3 content. 


The important technical applications of pure 
a-Al,O3 crystals, and especially ruby (i.e. , cor- 
undum with chromium impurities) give rise to con- 
tinued interest in their physical properties. 

Despite various approaches made toward the 
study of corundum crystals, no precise values of 
their lattice parameters as a function of tem- 
perature and impurity content have yet been ob- 
tained [1-4]. A recent paper by Tillo [5], more- 
over, contains the doubtful assertion that the intro- 
duction ofchromium into a-Al,O3 only begins to af- 
fect the lattice parameters of corundum after the 


chromium content reaches a critical value of 8 at.%. 


We undertook the work of determining the most 
important lattice constants of pure @-A1,O3 as ex- 
actly as possible, studying their variation with 
chromium content, and finding the thermal-expan- 
sion coefficients of pure corundum and ruby. The 
single-crystal ruby samples used (kindly provided 
for our measurements by 8S. V. Grum-Grzhimailo) 
were described in [6] together with data on their 
chemical analysis. Polycrystalline samples are 
described in [1,2] and their chemical analysis is 
given in [7]. 

The lattice constants of corundum were calcu- 
lated from photographs obtained in the RKU-114 
camera. Because of some uncertainty in indexing 
the diffraction lines, some control measurements 
were made, using oscillation photographs of a 
single crystal in the RKV-86 camera with an at- 


tachment for taking epigrams. This made it easy to 
index the RKU-114 photographs. The constant a 
was calculated from the 330 doublet and c from the 
21.14 doublet. The results in kX are given in Table 
1 together with those of other authors. The con- 
stants (a. and c) for ruby were obtained in a similar 
way and appear in Table 2. The data of Tillo, which 
have a special part to play, are compared with our 
data in Fig.1. 

The monotonic form of our graph suggests that 
Tillo's results [5] for small chromium content are 
wrong. The curious shape of Tillo's curve is prob- 
ably due to the low sensitivity of his apparatus. 

The thermal-expansion coefficients of pure 
corundum and rubies were measured in a"Unicam" 
camera over the temperature range 20-1000°C. The 
stability of the temperature over the sample was 
0.4 deg/mm for a sample length of 15 mm. This 
enabled the lattice constants to be measured with 
an accuracy of 0.001 kX. The results of the meas- 
urements are given in Table 3. 

The linear-expansion coefficients 8; and By 
were calculated from the data of Table 3 by using 
empirical formulas for the experimental relation- 
ships a(t) and c(t). 

The formulas have the form 


a(t) = 4.7474 + 2.8-107 t + 9.7-1079 22, (1) 


c(t) = 12.9554 ++ 8.01055 £ 4- 2.3.-10-8 72 (2) 
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TABLE 1* 
a a 
No, | Year Authors a, kx @, kX 
1 1920 |Davey and Hofman 4.86 13.23 
2 4922 Brag 4.75 12.974 
e. 4927 Harrington 4.761 12.941 
4 1928 |Zachariasen 4.76 13.00 
5 1930 |Passerini 744 12.956 
6 1953 |Swanson and Fuyat 4,749 12.965 at 26°¢ 
ii 1956 |Present work 4,748 PALdo at. 2076 
*Data 1-6 in this table are taken from [4]. 
TABLE 2 
% CreO5 a, kx Ci, TAG 
() 4.7477+4 12.957-+1 
1.5 [9] 4.749044 | 12,9621 
2.9 » 4.752145 12.967+1 
9.6 [12] | 4.764345 | 12.9942 
23.2 » 4.791446 | 13.0712 
48.7 » 4.836546 | 13.1984-2 
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Fig. 2, Lattice spacing a asa 
function of temperature, The 
curve follows formula (1): the 
points give our experimental 
data, 
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Fig. 1. Lattice spacing a asa 
function of CrgO3, wt.%; points 
© represent Tillo's data [5]; @ 
represent our own. 
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Fig. 3, Lattice spacing c asa 
function of temperature, The 
curve follows formula (2); the 
points give our experimental 
data. 
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TABLE 3 
pet a, kX Gs AX 
a-AlyO3 
20 4.748 +1 12,957-+3 
320 4,757+1 12.982+3 
0) 4.762-F 1 12,995+3 
610 4.763-+ | 13.0114 
770 4.775+1 13.0284 
980) 4.7841 13 056-4 
a-Al,Og +9,6%  C€rO3 [b2] 
- 
20 4.765+1 12.9943 
570 4. 784+ 1 13 ,062+-4 
TABLE 4 
Coptiieionss| INO! LS aie 
| 


B, (500°C) | 
8, (500°C) | 7. 


The agreement. between the formulas and the 
experimental results is indicated in Figs. 2 and 3. 
By using the formulas we obtain the following line- 
ar-expansion coefficients: 


B, = 5.85-10°8 + 2.05-10°° 2, 
8, = 6:23:10 — 7 840 fF, 


Measurements of a(t) and c(t) for ruby contain- 
ing 9.6% chromium showed that the linear-expan- 
sion coefficients were larger than those of pure 
corundum. Table 4 compares the average values of 
B14 and f|| for the temperature range 20 to 600°C 
for corundum and ruby. 

Comparison of the linear-expansion coefficients 
found for pure corundum with the results of other 
authors [8-11] shows that our values agree closely 
with dilatometric measurements [8-10] but disagree 
violently with Sharma [11]. 

We should remember that the linear-expansion 
coefficient 6 I depends much more strongly on the 
presence of chromium in @-Al,O3 than does 81. 
Hence, we may conclude that the isomorphic intro- 
duction of Cr in place of Al in the corundum lattice 
primarily tends to weaken the ionic bonds in the 
crystal in the c direction. 

The authors are extremely grateful to S. V. 
Grum-Grzhimailo for drawing their attention to the 
problems considered in this paper and for useful 
comments in discussing the results. 


N. A. SHAL'NIKOVA AND I. A. YAKOVLEV 
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Amorphous films of GaSe prepared by vacuum deposition on celluloid substrates were 
studied. An experimental curve of intensities normalized by the atomic-factor curve 
was obtained. In view of the similarity between the atomic numbers, it was assumed 


that the substance consisted of atoms of one type. 


The radial distribution curve was 


calculated. Interatomic distances r; = 2.40 and rg = 3.75 A were found; the numbers 


of nearest atoms were 6 and 12. 


These data show that in the first coordination groups 


of amorphous GaSe there is no correspondence with the lattice of the crystal phase. 


Binary and ternary systems containing the ele- 
ments Tl, Sb, Se, Ga, As, Bi, etc., formaseries 
of semiconducting compounds, some of which,under 
certain circumstances, are amorphous. The physi- 
cal properties of these (conductivity, activation 
energy, thermal emf, etc.) are characteristic of 
crystalline substances [1]. 

The aim of the present work was to determine 
the short-range order, which could have a substan- 
tial bearing on some of the physical properties of 
GaSe in the amorphous state. 

By studying the structure of GaSe films vacu- 
um-deposited onto a celluloid substrate, it was 
found that the GaSe condensed in the form of anex- 
tremely stable amorphous phase. The films re- 
mained amorphous after heating in vacuum for 
many hours. Crystallization took place more easily 
in the thicker films. In order to obtain crystalline 
samples it was essential to preheat the celluloid 
film to 120-130°C. 

For electron-diffraction examination we used 
GaSe prepared by fusing the elements (taken in 
stoichiometric ratio) in an evacuated quartz am- 
poule and subsequently cooling to room temperature 
[al 

Large pieces of this alloy were placed in acone- 
shaped tungsten spiral and vacuum-evaporated onto 
hot or cold celluloid films. The diffraction picture 
obtained from this specimen consisted of four dif- 
fuse rings (see Table 1). 
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In order to obtain an experimental intensity 
curve the multiple-exposure method was used for 
the diffraction pictures [2]. 

Six electron-diffraction photographs were ob- 
tained with multiple exposures and microphoto- 
metric curves were obtained on the MF-4 
photometer. After constructing blackening curves 
and determining intensities for various values of 
(sin’)/A from these, we obtained an experimental 
intensity curve for amorphous GaSe in arbitrary 
units. 

The application of the integral-analysis method 
in electron diffraction to complex substances includ- 
ing atoms of various kinds demands special theo- 
retical treatment. The case of GaSe may, however, 
be simplified. Remembering the proximity of the 
respective atomic numbers (Ga, 31; Se, 34) and 
scattering powers, we may regard this compound as 
consisting of atoms of one sort. 

The experimental intensity curve must be nor- 
malized by a curve representing the sum of coher- 
ent and incoherent scattering 


TABLE 1 
Ring 


No. Intensity 


(sin @)/A 


Weak 

Very very strong 
Weak 

Very strong 


0.4110 
0.166 
0.276 
0,430 


FN ces 


424 Ig bc 


S18) Gq + (28 15")se 


fea ale se 


(1) 


) 


where fGa,Se is the atomic scattering factor of elec- 
trons for Ga and Se; Z is the atomic number; ZS is 
is the incoherent-scattering function; s = (4sin¥ YA; 
23 is the scattering angle; » is the electron wave- 
length. 

In expression (1), the first term represents co- 
herent scattering and the second (calculated from 
data of Bewilogua [8]) incoherent scattering. The 
curves were joined at (sin ®)/A =~ 0.8. After nor- 
malization and subtraction of the incoherent scat- 
tering from the intensity curve (Fig.1), the quan- 
tities under the integral sign in the expression 


5 Ps 2ruene é 
Acro (r) = 4nr?o9+ a \ si (s) sin srds 
0 
could be determined. 
Here, r is the distance in A, p(x) is the atomic- 
density function, and py is the mean atomic density 
of the given substance 


re 5.03 g/em* 
M0 148.68-1.65- 10-24 


= 0.0205 at, /A’, 


For calculating the integral we used the simpli- 
fied method proposed by Harris, Wood, and Ritter 
[4]. The integration reduces to the summation of a 
large number of terms. The whole of the region s 
used is divided into small intervals, e.g., As = 
0.02. The integrand comprises the product of two 
factors si(s) and sinsr; the si(s) is a curve ina 
range varying, for example, from 0 to 10.05,which 


Texp 
| 6 
2 
Bae 
mm 


2 


Y 


wq 
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determines the practical limits of integration. Ifwe 
choose r so that for any small value of s the prod- 
uct rs = 27, then sinsr, sin2sr, sin3sr,... will 
be zero and at intermediate points the sines will 
have recurring positive and negative values. For 
example, if r= 5.0671, then for s = 0, 1.24, 2.48, 
...,9.92, the product rs is a multiple of 27. The 
whole range of s considered will have 8 cycles. 
However, we have only to calculate 15 sine values 
in a af period. The multiplication by si(s) reduces 
to addition of coefficients for the same values of the 
sines. 

It is still simpler to prepare strips represent- 
ing the product of sines by si(s) where si(s) varies, 
for example, from 1 to 100 in certain small inter- 
vals. 

Having calculated the integral in this way, we 
obtained a radial distribution curve (Fig.2). The 
first two maxima lie at rj = 2.40 and ry = 38.75 A. 
The area under these corresponds to approximately 
6 and 12 atoms. In order to interpret these results 
we turn naturally to the lattice of crystalline GaSe, 
the structure of which we studied earlier together 
with that of the amorphous type [5]. The GaSe lat- 
tice has distances Se—Se = 2.33 and Ga—Se = 2.48 
A, averaging 2.405 A. We may suppose that both 
these distances are reflected in the first maximum. 
In the GaSe lattice, however, each Ga or Se atom 
has three Se or Ga atoms around it at distances of 
2.48 A. The interatomic distance 3.75 A is the lat- 
tice spacing a. In the plane of the film there are 
six atoms of the same type at this distance from 
any one atom. Furthermore, there are three Ga 
atoms around each Ga at a distance of 3.84 A. We 
may suppose that the second maximum reflects both 
these distances. Thus, according to our data, there 
is no correspondence between the first coordination 
groups of the amorphous phase and the crystalline- 
phase lattice. Perhaps we should not have expected 
such correspondence, since it may well be that a 


Fig. 2 
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change of coordination takes place on passing from 
the amorphous to the crystalline state without any 
change in the composition. 
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GaSe films obtained by vacuum evaporation onto celluloid substrates and subsequently 
annealed were subjected to structural study. The hexagonal unit cell had a = 3.74, 
c =15.89 kX. The space group was Déh —C6mme;the atoms lay in position 4(f ), 


Wea = 103 “ee = OI. 


The structure was layerlike. 


The layers were composed 


of pyramids turned with the (Se) vertices toward one another, the bases being Ga 


atoms. 


The absence of many reflections from the electron-diffraction patterns and 


the diffuseness of others indicated defects in the lattice structure caused by the dis- 
placement of layers perpendicularly to the c axis by a vector distance 4x +Ay. 


The crystal structure of phases belonging to 
the Ga—Se system has been studied several times 
by x rays [1]. Apparently there are one or two 
phases corresponding to the ideal composition GaSe 
in the concentration range 48 to 52 at.% Se. Accord- 
ing to x-raydata [1], a rhombohedral structure pre- 
dominates for Se > 50 at.% and a hexagonal struc- 
ture for Se < 50 at.%. Samples were prepared by 
fusing the elements in vacuum; the processes of 
ordering and crystallizing take place very slowly 
(even at high temperatures) and some of the maxima 
on the x-ray diffraction pictures are diffuse. The 
two structures are very similar, and in powder 
photographs the occurrence of the rhombohedral 
phase is indicated by splitting of the (0113) line of 
the hexagonal lattice. 

The hexagonal-lattice spacings are a = 3.73; 
and c = 15.88, kX. There are four GaSe molecules 
in the unit cell. In describing the structure, the 
authors of [1] presented the concept of lattice de- 
fects consisting of sections with rhombohedral 
stacking forming part of the hexagonal phase. The 
space group Ciy = CG was selected; this enabled 
the coordinates of two molecules to be described by 
setting them in positions 2(g) and 2(i), the other two 
molecules requiring translations 7, 44, 4. From 
the data presented we find that a considerable num- 
ber of reflections which would have been expected 
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to occur are in fact absent from the x-ray diffrac- 
tion photographs. 

We notice the similarity between the unit-cell 
dimensions given for the hexagonal phase and those 
of GaS [2]: @ = 3.57, and c =15.4, kX. The GaS 
structure is described by the group Dén: the atoms 
of Ga and S occupy position 4(f) with parameters 
ZGq = 0.17 and zg = 0.60. From the data on GaSe 
we see [2] that the x-ray diffraction photographs of 
the hexagonal phase show no reflections contradict- 
ing the extinction conditions for Deh: Hence, we 
may expect the structures of GaS and GaSe to be 
very similar. 

In the present investigation we used samples of 
GaSe! prepared by fusion in a quartz ampoule under 
vacuum, taking the elements in stoichiometric ratio 
[3]. The specimens for electron-diffraction study 
were obtained by evaporating the GaSe onto a cellu- 
loid substrate. The first experiments led to the 
formation of an amorphous deposit. In order to ob- 
tain a crystalline film, the celluloid substrate was 
preheated to about 120°C. Crystallization took 
place more rapidly (3-4 h, subsequent to heating for 
thicker films. Thin films remained amorphous even 
after 10 h. Specimens for study in transmission 


“Institute of Physics, Polish Academy of Sciences, Warsaw. 
1 The compound was obtained from B.T, Kolomiets and N. A. 
Goryunova, Leningrad Physicotechnical Institute. 
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Fig, 1. Electron-diffraction photograph of GaSe. 


were also prepared from a suspension in water.The 
orientation in these was worse; hence, for obtaining 
oblique-texture diffraction patterns, films prepared 
by sublimation were used. From these diffraction 
patterns (Fig.1) we see that the processes of order- 
ing and crystallization failed to pass to completion, 
since two types of reflection occur: sharp and dif- 
fuse. The complete set of reflections corresponds 
to spacings of a = 3.74 and c = 15. 889 kX of the 
hexagonal unit cell, in agreement with the x-ray 
data. Diffuse reflections only occur on ellipses with 
h —k # 3n, the reflections of hki0 type on these hav- 
ing normal sharpness. 

In the first stages we made an attempt to cal- 
culate the structure from the sharp reflections only, 
these corresponding to spacings a = 3.74 andc = 
7.945 kX. The results obtained from the complete 
set of reflections show that, by analyzing the data 
with the sharp reflections only, we obtain the struc- 
ture of a single "packet" of the type making up the 
GaSe lattice (see below). 

The electron-diffraction pictures show the fol- 
lowing extinctions: hkil exist in all orders, 000/ 
only for 1 = 2n, hh2h/ only for 1 = 2n; of the hkil 
with h — k= 8n, only those for which / = 2n exist. 
The number of extinctions exceeds that which might 
be expected for any of the hexagonal space groups. 
Some of these extinctions may be connected with 
special positions of the atoms and perhaps other 
features of lattice structure; they are not contra- 
dictory to a number of groups, including space 
group Dene By analogy with GaS, we may site the 
Ga and Se atoms in the 4(f) position: 


*/s, */s 2, "Is “a2, */s, iis (2 ee fe ops ‘sz. 


We were able to locate 54 independent reflec- 
tions on the photographs (including 7 0001). 

The intensities of the sharp reflections were 
estimated visually and those of the diffuse reflec- 
tions on an objective nonrecording MF-2 micro- 
photometer. The intensities of the diffuse and 
sharp reflections were compared by measuring the 
half-widths of the maxima and equating these to the 
half-widths of the sharp reflections by introducing 
an appropriate factor. 


Fig. 2. Cross section of ® series (potential) of 
the GaSe structure along the x0z plane. 
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TABLE 1 
PI rere nn ae Se 
BUG: hkil @D exp @D theor ot hkil @D exp Dd theor 
4 1010 2.32 4.83 24 303.12 (), 86 0.83 
2 1013 4.86 3.95 25 303.14 1.16 1.40 
3 1015 (Lag 1.75 26 2240 2.68 2.02 
4 1017 4.04 A 33 27 2244 1.83 0.99 
5 1420 5.07 5.16 28 2248 1.20 0.53 
6 1122 ().57 (), 84 29 224.10 (),87 0,92 
a 1124 3,56 ven 30 224.12 0.97 0.73 
8 1128 0.54 Lu 3t 22414 0.99 4,23 
i) 112.40 1.66 1.7) 32 3140 1.36 0.95 
410 ian 1.76 1.24 33 3143 0.84 0.91 
14 412.14 1.86 1.95 34 3145 0.46 0.50 
{2 2020 1.95 ot 35 3147 CBT 4.45 
43 2023 48, 2.04 36 4040 143 0.80 
{4 2025 0.89 1.08 37 4047 435 4.24 
15 2027 2.30 2.86 38 3250 1.00 0.70 
6 2130 1.76 1.55 39 3257 4.64 1.09 
17 2133 1.66 1.45 AV A150 4.23 eo 
18 2135 0.54 (0,78 Al AV5A 1.07 0.65 
19 2137 t.54 2.14 42 415.10 0.85 0.65 
20 243.11 0.57 0.76 43 415.12 4.43 1.96 
24 3030 2.93 Dh its) 44 445.14 1.16 0.94 
> 3034 2236 123) 45 5050 1.09 0.5% 
23 303.10 (0,83 1 08 46 3360 1.14 1,07 
AT 4260 4.5953 (52 


From experimental values of 6? we calculated 
the cross section of the 6? series in the (00z) di- 
rection. The resultant interatomic vectors enabledus 
to site the atoms in the unit cell and make a prelimin- 
ary estimate of the parameters. These were refined 
by calculating the series for the potential cross 
section in the x0z plane (Fig.2). In order to deter- 
mine the signs of the structure amplitudes we used 
parameter values Zge = 0.17, zG@qg =—0.10. The 
parameters found from the synthesis were zSe = 
0.177 and zGa = —0.10. The R factor obtained from 
47 reflections (without 0001) was 25.7%. We see 
from the table that the reflection 1010 gives the 
greatest deviation from theory, which may be ex- 
plained by the difficulty of estimating intensities so 
close to the central beam. If this reflection is ex- 
cluded, we have R= 23.0%. 

The structure of GaSe may be described in the 
following way. As indicated by Figs. 3a and 4, the 
alternation of the plane hexagonal lattices of Ga and 
Se atoms along the c axis justifies us in referring 
to a layer-type structure formed by "packets" with 
an interval of 3.18 A between them (Fig. 3b). Inside 
each packet, pyramids with their bases at the boun- 
daries of the packets are turned with their vertices 
inward (Fig.4). Each Se atom has three Ga atoms 
at 2.48 A and one Se at 2.325 A as its nearest neigh- 
bors. The very similar scattering powers of the 


two kinds of atoms prevent us from choosing be- 
tween the form shown in Fig. 4 and another in which 
Ga and Se are interchanged from experimental data 
only. The arrangement shown in Figs. 3a,b is pre- 
ferable, since the distance between atoms a and a' 
is close to the minimum Se—Se distance in the Se 
lattice (2.32 A) and differs from the corresponding 
Ga-—Ga distance in the Ga lattice (2.45 A). We con- 
sider that the difference between the chemical func- 
tions of Ga and Se excludes the possibility of these 
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Fig. 3, Scheme of GaSe structure: (a) spatial 
arrangement of atoms; (b) section of the struc- 
ture in a (1120) plane. 
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Fig. 4. Coordination of atoms in the GaSe 
structure, 


atoms having a disordered arrangement in the lat- 
tice. 

We failed to observe the rhombohedral struc- 
ture found in the x-ray work, although all the hex- 
agonal-phase reflections for which h —k = 8n satis- 
fied the rhombohedral condition h—k+l = 3n. Re- 
flections with h—k # 3n, however, did not satisfy 
this condition. 

The structure described may be regarded as 
ideal. The model was obtained without allowing for 
the absence of a large number of reflections for 
which the calculated intensity was nonzero, and al- 
so without allowing for the diffuse nature of a cer- 
tain group of reflections. 

As mentioned earlier, all reflections with l = 
2n are absent from ellipses with h — k # 3n, and all 
present are diffuse. A possible explanation for 
these peculiarities is based on a theory relating to 
disorder in the structure of hexagonal cobalt [4]. In 
hexagonal cobalt the structural distortion is con- 
nected with the displacement of some layers by a 
vector length 44x +"4y. Asa result of this displace- 
ment, reflections with h — k = 3n and / ~ 2n drop 
out, while for J = 2n they have normal intensity and 
sharpness. Reflections with h—k ~ 3n become dif- 
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fuse and elongated along the c* direction; the dif- 
fuseness of reflections with / = 2n is three times 
greater than that of reflections with / ~ 2n, which 
explains the absence of the former from the elec- 
tron-diffraction patterns. All hki0 reflections re- 
main sharp. These features of the electron-diffrac- 
tion picture correctly reflect those found for the 
case of GaSe. 

The shift of 44x +*Ay in the structure of GaSe is 
entirely possible in view of its layerlike arrange- 
ment. Thus, the structure of our GaSe specimens 
cannot be characterized by any particular space 
group in view of the absence of strict periodicity 
along the c axis, i.e., periodicity in the positioning 
of the layers. At the same time, the projection of 
the assembly of superposed layers may be charac- 
terized by a plane translational group C3 — ps. 

It is interesting to compare the disorder of the 
GaSe structure here described with that found in 
layered lattices of Cdl,, CdBr,, etc. [5]. In contrast 
to the variable structure in Cdl, and others, GaSe 
does not apparently contain inclusions with a defin- 
ite spatial periodicity along the three directions. 
Packets of structure displaced from their "regular" 
position by 4x +Ay are distributed at random over 
the thickness of the crystal (along the x axis). 
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The structure of sublimed germanium films was studied as a function of the character 
of the substrates, their temperature, the rate of deposition, and the thickness of the 
films. Various types of crystallite orientation were found; these were associated both 
with the presence of the plane substrate and the rectilinear propagation of the atomic 
beams. Highly oriented films (single-crystal type) could be obtained by evaporating on- 


to germanium single crystals. 


It is well known [1] that when germanium is 
evaporated in vacuum onto a substrate held at a 
temperature below 370°C an amorphous film is 
formed; in order to obtain crystalline films, the 
substrate temperature must be raised above 370°C. 
Very different results were obtained by the author 
of a recently published paper [2], in which it was 
asserted that, on vacuum-depositing germanium on- 
to cleavage faces of zinc blende, crystalline films 
were formed even at 250°C. The reduction of the 
critical temperature to 250°C was explained in [2] 
as being due to the similarity between the lattice 
spacings of Ge and ZnS (5.63 and 5.42 A) and be- 
tween their respective chemical bonds. 

In the present investigation we studied the struc- 
ture of thin germanium films as a function of the 
nature and temperature of the substrate, the rate of 
deposition, and the thickness of the films. The ger- 
manium was vacuum-evaporated from tungsten 
evaporators heated by electric currents. Graphite 
evaporators heated by electron bombardment were 
also used. 

For the condensing substrates we used glass 
and quartz plates, corundum, carborundum, mica 
sheet, and cleavage faces of rocksalt, calcite,fluor- 
ite, and silicon. We also made a large number of 
experiments involving deposition on single-crystal 
germanium surfaces after grinding and etching. 

The substrate temperature during condensation 
was measured by an iron—constantan thermocouple 
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and in individual experiments varied from room 
temperature to 950°C. The film thickness was de- 
termined by weighing the substrates before and 
after deposition; it varied from 0.01 to 20 uw. The 
evaporation rate was determined from the film 
thickness and the time required to deposit the films. 
In individual experiments the evaporation rate 
varied from 107! to 2 - 107° g/cm? « sec, i.e., by 
150 to 200 times. 

In accordance with published data [1], on evap- 
orating the germanium onto cold (room tempera- 
ture) substrates, we obtained amorphous films. 
Crystalline films were obtained on evaporating ger- 
manium onto substrates heated above 370°C. 

Let us now consider the results obtained by 
studying the structure of the films. 

For rapid sublimation of germanium onto glass, 
polished quartz, or corundum (heated to 400°C), 
polycrystalline films were formed (Fig.1). Subse- 
quent annealing of these samples at 400 to 900°C for 
1 to 6 h had no effect on the structure of the films. 
In all the electron-diffraction photographs obtained 
with polycrystalline films, we observed the 222 
ring, which is forbidden by the structure factor; the 
presence of this was noted earlier [3] and associated 
with dynamic effects. 

It is well known [4] that thin films often develop 
texture, this being associated with the presence of 
the plane substrate and nonuniform development of 
the crystal faces. In the case of metal films ob- 
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Fig. 1. Electron-diffraction pic- 
ture of a polycrystalline Ge film. 


tained by sublimation, the formation of texture re- 
quires, in some cases, that the sublimation should 
be effected onto a heated substrate and in others 
that the samples should be annealed at some defin- 
ite temperature after sublimation [5]. 

In our experiments with germanium, textures 
of various kinds developed during rapid evaporation 
onto substrates heated above 500°C. Thus, for 
example, when germanium was deposited on po- 
lished graphite plates, calcite crystals, and silicon 
carbide heated to 550 or 600°C, textures compris- 
ing germanium crystallites oriented with the (110) 
and (100) faces parallel to the substrate (the first 
of these predominating) were obtained; when the 
germanium was evaporated onto glass at 520°C, 
texture with orientations (111) and (100) developed, 
while for evaporation onto polished corundum at 
700 to 800°C the (100) faces lay parallel to the sub- 
strate. 

Later, on studying films 10 to 20 uw thick, we 
found that the germanium texture was connected not 
only with the presence of the plane substrate, but 
also with the fact that, in the course of vacuum 
evaporation, there was a unidirectional flow of the 
material supplying the crystallites, owing to the 
rectilinear propagation of the atomic beams (so- 
called feeding texture). In contrast to the types of 
orientation mentioned above, in this case the crys- 
tallites set with their (110) face parallel to the sub- 
strate, or perpendicular to the direction of the 
atomic beam. 

The variety of textures formed on substrates 
by evaporated germanium is apparently connected 
with the isometric habit of the crystallites, by vir- 
tue of which the disposition of one or another of the 
faces parallel to the substrate may depend on anex- 
tremely small difference in the conditions of deposi- 
tion. 

In the case of feeding textures, only one un- 
varying type of orientation was observed; judging 
from the lengths of the arcs on the diffraction pic- 
tures (Fig.2), the orientation here was more per- 
fect than in the case of textures associated with the 
plane substrate. 


OF GERMANIUM FILMS 


431 


Fig. 2. Electron-diffraction pic- 
ture of Ge texture, 


Fig. 3, Point electron-diffrac- 
tion picture of a Ge film ob- 
tained by deposition on calcite. 


Apart from photographs of polycrystalline 
samples and textures, we also obtained point elec- 
tron-diffraction pictures indicating the presence 
and azimuthal orientation of germanium crystals 
relative to the substrates. The most perfect orient- 
ation occurred when depositing germanium on cleav- 
age faces of calcite and mica (Fig. 3). 

By solving the point electron-diffraction pic- 
tures we found that, in the case of calcite, the ger- 
manium crystallites were set with the (111) face 
parallel to the (1011) face of the calcite, for a sub- 
strate temperature of about 520 to 600°C, andsome- 
times with the (100) face so oriented. In the case 
of mica, the germanium crystals were also setwith 
the (111) face on the substrate, or more rarely the 
(110) face. Diffraction photographs of this type al- 
ways Showed weak Debye rings as well as point re- 
flections, indicating that some of the crystallites 
were not oriented. Point diffraction pictures were 
obtained for films up to 10~° cm thick. Thicker films 
only constituted textures and polycrystalline ag- 
gregates. We also evaporated germanium onto ger- 
manium single crystals, ground and subsequently 
etched. These experiments were aimed at eluci- 
dating the possibility of growing a germanium _ingle 
crystal by crystallization from the vapor phase in 
vacuum. Apart from its theoretical interest, this 
question has considerable practical importance. It 
is well known that a layer of germanium produced 
by vacuum evaporation has hole-type conductivity 
[6]; hence, by evaporating onto an n-type single 
crystal we might well expect to produce crystals 
with p-n junctions. In this respect, it is particular- 
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Fig. 4. Electron-diffraction pic- 
ture of an etched germanium 
single crystal. 


Fig. 5. Electron-diffraction pic- 
tures of germanium films ob- 
tained by deposition on the etched 
surface of a single crystal. 


ly interesting to study the structure of the layer 
growing on the single crystals and to determine the 
conditions under which it will acquire a single-crys- 
tal nature. 

The crystals were etched in hydrogen peroxide. 
Figures 4 and 5 show electron-diffraction photo- 
graphs obtained from one of the single crystals be- 
fore and after depositing a 2-u layer. The layer 
was deposited on the surface of a crystal heated to 
520°C, the 2-y layer being deposited in 30 min 
(slow evaporation). The fact that Fig.5 contains 
reflections forming a regular lattice indicates that 
under these conditions a layer with the structure of 
a mosaic single crystal forms on the original single 
crystal. The cross shape of some reflections indi- 
cates that the growing layer includes crystallites 
with octahedral faces, possessing insufficient re- 
solving power [5]. 

Figures 6 and 7 present electron-diffraction 
photographs obtained from another single crystal 
before and after depositing a 2-u layer. The crys- 
tals were heated to 520°C and the deposition took 3 
min (rapid evaporation). The presence of Kikuchi 
lines on the diffraction pattern of Fig. 7 indicates 
the formation of a single-crystal layer with a more 
perfect structure than that of the earlier case. 
Hence, increasing the evaporation rate improves 
the structure of the deposited layer. Increasing the 
substrate temperature leads to the same result. 

As indicated earlier, the author of [2] observed 
crystalline germanium films on evaporating ger- 
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Fig. 6. Electron-diffraction pic- 
ture of an etched germanium 
single crystal. 


Fig. 7. Electron-diffraction pic- 
tures of germanium films ob- 
tained by deposition on the etched 
surface of a single crystal. 


manium onto crystal faces of zinc blende heated to 
250 rather than 370°C, and explained this as being 
due to the similarity between the lattice spacings 
and chemical bonds in the two structures. This is 
rather unexpected. It is certainly well known that 
similarity between the spacings and bonds leads to 
the orientation of the evaporated crystallites on the 
crystalline substrate. As regards crystallization 
from the vapor phase during vacuum sublimation, 
so far as we know the formation of a crystalline or 
amorphous phase is determined by the substrate 
temperature and the time of holding the sample at 
some temperature rather than by the similarity of 
the lattice spacings or identity of chemical bonds. 
Elementary consideration of the crystallization pro- 
cess indicates that the formation of a crystalline or 
amorphous phase depends mainly on the mobility of 
the atoms; if the atoms on the substrate have low 
mobility and cannot assemble into aggregates (high 
viscosity), an amorphous phase is formed; if, how- 
ever, the mobility of the atoms is sufficient for the 
formation of nuclei, then a crystalline phase is 
formed. From this point of view it is quite clear 
that raising the substrate temperature increases 
the kinetic energy and mobility of the atoms, and 
thus favors the formation of a crystalline phase. 
As regards the closeness of the lattice spacings 
and identity of the bonds, it is not really quite clear 
bow these can facilitate the formation of a crystal- 
line phase. If identity of bonds and similarity of lat- 
tice spacings facilitated the formation of a crystal- 
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line phase, then, clearly, by substituting germani- 
um crystals for ZnS (depositing germanium on ger- 
manium) we could still further lower the tempera- 
ture for the formation of the crystalline films. 

In view of this, we made some special experi- 
ments and found that, on depositing germanium on 
germanium single crystals heated to 250 or even 
370°C, only amorphous films were formed. Crys- 
talline films formed above 370°C. Thus we must 
regard the results of [2] as erroneous. 


LITERATURE CITED 


1. G.Haas, "Studies on the structure and be- 
havior of thin evaporated films of silver, 


aluminum, silicon, and germanium," Phys. 
Rev. , 72, 174 (1947). 

A. Segmiiller, "Texturuntersuchungen an 
diinnen Zinn-und Germaniumschichten," Z. 
Krist. ,107, 1/2, 18-35 (1956). 

R. Heidenreich, "Theory of 'forbidden' (222) 
electron reflection in the diamond structure," 
Phys. Rev. ,77, 2, 271-284 (1950). 

Baise Vainshtein, "Electron-diffraction tex- 
ture photographs," Tr. Inst. Kristallogr. ,6, 
173-193 (1951). 

Z.G. Pinsker, Electron Diffraction [in 
Russian] (1949). 

J. Thornhill and K. Lark-Horovitz, "Films of 
germanium," Phys. Rev. , 82, 762-763 (1952). 


SOVIET PHYSICS —- CRYSTATLOGRAPRHY VOL 1,7 NOs t5 SEPT .-~OCT., L206 


MEASUREMENT AND RECALCULATION OF THE 
DYNAMIC ELASTIC CONSTANTS FOR QUARTZ 
V. G. Zubov and M. M. Firsova 


M.V. Lomonosov State University, Moscow 
Translated from Kristallografiya, Vol. 1, No. 5, 
pp. 546-554, September-October, 1956 

Original article submitted July 13, 1956 


The dynamic constants have been measured by Bergmann's method, with independent 
and direct determinations of c1;, C33, C44, and Cgg. It is shown that Lawson's piezo- 

electric corrections are needed. Discrepancies are observed between the C44 and cy, 
determined in different ways, but the Laval-Raman theory eliminates the discrepan- 
cies and gives a fuller description of the dynamic behavior. 


Introduction plates; the other constants were found by calcula- 
tion. In particular, their values for cg3 and cy3 are 
very unreliable, as very complicated equations are 
involved. Mason deduced cCegg, C4y, and cy, from the 
thickness modes ofA T-cut plates, the others being 
calculated from the sjk. Vogt's theory gives equa- 
tions different from those of the Laval-Raman 
theory, so these values cannot be used. 

Nomoto used Bergmann's method to deduce all 
the constants from the diffraction figures given by 
ultrasound in a quartz cube. In this way he found 
C33 independently, as well as cy, and Cg. 


Much attention has recently been given to the 
limitations of Vogt's theory of the dynamic elastic 
behavior of crystals and to the construction of a 
new atomistic theory. Laval and Raman have made 
major advances in this respect. 

Laval [1], in 1951, and Raman [2], in 1955, 
employed very general concepts on particle inter- 
actioninrather different demonstrations that Vogt's 
theory is not applicable to the general case of dy- 
namic processes in crystals, which for complete 
description require 45 independent constants in- 
stead of Vogt's 21. The Laval-Raman theory also The errors in the directly measured constants 
predicted discrepancies between static and dynamic do not exceed 0.5%, whereas those in the calculated 
measurements of the adiabatic constants; it implied eS are 0.5-1%, except for Nomoto's cyy and 


a markedly different relation between the elastic Atanasoff and Hart's c33 and cy3, which give over 1%. 
moduli s;, and the elastic constant ¢jk. Table 1 lists the most reliable effective values 
The new theory naturally demands a fresh for the Cik. 
examination of the existing experimental evidence The values of Table 1 take no account of the 
and much more careful new direct measurements differences in the experimental conditions; no piezo- 
of elastic constants, especially for quartz, a ma- electric corrections have been applied. Atanasoff 
terial of great technical importance. and Hart obtained two different values for cy, by 
It has been shown [3] that there is no agree- measurement of natural frequencies with different 
ment on the numerical values for the elastic con- circuits; Lawson's [7] corrections eliminated the 
stants of quartz; direct measurements independent discrepancy. These corrections are required not 
of results for other constants are available onlyfor  OMly for resonance in plates, but also in Bergmann's 
C11 and Cgg, even now. method. The most reliable values, as thus cor- 
The most careful and precise dynamic deter- rected, are then as given in Table 2. 
minations of the adiabatic constants are those of Table 2 also lists adiabatic values deduced 
Atanasoff and Hart [4], Mason [5], and Nomoto [6]. from Vogt's static measurements; these deviate 
Atanasoff and Hart found only cy, and cee directly by widely from the dynamic results in nearly allcases, 
measurement of vibrational frequencies for quartz and the discrepancies are not explicable on Vogt's 
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TABLE 1 
so ah a a ret A ee 


en, fore dyne/cm? 
Source Se ee 


Atanasoff and Hart Sipe ip 
Mason — = 
Nomoto * 


o7. 19 


Ai ST 2 oes. 40 
58.65 d 


40.5 {8.2 
— Lion 


5) 


oy eee 


*From the original paper [6]; other values are usually given in the 
literature, as a result of incorrect normalization of Nomoto's re- 
sults, Bhagavantam's values could not be included, as we have 
been unable to obtain the original paper. 


TABLE 2 


Cik* io; dyne/cm? 
Source = 


Atanasoff and Hart 86.75 = 08.0 39.9 17.97 
Mason — — 58.6 401.3 18.2 
Nomoto * 86.5 LOS TeG — hoe) 
Vogt 85.4 | 405.¢| 57.1 | 30.1 | 16.9 


*Nomoto's cy, has not been corrected, because the error of meas- 
urement exceeds the correction. 


theory, whereas they agree with the Laval-Raman 
theory. There is also fairly poor agreement over 
the dynamic cy, and cy4,, which cannot be explained 
in Vogt's theory. 

Here we describe new independent measure- 
ments of some of the constants, which were made 
in order to clear up these discrepancies. 


Methods 


We used Bergmann's method [8], because it is 
far more flexible than the purely electronic meth- 
ods. A particular advantage is that all the con- 
stants can be determined on the same specimen; 
further, it is possible to use in the calculations 
modes that do not allow of direct piezoelectric ex- 
citation, and that cannot be detected by electronic 
methods. This gives independent access to C33, C44, 
and Cy, aS well as ci; and c¢gg; in addition, the equa- 
tions are much simplified, and it is possible to 
check that the orientation is correct by reference 
to invariant relations. 

Figure 1 shows the system used to observe the 
diffraction by the ultrasound, which produces very 
complex diffraction patterns in the focal plane of 
O,, which reflect the symmetry of the crystal and 
the presence of three types of waves along each di- 
rection. 

Figure 2 shows some of the patterns recorded. 
Each pair of spots here corresponds to a plane wave 
whose wavelength A may be calculated from 
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Fig. 1, The apparatus: L) light source; K) condenser lens; F) 
filter; D) pinhole; O,) collimator lens; A) specimen; O,) lens; 
B) screen or plate in focal plane of lens, 


Fig. 2. Diffraction patterns for various di- 
rections of observation: a) along Z axis; 
b) along Y axis; c) along X axis; d) at 45° 
to the X and Z axes. 


1 K 
Pig == ll 0 Ne 
in which ry is the distance from the central maxi- 
mum to that of order n, F is the focal length of the 
lens, and A is the wavelength of the light. 

The A are used with independent measurements 
of the frequency and density to find the effective 
elastic constants corresponding to these waves, and 
hence the principal constants. 

Atanasoff and Hart derived [4] the following 
equations for the thickness modes from Vogt's 
theory: 

a) waves along the X axis: 


Cy, —-k =0, (X,) 
(Gay h) (Ga k) eC == () (X», X3) 
b) waves along the Y axis: 
(G4) (Ca *) Grae, (Y1,Y2) 
Cog — k = 0, (Y3) 
c) waves along the Z axis: 
Cy; —k = 0, (Z1) 
Cy _ k ==) 0, (Zy) 


(X,), (Y,), and (Z,) give the effective coefficients 
for quasi-longitudinal waves; the others are for 
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quasi-transverse waves. We have used these equa- 
tions with our data. 

The light source was a mercury arc type PRK 
or SVDSh, which was filtered to isolate } = 5461 A; 
Lenzos interference filters were used, the band 
width and peak transmission being checked from 
spectrograms. O, had a focal length of 167.0 cm; 
O, had 50 cm. The pinhole was 0.01 mm in diam- 
eter. The patterns were measured to 2 - 10-* mm 
with an IZA-2 comparator. 

These readings were checked from micro- 
photometer traces; corrections were made in cer- 
tain cases by measuring the distances between 
higher-order spots. 

The waves were excited either mechanically 
(with a quartz piezoelectric plate) or electrically 
(directly in the specimen); both gave the same 
values for the effective constants. The power was 
provided by the usual three-point oscillator: output 
power 80 W, peak output 3000 V, frequency range 
8-12 Mc.! These frequencies correspond to the 
30-40th harmonics in the specimen, so edge effects 
are small. The frequencies were measured to 
0.02% with a heterodyne wavemeter. 

The specimens were made at the Institute of 
Crystallography from fine natural quartz as cubes 
20 X 20 X 20 mm oriented on the principal axes; 
they were carefully checked for orientation and ab- 
sence of twins.” 

The diffraction patterns also provided a test 
for twinning, while the invariants [3] provided a 
check on the orientation. 

The constants should be obtained with errors 
not exceeding 0.2% under these conditions. 

All results quoted for velocities and elastic 
constants at room temperature are derived from 
400-600 independent measurements; these numer- 
ous data provided a good test of the performance, 
and the individual values seldom deviated by more 
than 0.5% from the mean, except when photographs 
of poor quality were used. 

The final values are means from six cubes and 
nine plates; here again the deviations from the mean 
did not exceed 0.5% for properly oriented speci- 
mens. 


Results and Discussion 


The measurements were made at 20°C for 


waves propagating along all the principal directions. 


Table 3 gives the effective elastic constants k for 
the various modes. 

Calculation from Vogt's theory gives C41, Ceg; 
C33, and Cy, directly from the k for (X;), (Y3), (Z;), 
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and (Z,), respectively; the coefficients for (X_) with 
(X3) or (Y;) with (Y,) may be taken in pairs to give 
C14 and cy, by Atanasoff and Hart's method, and al- 
so to test the results for consistency. Table 4 pre- 
sents this test. 

Tables 1 and 4 show that our cj; and Cgg agree 
satisfactorily with the direct and independent de- 
terminations by others. 

Tables 1 and 4 also show that the cy, calculated 
in the various ways are in conflict, as in Atanasoff 
and Hart's work. This conflict, together with the 
agreement over C1, and Cg, for identical electrical 
conditions, would indicate that results from Berg- 
mann's method require piezoelectric corrections 
such as Lawson applied to Atanasoff and Hart's re- 
sults. Table 5 gives our results for 20°C and fixed 
electric field as corrected in this way. 

These corrections give good agreement be- 
tween the cy, as found in the various ways, but left 
unaltered the spread in cy,. 

Comparison of Tables 2 and 5 leads to the fol- 
lowing conclusions. 

1. Our results show that static measurements 
give a large and persistent deviation from dynamic 
ones for all the adiabatic constants except c33. Our 
data on these disagree with Vogt's and with all 
others. 

2. Our Cg3 agree closely with Nomoto's (also 
determined directly) and do not confirm the c33 = 
106.8 suggested by Atanasoff and Hart, and also by 
Cady [9]. The complicated calculation evidently in- 
troduces some gross error; it is better to reject 
this value and to use C33 = 105.5 as the most prob- 
able. The static and dynamic values here are in 
agreement. 

3. Our cy, (found directly and independently) 
agrees well with Mason's value, whereas calcula- 
tion via Atanasoff and Hart's equations gives good 
agreement with their value, the difference between 
the calculated and directly measured c4, exceeding 
the errors of measurement. The difference be- 


TABLE 3 


Mode 


cx | (X2) | (Xs) (Ys) | (Z,) | (Z2) 


(Y)) | (Y2) 


Effective elastic 
constants k - 10729 
dyne/cm? 


87,38] 68.0} 28.92] 94.96 | 49.47] 40,37| 105.4] 58.4% 


1We are indebted to E. N. Volkova for assembling all the electron- 
ic equipment, 


2We are indebted to A. B, Gil'varga for excellent work on the 
preparation and testing of the specimens. 
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TABLE 4 
Cu Cas Coe Cys Cig 
Independent and 
direct determination 87.38 105.4 40.37 58.4 — 
From (X2} and {Xs} — _ -- 16,55 (7.8 
From (Y;) and (CY, - — -- OOD 16,95 
TABLE 5 
| Cu C33 Cae Cy4 Cys 
Independent and 
direct determination 86.6 105.4 39.6 08.4 -- 
From (X,) and (X4) = = — DT ap (74 
From (Y,) and (Yq) i = ee 97.89 17.6 
Mean — — =: 76 Wied 
TABLE 6 
Ny Nga Nes N55 Nee Nis Thy 
Our data BG | MO ISA) TS SOG TES 
Atanasoff and Hart 86.75 — == SOLO SOO 8x4 A797 
Mason — — 98.6 | — | 40%35) 48.2 = 
Nomoto 86.5 105.7 -— OO) — = Alife 
Mean 86.6 105, oer p Oly MS ESE Bet yy 


tween our cy, and Mason's also exceeds the error 
of measurement. 

We thus have here a new conflict that cannot be 
explained on Vogt's theory, as for the conflict be- 
tween the static and dynamic measurements. 

The Laval-Raman theory explains the discrep- 
ancies for cy, and cy; it replaces Vogt's matrix for 
the elastic constants by the form [10]: 


Nyy Ryo yy yy O 0) Nig WV ) 
Nyy Ny ) 
N3y 0 0 () Q) 9) () 


The effective constants for (X_), (Xs), (Y;,), and 
(Y,) are dependent on nss and nz [3], while those for 
(Z.) and for shear waves in AT sections are depend- 
ent only on ny, and ny. Then Atanasoff and Hart's 
method for cy, and cy, in fact gives ns, and njqz; di- 
rect measurement of cy, and cy, or use of Mason's 
method, gives ny, and ny. The conflict is thus to 
be expected on the Laval-Rarnan theory. 


The various previous results can then be inter- 
preted to give the constants at room temperature 
and constant electric field (Table 6). 

This table shows that the Laval-Raman theory 
gives much better agreement between results than 
does Vogt's theory; the predicted difference be- 
tween ny, and nss, or between ny, and ny7, is actual- 
ly observed. We have also shown that the differ- 
ences persist at all temperatures up to the transi- 
tion point. The better agreement and lack of con- 
flict show that the Laval-Raman theory gives a ful- 
ler and more exact description of the dynamic elas- 
tic behavior. 

Raman drew the same conclusion from the data 
on cubic crystals. 

The table also shows that the Laval constants 
Mads N55, Nyg, and ny7 must replace Cy, and cy, when 
the usual equations are used to find the natural fre- 
quencies of quartz plates in any section. 


Conclusions 


Bergmann's method provides direct and inde- 
pendent measurement of C33 and c4y4, as well as cy; 
and cgg, the errors being less than 0.5%. The re- 
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sults need Lawson's piezoelectric corrections, on 
account of differences in the conditions for the vari- 
ous elastic constants, which is clear from a com- 
parison of the results with those obtained in other 
ways. 

Our results confirm that the static adiabatic 
constants of quartz inevitably deviate from the dy- 
namic ones; they also reveal a conflict (not pre- 
dicted by Vogt's theory) between the cy, and cy, de- 
rived from different natural vibrational modes of 
quartz. 

The Laval-Raman theory is applied to remove, 
these conflicts and to give good agreement between 
the various sets of results. It is concluded from 
these successes that the Laval-Raman theory gives 
a fuller and more correct description of the dyna- 
mic elastic behavior of crystals. 

We are indebted to A.V. Shubnikov and I.S. 
Zheludev for direction and discussion, and to Pro- 
fessor Laval of the College de France for providing 
access to his and Le Corré's data. 
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Tests do not confirm Bergmann's result for the ratio of peak intensities in diffraction 
of light by longitudinal and shear waves; a new effect is observed in the diffraction at 
fast shear waves propagating along the X axis of quartz. 


A crystal can carry three waves along any di- 
rection: one quasi-longitudinal (compression wave) 
and two quasi-transverse (shear waves). The dis- 
placements occur in three mutually perpendicular 
directions, and the velocities are usually different. 

An assertion often encountered in the literature 
is that the diffraction at shear waves is much 
weaker than that at longitudinal ones; Bergmann [1] 
has repeatedly asserted that the diffraction spots 
corresponding to shear waves must always be 
weaker than those corresponding to compression 
waves, and also that it is almost always impossible 
to detect higher-order diffraction peaks from shear 
waves. 

But we observed (in a study of the elasticity of 
quartz via diffraction at ultrasonic waves [2,3]) 
that this is far from being the case; the intensity 
ratio in the peaks is very much dependent on the 
resonance and excitation conditions. Mechanical ex- 
citation (from another quartz crystal) usually pro- 
duced weaker spots from the shear waves, butelec- 
trical excitation enables one to obtain any intensity 
ratio. Slow shear waves along the X and Y axes, 
and shear waves in 45° XZ sections, are most 
strongly excited and give the best patterns, reso- 
nance, and excitation conditions being the same. 
These waves produce higher-order diffraction 
peaks as easily as do compression waves; ones up 
to the order of 16 are readily observed even in 
fairly weak fields. 

Figure 1 shows apatternwith high-order peaks 
produced by slow shear waves propagating along 
the X and Y axes. Figure 2 shows a similar pat- 
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tern from compression waves propagating along the 
X axis, which required stronger fields than did the 
pattern of Fig. 1! 

This erroneous assertion may have arisen be- 
cause coupling to shear waves is very difficult in 
mechanical excitation. This difficulty would lead to 
unequal amplitudes and hence to a systematic dif- 
ference in the ratio of the intensities. 

Slow shear waves along the X axis show an ef- 
fect that does not fit the Raman-Nata treatment of 
diffraction intensities. 

Figure 3 shows the pattern seen along the Z 
axis in quartz when that axis coincides exactly with 
the incident beam. Slow shear waves along the X 
axis then produce peaks of even order very much 
stronger than those of odd order; the even ones are 
so strong that secondary peaks resembling the basic 
pattern are formed around them. A deviation of 
only 1° from the correct setting of the Z axis causes 
a pronounced change in the intensity ratios; the first- 
order peak becomes strongest, the intensity falling 
monotonically toward the higher orders. This fall 
is the same on both sides of the zero order, so that 
there is no unilateral extinction of peaks. 


It is difficult to explain the diffraction seen 
along the Z axis for fast shear waves propagating 
along the X axis. 

Figure 3 shows that the peaks corresponding to 
these are scarcely visible in the main diffraction 
pattern seen along Z, although there are gaps at 


1we are indebted to photographic technician V. A. Vinogradov for 
processing the color photographs. 
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the appropriate points in the second rings of the 
main and higher-order figures. Figure 3 also shows 
two diffraction points within the first ring of the 
main figure, and also the same within the firstrings 
of the higher-order figures. The distance between 
points in adjacent figures show that they correspond 
to fast shear waves along X and arise from diffrac- 
tion of the beams corresponding to the odd peaks 
produced by the slow shear waves along X, not 
from diffraction of the main beam. 

Figure 3 thus shows that diffraction along Z 
arises from fast shear waves along the X axis in 
quartz only when the slow shear waves along the 
same axis are strongly excited; and even then itoc- 
curs in the beams corresponding to the odd peaks 
produced by the slow shear waves. 


These effects are not explicable by existing ap- 
proximate theories of diffraction at ultrasonic 
waves; they require a more rigorous and complete 
treatment of the orientation effects produced by ex- 
citation of shear waves. 
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The photoelectret charge has been measured as a function of intensity and time of 
illumination during polarization. The lifetime of the residual polarization is discussed. 


Introduction 


Polarization in photoconducting crystals was 
studied by Hilsch and Pohl [1], after it had been 
shown [2] that photopolarization in crystals is as- 
sociated with the accumulation of a space charge 
due (at least in part) to electron trapping. Later 
[3,4] additively colored alkali halides were used in 
measurements of the spectral distribution of the 
depolarization current. These curves had two peaks, 
which correspond to transition into the conduction 
band from the F centers (the main short-wave peak) 
and from the F' centers, the latter associated with 
defects. The polarization during photoconduction 
was due to the transitions from the F centers into 
the conduction band, where the electrons are dis- 
placed by the external field and eventually become 
trapped at F' centers. Photodepolarization is due 
to transition into the conduction band from F and 
F' centers, the electrons then migrating in the field 
of the space charge. Similar deductions have been 
drawn [5] for single crystals of sulfur. The photo- 
polarization of sulfur has also been examined by 
Kurrelmeyer [6] and Vorob'ev [7]; amorphous se- 
lenium shows a similar effect [8]. 

Nadzhakov, in 1937, discovered fixed internal 
photopolarization in polycrystalline sulfur [9] and 
examined the lifetime of the state in the dark. He 
found that a polycrystalline plate retained its po- 
larization unchanged for at least 15 h, in which it 
resembled the electrets discovered by Eguchi in 
1921. The heterocharge on an Eguchi electret is 
due to the temperature dependence of the ion mo- 
bility and dipole relaxation time. 
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Nadzhakov's electrets are evidently produced 
via photoconduction and trapping of displaced 
charges at localized levels. Nadzhakov proposed 
[9] the name photoelectrets, Eguchi electrets then 
being thermoelectrets. The methods devised for 
charge measurement on electrets could then be ap- 
plied to photopolarization. 


Nadzhakov has since made further measure- 
ments of the lifetime of the polarization in poly- 
crystalline sulfur and of the effects of temperature 
on the polarization [10,11]. He showed [10] that 
the purity and thermal history of the material are 
important. The polarization persisted almost unal- 
tered for 37 h in pure polycrystalline sulfur (maker 
Kahlbaum), whereas sulfur of commercial grade 
showed a considerable fall in this time. In addi- 
tion, he found that single crystals showed a more 
rapid fall. Brief illumination during polarization 
revealed that the photopolarization proper is ac- 
companied by dark polarization (which arises when 
a field is applied in the dark); the latter is depend- 
ent on the dark resistance and eventually decays to 
zero, whereas the photopolarization cannot be com- 
pletely eliminated by internal dark conduction. 


A recent paper [12] deals with the persistent 
internal polarization produced by visible light, ul- 
traviolet, x rays, and fast electrons; italso reports 
the relation of charge to polarization conditions and 
seeks to relate the polarization to photoconduction 
and luminescence. However, no data are presented 
on the lifetime of the polarization (electret state 
proper). Another recent paper [13] merely presents 
essentially Nadzhakov's results of 1937. 
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Here we report preliminary results for the 
photoelectret state in single crystals, as the per- 
sistent (or slowly decaying) internal photopolariza- 
tion in a single crystal appears to have a bearing 
on crystal physics generally. It is also of interest 
to compare the lifetime and production conditions 
for single crystals with those for polycrystalline 
aggregates. 

The persistent internal polarization in a poly- 
crystalline material is equivalent to a polar texture 
of symmetry © - m (Shubnikov); Nadzhakov showed 
that such a texture is stable in pure polycrystalline 
sulfur. Residual polarization in a single crystal is 
equivalent to a special polar direction, which in 
general is equivalent to alteration of the initial sym- 
metry to the symmetry of one of the subgroups of 
o-m. Itis of interest to compare the magnitude 
and lifetime of the polarization for nonpolar, polar, 
and special polar directions, particularly as re- 
gards the anisotropy of these for a given crystal. It 
is to be expected that the production conditions and 
magnitude will vary with the direction. The life- 
time of the internal polarization is the lifetime of 
the altered symmetry. 

Here we give results only on the production 
conditions and lifetime for photopolarization in sul- 
fur single crystals. 


Methods 


We used plates about 1 cm? in area and 1-1,5 
mm thick, cut from natural crystals of orthorhom- 
bic sulfur perpendicular to the acute bisector of the 
optic axes; in a few cases we used cubes whose 
faces were perpendicular to the diad axes. The 
specimens were brought to a mirror finish by po- 
lishing on silk wetted with kerosene. 

The measurements were made with the vacu- 
um-tube electrometer described in [14]; the meas- 
uring cell is shown in Fig. 1. The body 1 has a 
hole to admit light, which is closed by the quartz 
plate 2 held by ring 3, which carries a brass grid, 
which acts as one electrode. The second electrode 
4 is of aluminum and is held in the amber insulator 
7 mounted in the top 6, which fits tightly in the body. 
The specimen 5 is gripped between the two elec- 
trodes by the spring 8; it has no guard-ring, because 
surface leakage in sulfur is very slight [7]. The po- 
larizing voltage is applied to the body of the cell 
and hence to the grid electrode; the aluminum elec- 
trode is connected to the grid of the electrometer 
tube. 

The cell was placed in a thermostat whose tem- 
perature was adjustable up to 100°C; this allowed 
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Fig. 1. The measuring cell. 


entry of a beam of light parallel to the polarizing 
field. The intensity at the specimen was measured 
with an exposure meter of known spectral response, 
for which purpose the cell was moved to one side. 

The thermostat and light source (PRK-4 mer- 
cury arc) were mounted on an optical bench. The 
high voltage was provided by an electronically sta- 
bilized supply. All measurements were made in 
terms of the steady-state reading. 

The voltage was usually applied at the same 
time as the illumination, the latter being turned off 
after a set time; the voltage was then removed and 
the grid electrode was grounded. The plate was 
left in the dark for a set time with the electrodes 
joined together; the charge was measured by de- 
polarization in response to fresh illumination (de- 
polarization current as a function of time). The 
polarization (numerically equal to the surface 
charge density) was usually evaluated from the ini- 
tial discharge current. The current could also be 
integrated to give the charge released; the polariz- 
ing current was also recorded as a function of time. 


Results 


Polycrystalline aggregates and single crystals 
were used, the latter 1 em? by 1.5 mm thick (see 
above); polycrystalline specimens of about the same 
size were made by melting these plates, with sub- 
sequent recrystallization by Nadzhakov's method 
(Lol: 

The specimen was polarized under white light 
for 10 min at 200 V; then the electrodes were 
grounded, the specimen being left for a set time in 
the dark. The charge was measured by depolariza- 
tion. 

Figure 2 shows the results as the ratio of the 
initial discharge current it at time t after the po- 
larization to ij, the same immediately after po- 
larization. The single crystal retained hardly any 
polarization after 50 h in the dark, whereas the 
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polycrystalline specimen retained about 40% after 
90 h. 

The effect is to be ascribed to differences in 
the conditions for electron trapping at localized 
levels. The activation energy governs the trapping 
time, and this is clearly much greater for the poly- 
crystalline specimen. The limiting (perhaps finite) 
value of the polarization for the polycrystalline ma- 
terial demands special consideration. 

The charge (polarization) of a single crystal is 
also dependent on the intensity and illumination 
time (Figs. 3-5). 

Figure 3 shows iy as a function of intensity dur- 
ing polarization for 10 min in white light at 200 V 
for a single crystal. The intensity could be varied 
by moving the source and by inserting metal grids 
(neutral filters). The discharge current was meas- 
ured 5 min after the end of the illumination period, 
the depolarizing intensity being always 3 - 107 W 
per cm’. This figure also shows the initial current 
during polarization as a function of intensity; the 


443 


two are linearly related, but the polarization shows 
saturation, which is apparent even at E = 5° ton: 
W/cm’. 

Figure 4 shows the relation of initial discharge 
current to intensity at low intensities; the polariza- 
tion produced at 2.4 : 10-° W/cm? was equivalent to 
ip = 1.3 ° 10-!! A, while the dark polarization pro- 
duced in the same time was equivalent to iy = 0.5 ° 
10“ A. The photopolarization effect thus differs 
appreciably from dark polarization even at this re- 
latively low intensity; but it is of little significance 
to attempt to divide the total polarization into light 
and dark parts as regards lifetime, because the 
over-all polarization of a single crystal eventually 
falls to zero. On the other hand, such a division is 
desirable for polycrystalline material, because the 
dark polarization decays to zero, but the photopo- 
larization is not lost completely [10]. 

Figure 5 shows the photopolarization (in terms 
of ip) as a function of illumination time. The single 
crystal was polarized with white light (3 - 104 W 
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Fig. 2. Residual polarization in: (a) single crystal of sulfur; (b) polycrystalline 
aggregate, 
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Fig. 3, Intensity dependence of: (a) initial depolarization current; (b) initial polarizing current. 
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Fig. 4. Intensity dependence of initial depolari- 
zation current for low intensities during polariza- 


tion, 


t, min 


Fig. 5. Initial depolarization current as a function of illumina- 
tion time during polarization, 


per cm? at 200 V, the latter always applied for 10 
min. Depolarization was by white light at 3 - 105° 
W/cm? after 5 min in the dark. There is clearly a 
saturation effect. 

Illumination for only 1 sec (with a camera 
shutter) gave a polarization corresponding to ig = 
0.8-107!! A, whereas dark polarization for 10 min 
gave only 0.5 - 1071! A: but shorter times (p(s - 
0.2 sec) gave a total polarization nearly equal to 
the dark component, so the photopolarization was 
undetectable in this way. The photopolarization 
produced by such short exposures may [10] be 
measured by leaving the specimen for a sufficient 
time in the dark; the residual effect is the photopo- 
larization alone. This method is not applicable to 
single crystals, because the total polarization falls 
rapidly to zero. 

Comparison of Fig. 3 with Fig. 5 reveals re- 
ciprocity failure for the internal polarization; a 


given Et (in which E is intensity and t is time) pro- 
duces a residual polarization dependent on the ex- 
posure time. 

The saturation effects observed for intensity 
and time are due to filling of the restricted number 
of traps associated with the polarization. The num- 
ber of such traps might be determined from the 
quantum yield if monochromatic light were used. 

We are indebted for discussions to Academi- 
cian A.V. Shubnikov, to G. Nadzhakov (Bulgarian 
Academy of Sciences), and to I.S. Zheludev.Assist- 
ance with the measurements from Yu.N. Martyshev 
and A.I. Delova is also gratefully acknowledged. 
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An optical and x-ray study has been made 


Nomeo 


of kink bands in Tl1Br—TII specimens having 


different crystallographic orientations. Complete agreement has been found between the 
optical pattern and the results of the x-ray study. The gradual nature of kink-band forma- 


tion has been investigated. 
are not always present in a kink band. 


We have found [1] that in the longitudinal bending 
of TIBr—TlI crystals, kink bands are sometimes 
formed in them. The term "kink" is used by us to 
denote a form of deformation resulting in the elbow- 
shaped bending of the specimens, when part of a 
specimen is displaced in relation to the other part, 
and a transition zone is formed between the dis- 
placed parts. We call this transition zone a kink 
band. In another of our investigations [2], a study 
was made of the conditions under which the defor- 
mation of crystals of thallium and cesium halides 
was accompanied by kink-band formation. It was 
found that kink bands were not formed when the crys- 
tals were subjected to uniform tension. They were 
formed only if tension was accompanied by macro- 
scopic bending. In the case of deformation by com- 
pression, the kinks were formed at the ends or in 
the middle part of the specimen if compression was 
accompanied by longitudinal bending. Thus, kink 
formation is always due to the presence of nonuni- 
formly distributed stresses. 

Optical investigations showed that the crystal 
lattice gradually changes its orientation over the 
width of the band connecting the crystal parts which 
have been displaced in relation to each other. The 
gradual change in orientation results from the for- 
mation of two opposing fans of wedge-shaped por- 
tions [1]. 

The present article deals with a more detailed 
study of the structure of this transition region (kink 
band). 
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Indications have been obtained to the effect that slip traces 


Specimens of TIBr—TII crystals, in which a 
sharply pronounced kink band had been formed on 
compression accompanied by longitudinal bending, 
served as objects ofthe investigation. The bands in 
specimens of three crystallographic orientations 
were studied. The principal axis of one specimen 
formed an angle of 26° with the [100] direction (slip 
direction); the principal axis of the second specimen 
coincided exactly with the [100] slip direction; in the 
third specimen, the principal axis was parallel to 
the [111] direction. 


Study of Specimen in which the Com- 
pression Axis Formed an Angle of 26° 
with the [100] Direction 


Photographs in polarized and natural transmit- 
ted light of this specimen are reproduced in [1]. A 
thin section, 0.15-mm thick, parallel to the (110) 
plane, was prepared from this specimen. The thin 
section was made in the following way. A plate 3 
mm thick was cut from the specimen by means ofan 
automatic fret saw (TIBr—TII crystals are difficult 
to dissolve and, therefore, it was not possible to 
cut them by means of a thread moistened with a sol- 
vent). For removing the surface layer, damaged in 
cutting, from one side, this plate was ground (with 
emery) on ground glass for about 1 mm. The sur- 
face of the plate was then polished, and the plate was 
stuck onto glass by the polished side to avoid bend- 
ing when subjecting the other side to a similar treat- 
ment. A number of Laue diffraction patterns were 
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Fig. 1. Laue pattern from non- 
deformed part of specimen. The 
indices of the spots are shown 
(exposure, 4 h). 
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obtained from this thin section, using wide x-ray 
beams (diameter 1.5 mm). The x-ray irradiation 
used was from a tube with a molybdenum target, vol- 
tage57kV andcurrent10mA. Exposure was 4to12h. 

First of all, the specimen was adjusted so that 
the x rays passed strictly along the [110] direction. 
Adjustment was made according to Laue patterns 
from an undeformed portion of the specimen (Fig.1). 
Laue patterns from the kink band were then taken. 
Figure 2 is a diagram showing the structure of part 
of a kink band, drawn from a photomicrograph of 
the thin section. The circles in the diagram mark 
the places for which Laue patterns were photo- 
graphed. The width of the band was such that the 
x-ray beam completely covered it. Laue patterns of 
the kink band showed splitting of the spots along the 
circumference into a number of reflections (Fig. 3). 
Radial broadening of the spots was not observed. 
This is evidence of the fact that the x-ray beam 
passed strictly along the axis of rotation of the lat- 
tice. Since the x-ray beam passed along the [110] 
direction, the axis of rotation of the lattice was the 
[110] direction. 

In the microscopic study of this specimen, it 
was observed that all the wedge-shaped regions 
were made up of slip traces, ! changing their direc- 


Fig. 2, Diagram of kink-band structure in a T1Br—TI 
crystal, drawn from a photomicrograph. The circles 1 
and 2 mark the places subjected to x-ray study with a 
wide x-ray beam (diameter 1.5 mm); the circles 3 and 
4 denote places from which Laue patterns were ob- 
tained with an x-ray beam 0,1 mm in diameter. 
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Fig. 3. Laue pattern from 
kink band (exposure, 12 h). 


Fig. 4. Photograph in natural transmitted 
light of part of kink band in a T1Br— TI 
crystal with impact figures formed in it. 
Demonstrates lattice rotation in the kink 
band, 


c d 


Fig. 5. Form of Laue pattern spots taken 
(a) from entire kink band (part 1, Fig. 2), 
set of 34° angles; (b) from edge of kink 
band (part 2, Fig. 2, exposure 10 h), set 
of 10° angles; (c) from boundary of two 
wedge-shaped regions (part 3, Fig. 2), 
and, (d) from an individual wedge- 
shaped region (part 4, Fig. 2). 


1By the term “slip traces" we understand bands visible in thallium 
halide crystals in transmitted polarized light, owing to double re- 
fraction caused by local internal stresses, and in natural transmitted 
light, owing to local variation of the index of refraction. 
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Fig. 6. Photomicrographs of two cuts of T1Br—TII crystal, compressed in 
the [100] direction. Impact figures produced inside and outside the kink 
band: (a) cut along the (110) plane; (b) cut along the (110) plane. 


tion on passing from one edge to another (Fig. 4). 
The plane of the thin section was perpendicular to 
the slip planes, so that judging by the arrangement 
of the slip traces, the axis of rotation of the lat- 
tice was situated in a slip plane, and was perpendi- 
cular to the slip direction. 

Figure 5 shows individual x-ray pattern spots 
taken from two portions of the specimen. The spot 
in Fig. 5a corresponds to the case where the beam 
entirely covered the whole of the band (part 1 in 
Fig. 2). The complete x-ray pattern for this case 
is shown in Fig. 3. According to measurements 
made on a gnomonic projection of this Laue pattern, 
the aggregate of the reflections corresponding to a 
spot on the x-ray pattern (Fig.1) corresponded to 
the aggregate of lattice rotations in the kink as a 
whole through an angle of 34°. When the second x- 
ray photograph was taken, one of the spots of which 
is reproduced in Fig. 5b (part 2 in Fig. 2), the x- 
ray beam was displaced from the center to close to 
the edge of the specimen, and covered only a few 
wedge-shaped regions. In this case, the spot had a 
spread of 10° and consisted of five separate reflec- 
tions. Each pair of adjacent reflections corre- 
sponded to two wedge-shaped regions rotated 
through approximately 2° in relation to each other. 

According to optical measurements, the lattice 
in adjacent wedge-shaped regions of this crystal 
was also rotated through an angle of 2° [1]. From 
a comparison of the results of the optical and x-ray 
measurements, it may be said that in this case the 
x-ray beam covered approximately five wedge- 
shaped regions. This conclusion is in agreement 
with the number of wedges covered by the beam 


(part 2) in Fig. 2. The degree of broadening of 
these spots was ascertained by superimposing two 
x-ray patterns: one taken from an undeformed part 
of the specimen (Fig.1), and the other from a kink 
band (Fig. 3). 


Furthermore, to obtain information on the 
state of the lattice inside the wedges and on their 
boundary, x-ray patterns of these parts were ob- 
tained by means of a narrow x-ray beam 0.1 mm 
in diameter. The narrow beam was produced by 
placing a lead diaphragm 0.1 mm in diameter at the 
exit of the x rays from the collimator. In taking 
these x-ray patterns, a tube with a tungsten target 
was used (voltage 57 kV, current 13 mA, exposure 
40 to 60 h). 


Figure 5 (c and d) shows spots taken from two 
of the Laue patterns. The spot in Fig. 5c corre- 
sponds to the boundary between two wedges (the 
beam covered parts of two adjacent wedges). Itis 
quite clear that this spot can be divided into two re- 
flections, with an angular spacing between them of 
~2°, Figure 5d shows a spot of an x-ray pattern 
taken from inside a wedge-shaped region. This spot 
is quite unsplit, and near it is a very faint radial 
diffuse tail, invisible in the reproduction. Similar 
faint tails are to be found near spots in x-ray pat- 
terns taken with a narrow beam, remote from the 
kink band in the undeformed region of the specimen. 


On the basis of the x-ray patterns reproduced 
(in Figs. 5a-d), it may be concluded that the split- 
ting of the spots is due to the presence of the aggre- 
gate of rotated regions, i.e., to macroscopic dis- 
tortion of the crystal. The slip traces present in 
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the interior of the wedges do not cause asterism of 
Laue spots. 

Lattice rotations in the kink band can be demon- 
strated very graphically by means of impact figures 
[3].. Figure 4 shows a photomicrograph of a portion 
of a kink band. The wedge-shaped parts and the 
displacements inside them are clearly to be seen 
(photography in nonpolarized light). In addition, two 
large impact figures, produced by means of a 
needle, are visible. One of these figures has been 
applied outside the kink, and the other inside it.The 
long axis of the figures coincides with the [110] di- 
rection. From the mutual arrangement of the im- 
pact figures, it may be seen that the lattice in the 
kink band has been rotated with reference to the lat- 
tice of the original crystal, and it is even possible 
to make an approximate estimate of the angle of ro- 
tation. 


Study of a Specimen with the Axis 
Parallel to the [110] Direction 


Figure 6 shows lattice rotations in a kink band 
by means of impact figures in another T1Br—TII 
specimen, compressed along the principal axis, 
parallel to the direction of the cube edge. Both 
photographs were taken in nonpolarized light. The 
plane of cut of the plate shown in Fig. 6a was per- 
pendicular to the axis of rotation of the crystal lat- 
tice, while the plane of cut of the plate shown in 
Fig. 6b was parallel to that axis and perpendicular 
to the first cut. Both plates were approximately 
parallel to planes of the rhombododecahedron. In 
the first figure, the impact figures both inside and 
outside the kink band (i.e., in the latter case, in 
the nondeformed region) have the same form but 
are rotated in relation to each other through an 
angle of 28°. The axis of rotation is normal to the 
polished surface of the cut, i.e., it is the [110] di- 
rection. In the second cut (Fig. 6b), the axis of ro- 
tation of the lattice lies in the plane of the cut, and 
impact figures corresponding to the kink band and 
to the nondeformed part of the crystal have differ- 
ent symmetries; whereas, outside the band, the im- 
pact figure has a two-fold axis of symmetry (corre- 
sponding to the crystal symmetry in the [110] direc- 
tion), inside the band, the impact figure has ap- 
proximately a three-fold axis of symmetry and cor- 
responds to a plane close to [111]. Inside the kink 
band, the lattice has been rotated in such a manner 
that in place of the rhombododecahedron plane there 
has emerged a plane close to the octahedron plane. 

The lattice in the kink band is not in the twin- 
ning position in relation to the original crystal, 
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since the interface between them is not a plane of 
symmetry of the conjugate lattices. This conclusion 
is based on x-ray patterns taken perpendicular to 
the cut, shown in Fig. 6a, from a part containing a 
kink-band boundary. The contact boundary was 
found to have the indices (241) in relation to the ini- 
tial lattice, and the indices (17.17.10) in relation to 
the rotated part of the crystal (kink band). 

Unlike the specimen shown in Fig. 4, in the 
present case there were no slip traces in the kink 
band. Subsequent experiments with ten specimens 
with this orientation confirmed the absence of slip 
in the band. 


Study of Specimen with Its Axis 
Parallel to the [111] Direction 


A general view of this specimen, in which the 
principal axis coincided with the [111] direction, in 
polarized transmitted light is shown in Fig. 7. On 
compression, the specimen underwent a smooth 
longitudinal bend, and assumed a form reminiscent 
of a cosine curve with a 3/4 period. On compres- 
sion, two kinds of bands were formed, one (1-1) of 
which passed through the entire cross section of the 
Specimen, while the other (2-2) failed to reach the 
other side of the specimen. Inside both kink bands 
there were many cracks, but in the second band 
(incompleted kink) cracks were considerably fewer. 


Fig. 7. Photograph in polarized light of 
T1Br—TII specimen compressed along an 
axis parallel to[111]. Asa result of the 
compression, two kink bands were formed, 
one (the bottom one) of which failed to 
pass through the entire specimen. 
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On examination of the structure of the incomplete 
kink, an impression was gained of the gradual na- 
ture of its development. Figure 8 shows Laue spots 
taken from three different parts of the incomplete 
kink. During exposure, the x-ray beam passed 
along the [110] direction, coinciding with the direc- 
tion of the axis of rotation of the lattice in the band, 
which is shown by splitting of the spots only on the 
circumference. The x-ray pattern spot (Fig. 8a) 
taken from the end of the band, corresponds to the 
initial stage of kink formation. This Laue spot 
shows that in the initial stage of deformation there 
were no considerable lattice rotations, the lattice 
distortions were quite small and produced only 
slight broadening of the Laue spots in the radial di- 
rection. The middle part along the length of the band 
corresponds to more intense deformation when con- 


siderable lattice rotations were already taking place. 


In the Laue pattern (Fig. 8b), taken from this part 
of the kink band, circumferential splitting of the 
spots into three spots with a final angle of rotation 
of 9° can be observed. The intensity of the addition- 
al spots, however, is very small. This shows that 
an insignificant part of the crystal volume was in 
the rotated position. In the initial part of the band 
of the incomplete kink, situated close to a lateral 
surface of the specimen, the deformation of the 
crystal was greater. The considerable intensity of 
the additional spots in the x-ray photograph shown 
in Fig. 8c indicates that a larger volume of thecrys- 
tal was in the rotated condition. In this case, the 
number of additional spots increased, as also did 
the range of the angles of rotation (18°). This lat- 
ter fact shows that at large deformations not only 
does the volume of rotated material increase, but 
also the number of rotated regions, and, conse- 
quently, also the set of discrete orientations of the 
crystal lattice inside the kink band. 

Thus, a study of the structure of an incom- 
plete kink band has shown that a kink band develops 
by spreading successively into newer and newer re- 
gions. At first, lattice rotation occurs in a small 
region of the crystal, and the lattice rotates through 
a small angle relative to the initial position. Fur- 
cher development of the band occurs by increase of 
the angle of rotation in the previously rotated re- 
gion and, in addition, rotation takes place in adja- 
cent parts situated on either side of the originally 
formed band. 

In both kink bands of this specimen, as in the 
specimen compressed along the [100] direction 
(Fig. 6), no slip traces were visible. 

We shall examine the arrangement of the kink 
bands in relation to the crystallographic axes in the 
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| Fig. 8. Laue spots correspond - 
ing to three places in an in- 
complete kink: (a) end; (b) 
middle; (c) beginning. 


three specimens studied. In the first of them (the 
principal axis of the specimen made an angle of 26° 
with the [100] direction), the kink band was situated 
at an angle of 48° to the principal axis, or at an 
angle of 106° to the [100] direction. In the second 
specimen (the principal axis was parallel to the 
[100] direction), the kink band was formed at an 
angle of 74° to the principal axis, coinciding with 
the slip direction. In the third specimen, the kink 
band made an angle of 31° with the [111] axis,or 85° 
with the [100] slip direction. On a basis of the com- 
parison of these data, it may be concluded that a 
kink band is always formed approximately perpen- 
dicularly to the planes and direction of slip, but not 
in the direction of action of the maximum cleavage 
stresses, as might be expected. 


Discussion of the Results 


An x-ray study of kink bands in TIBr—TII crys- 
tals has shown that a kink band consists of an aggre- 
gate of wedge-shaped regions rotated in relation to 
each other. 

The rotation of the lattice in individual (adja- 
cent) wedge-shaped regions does not occur gradual- 
ly but abruptly, with a jump. The results of the x- 
ray study are in agreement with an optical study and 
with a study of the rotation by means of impact 
figures. Such abrupt lattice rotations in the case of 
kinking have not been observed in other crystals, 
for example in zinc [4], cadmium [5], CsI [6]. The 
authors of articles on the investigation of these crys- 
tals consider that there is smooth bending of slip 
planes in kink bands. This view, however, is not 
shared by other investigators. Gilman [7] affirms 
that in kink bands of zinc crystals, fracture of slip 
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planes occurs, not their bending. Moore has con- 
firmed this experimentally [8]. He measured the 
slope of the basal planes in relation to the original 
crystal in a twin and in the accommodation bands 
adjacent to the twin. He found that the accommoda- 
tion bands, representing transition zone, like kink 
bands, consisted of one or more undistorted re- 
gions, rotated through angles of 10'-30' in relation 
to each other and to the original crystal. 

We consider that on deformation, kink forma- 
tion does not occur by bending of the lattice but by 
discrete rotations. The magnitude of an elementary 
rotation depends on the plastic properties of the in- 
vestigated crystal, and also on the experimental 
conditions (crystallographic orientation of the spe- 
cimen, rate of deformation). 

It should be noted, however, that in investigat- 
ing the structure of kinks, it is necessary to take 
into account the time elapsing between the moment 
of producing the kink and its examination. Grounds 
for such a conclusion may be provided by Honey- 
combe's experiments [9] in explanation of the influ- 
ence of recovery on kink-band structure in cadmium 
and aluminum. From Laue patterns of kink bands 
in cadmium at different stages of recovery, Honey- 
combe discovered that closely situated faint spots 
in the course of time merged into a number of 
strong spots situated fairly remote from each other. 
In aluminum single crystals, this phenomenon was 
not observed at room temperature, but heating to 
600°C was necessary for merging of the spots. 
After high-temperature annealing of deformed spe- 
cimens of cadmium, intensive recrystallization oc- 
curred and the entire kink band was transformed in- 
to two grains. 

The specimens we studied had been kept for a 
long time (several months) before x-ray examina- 
tion. It is possible that during this period their re- 
covery occurred, and small slightly rotated wedge- 
shaped regions united to form larger portions, more 
considerably rotated in relation to each other, there- 
by increasing somewhat the abruptness of the rota- 
tions. 

Furthermore, in all the foreign articles on the 
study of kink bands in metals, it is pointed out that 
the kink boundary bisects the angle between the slip 
traces situated on either side of this boundary. 
Orowan [5], in outlining kink formation, states that 
the lattices in such regions are not in a mutually 
twinning position, although the plane of contact is 
common to both of them. In contradistinction to 
this, however, Kolontsova etal. found that in CsI 
crystals, in one part of a kink band, the lattice of 
one of the wedge-shaped regions was rotated in re- 
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lation. to the original crystal in accordance with the 
twinning law. In the twinned components, the (113) 
plane was common to both. 

Stepanov [10] points out that in the flaking of 
rocksalt crystals, the boundaries of the flakes are 
not crystallographically equivalent in adjacent ro- 
tated regions (flakes). Our investigations of a kink 
band in TIBr —TII [1] gave a similar result. Re- 
peated measurement of the angle between the boun- 
daries of the wedge-shaped regions and the slip 
traces inside them in one of the specimens showed 
that one of the boundaries of these regions was per- 
pendicular to the displacements, while the other 
made with the first boundary an angle of 2°, being 
perpendicular to the displacements present in the 
adjacent wedge-shaped region. 

It may be concluded from this that the rotation 
of the lattice may not only be symmetrical, as oc- 
curs in a CsI crystal, but also nonsymmetrical, 
when the adjacent rotated regions do not represent 
a twin. 

The presence of slip traces in the wedge-shaped 
parts of a kink band may lead one to conjecture that 
nonsymmetrical rotation is performed by reason of 
slip accompanied by rotation of slip planes around a 
direction perpendicular to the direction of slip and 
situated in theslip plane. 

Currently, therefore, the following proposals 
exist concerning the mechanism of kink-band forma- 
tion: 

1. Slip accompanied by bending of the slip 
planes; 

2. Slip accompanied by rotation of the slip 
planes; 

3. Twinning on an irrational plane [11], i.e., 
symmetrical rotation; 

4. Nonsymmetrical rotation. 

We do not support the first of these four pro- 
posals. The participation of twinning along irra- 
tional faces [6] and nonsymmetrical reorientation of 
the lattice in kink-band formation we regard as de- 
monstrated. Twinning on irrational faces and non- 
symmetrical reorientation of the lattice represent 
different stages of the same phenomenon. 

Figure 4 shows that slip traces exist inside 
each wedge. We have observed a similar phenom- 
enon in the compression of cesium bromide and 
iodide crystals. In kink bands formed in these 
crystals, in all the orientations we have investi- 
gated, slip traces (the same slip elements as in the 
case of TlBr —TII) were always present. 

It should be noted, however, that we have ob- 
served slip traces in specimens examined either 
after completion of the deformation process (Fig. 4), 
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during the deformation process by filming at a rate 
of 24 frames per second, or when kink bands were 
fully formed in the time between two frames (in a 
time <0.04 sec) [2]. On the basis of these data, we 
were unable to affirm that lattice reorientation in 
the wedge-shaped regions and slip necessarily oc- 
cur simultaneously. Slip could develop even after 
reorientation of the lattice. The photograph in Fig. 
6 and the additional investigations of nine speci- 
mens (the axis of compression of six specimens 
was parallel to [100], the axis of compression of 
three specimens was at a small angle to the [100] 
direction), in which we also failed to observe slip 
traces on examination in polarized light, support 
the possibility of lattice rotation occurring without 
the participation of slip. 


Summary 


1. Kink bands in three specimens of T1Br—TII 
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having different crystallographic orientations have 
been studied. 

2. Complete correspondence between the op- 
tical pattern of the kink-band structure and the form 
of the Laue spots has been shown; kink bands are an 
aggregate of regions rotated in relation to each 
other. 

3. The axis of rotation is the [110] direction. 

4, Kink bands are developed approximately 
perpendicularly to the slip planes. 

5. The gradual nature of kink-band formation 
has been studied in the case of an incompleted kink. 

6. It is shown that kink formation is not al- 
ways accompanied by slip. Kinks may be formed 
by means of nonsymmetrical lattice rotations. 

The radiography was carried out in the X-Ray 
Laboratory of the Institute of Crystallography. We 
are extremely grateful to Academician N.V. Belov 
for the facility to work in the laboratory supervised 
by him. 
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The phase diagram of 99.995% Bi has been examined at pressures up to 30,000 kg/cm? 
and at temperatures from normal up to the melting points of the phases. The 6 phase 
melts at 184°C between 17,300 and 22,400 kg/ em’, so 6-Bi melts without change of 
volume, whereas the melting point of y-Bi increases with pressure, with a mean coef- 
ficient of 0.0072 deg-cm?/ keg. The af and 6=y transitions are very fast but show 
substantial hysteresis at all the temperatures used. The ratio of the heats of transi- 


tion and the volume changes in these transitions have been determined by a thermo- 


metric method. 


Bridgman [1] first detected the polymorphic 
transitions in bismuth under hydrostatic pressure; 
he used the displacement of the plunger. The first 
(monoclinic) modification is stable up to 25,900 kg 
per cm? at 20°C; the second is stable over the range 
25,900 to 27,760 kg/cm’; and the third is stable at 
and above 27,760 kg/cm’. The lattice structures 
of the latter two modifications are unknown. These 
three modifications may be denoted in the standard 
way as the a, 6 , andy modifications. Bridgman 
recorded the melting point of @-Bi as a function of 
pressure, the @-f6 and 6-y equilibrium curves 
over the range —50 to 170°C, the changes in spe- 
cific volume in the a = 8 and 6 =y transitions, 
the heats of transition, and the coordinates of the 
triple points a—B-liq, B-y-—liq, a—fB-y. 

The transition pressures of bismuth at room 
temperature are presently the principal reference 
points in the calibration of manganin resistance 
gauges. Bridgman determined the transition pa- 
rameters of bismuth several times [2-5], but only 
at room temperature. The transition points were 
detected from the displacement of the plunger and 
from the resistance, but the latter does not yield 
important parameters such as the heat of transi- 
tion or the volume change, while the piston dis- 
placement does not yield the volume change for 
each transition separately, »ecause the two transi- 
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tion points are close together. In his later papers 
[2,3] he gives the sums of the volume changes for 
the two. Other disadvantages of his methods are 
that they do not give access to temperatures above 
200-250°C, nor can the transition rates be deduced. 

The liquidus line for bismuth was followed up 
to 22,200 kg/cm? in 1953 at this institute by differ- 
ential thermal methods [6]. 

We have now made a detailed study of the phase 
diagram at pressures up to 30,000 kg/ em? on bis- 
muth of 99.995% purity!; the heats and volume 
changes in the a = 6 and 6 = y transitions have al- 
so been determined. The hydrostatic pressure was 
provided by a mixture of isopentane with n-pentane 
in a pressure amplifier [7]; a manganin device 
reading to +100 kg/cm? being used in conjunction 
with an iron—nichrome thermocouple reading to 
+1.5°. Full details of the methods have been given 
elsewhere [8]. 

The heating and cooling rates were about 1 deg 
per sec. The transition points were readily de- 
tected with the differential and single thermo- 
couples, on account of the fairly high rates of 
transition and latent heats. Both transitions show 
appreciable thermal hysteresis; Table 1 gives data 


1The main impurities in the Bi were: 0.003% Pb, 0.001% Fe, and 
0.0003% Hg. 
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Fig. 1. Phase diagram for bismuth. O— 
Our results; x — Bridgman's results, 


on the @ ~8 and 6 > transitions at various pres- 
sures. 

The hysteresis in the @ = 8 transition in- 
creases from 3° for Tay = 163. 5° to 15° at Tay = 
65.5°, in which Tg, is the average of the tempera- 
tures for the a > B and 6— dc transitions at that 
pressure. 

Figure 1 shows the P—T phase diagram de- 
rived from the data results, the points being Tay, 
and each being the mean of at least six measure- 
ments. The new melting curve of @-Bi agrees 
closely with the older one, as do the @-f and B-y 
transition curves. The melting point of 6-Bi is in- 
dependent of pressure within the errors of experi- 
ment, so 6-Bi melts without change of volume. The 
melting point of y-Bi rises from 184°C at 22,400 
kg/cm? to 239° at 30,000 kg/cm?, the mean slope 
being 0.0072 deg-cm?/ kg. 

The latent heats of the transitions were deter - 
mined as follows. The bismuth specimen 1 (weight 
about 6 g) was placed in the glass tube 2 (Fig. 2) in- 
side the container 3, which was closed by the plug 
6. The specimen and container had recesses to 
take the junctions 4 and 5 of the differential thermo- 
couple. The container was placed in the pressure 
amplifier 7, and the pressure was raised to 25,000 
kg/cm? (below the first transition point). It was 
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Fig. 2. Method of determining heats 
of transition in the pressure vessel. 


then raised from 25,000 to 28,000 kg/cm? over 

of the apparatus being kept constant at 20°C. The 
temperature of the specimen was recorded with a 
Kurnakov system. The ratios of the latent heats of 
the a =f and 6B = y transitions were deduced from 
the ratio of the areas on the dta diagram. The 
heats and volume changes were deduced from four 
equations in four unknowns: 


fee AV, 
ahs Ome 
4 far AV, 
9 ear ca Om 
AV, + AV, = — 0.0085 cm*/g [2,3] 
Q; : Qs — ts 


Here, AV, is the change in specific volume in the 
af transition, AV, is the same for the B = y 
transition, (dT/dp),; and (dT/dp), are the slopes of 
the a-6 and 8-y transition curves at 20°C, and 
Q,:Q. = nis the ratio of the heats of transition. 
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Fig. 3. Dta curves for the a= 68 and 8 = y transitions, 
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TABLE 2 
«-68—liq | @-—yY-- liq a — f Y 
Q, Q2, AV), AV», oe aes “¢ 
cal/g \cal/g | cm*/g | cm3/g Bees S cae 2 re 
ae | Sila Pl a le P| & 
Our data Lain —0.71/—0.0055 |—0.0030 | 17300 | 184 | 22400 | 184 | 32200 |--108 
Bridgman [1] i .62/—0.68]/—0.0047 |—0.0035 | 17300 | 183 | 22400 | 185 | 32300 |—110 


In all,we measured the heats six times; in four 
cases the specimen was held in a lucite cylinder, 
and in two in a copper cylinder. Figure 3 shows a 
typical recording made in one run, in which A and 
B denote, respectively, thea—~f and 6 > y transi- 
tions (pressure increasing), while a and b denote 
the reverse transitions (pressure decreasing). 
Parts 1-2 correspond to absorption (release) of 
heat, while parts 2-3 correspond to equilibration 
after the transition. The shape of part 1-2 is 
governed by the rate of change of pressure. The 
density of 6-Bi is higher than that of @-Bi, so the 
volume is lower. An applied pressure, only slightly 
in excess of the equilibrium value ,causes the 
volume to decrease until the equilibrium pressure 
is reached, at which point the transition stops. The 
a@ and 6 forms then coexist. Complete transition 
from @-Bi to the 6 form requires continuing in- 
crease in pressure. Part 1-2 (a-f transition) is 
of low slope when the rate of change of pressure is 
small; the transition occurs at a fixed pressure. 
High rates cause part 1-2 to be nearly vertical; the 
transition occurs in less than 1 sec. The exact 
time can be evaluated from the recording speeds 
and the known lag in the system. The slope of part 
2-3 is largely independent of the rate of change. 
These remarks apply also to the B > y, y > B, 
and B > co transitions. 

The average Q,; :Q, is 2.51, with a maximum 
range of +2%, in spite of the marked variation in 
conditions (container materials and rate of pres- 
sure variation). Table 2 gives our results, and al- 
so Bridgman's. 

Both transitions at 20°C involve considerable 
pressure hysteresis: some 1000 kg/cm? for a = 6 
and some 900 kg/cm? for 8 = y. This hysteresis 
is not dependent on the rate of change. The transi-~ 
tion point has been taken as the mean of the pres- 
sures for the forward and reverse transitions. This 
hysteresis has not been reported before; it remains 
an open question whether the effect is inherent in 
pure bismuth or is due to impurities. 

The heats of the a = £6 and 6 = y transitions 
may be used to deduce similar heats for other sub- 
stances at high pressures. 


Conclusions 


1. Dta has been applied over the range 1 to 
30,000 kg/cm? and temperatures from normal up 
to the melting points of the phases. 

2. The melting point as a function of pressure 
has been determined for y -bismuth; the mean coef- 
ficient is 0.0072 deg-cm?/ kg. 

3. The ratio of the latent heats of the a= 6 
andB = y transitions at 20°C has been deduced, to- 
gether with the changes in specific volume. 

4, There is appreciable hysteresis in thea = 
6 and 6 =—y transitions at heating and cooling rates 
around 1 deg/ sec; the rates of both transitions are 
high at the temperatures used. 
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Cinemicrography has been applied to unstable dendritic crystals; an explanation is offered 


for the observed decomposition and conversion to a granular structure. 


The time for de- 


tachment of a branch has been measured (to about 20%) as a function of the initial thick- 


ness at the point of detachment. 


Dendritic crystals are unstable and eventually 
produce a system of isometric grains, but the pro- 
cess can be very slow. The details are also not 
visible to the unaided eye. 

Time-delay photomicrography has been applied 
to the stages from dendrites to rounded grains in 
ammonium chloride at constant temperature; the 
process is driven by the tendency of the surface 
free energy to become minimal. The experiments 
were essentially a repetition of those of Saratovkin 
[1] with NH,Cl between glass slides. Aqueous 
NH,Cl was chosen because it readily provides the 
high supersaturations needed for the growth of den- 
dritic crystals. 

A drop of the saturated solution was placed be- 
tween two glass slides in an air thermostat at 70°C. 
The excess of solution was squeezed out, and the 
specimen was transferred to the microscope stage, 
where it cooled to room temperature. Water loss 
was suppressed with a layer of grease. The cool- 
ing was accelerated (to produce high supersatura- 
tion) by directing a jet of air onto the specimen. 
Dendritic crystals usually grew inwards from the 
edges; the camera was Started at this point, the ini- 
tial stage being recorded at 24 frames per second. 
At the end of this stage the rate was reduced to lor 
2 frames per second (early period of conversion of 
the dendrites), with the intervals increasing gradu- 
ally toward the end of the experiment (to intervals 
of 1.5-2 min at 5-7 h). The temperature was kept 
constant. 

The conversion of the dendrites began very 
soon after they had grown, at which point the solu- 
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tion was no longer supersaturated, so the solidwas 
in equilibrium with the solution, speaking general- 
ly; but material is taken up from one part of a den- 
drite and is deposited on other parts, on accountof 
the differences in equilibrium concentration as be- 
tween concave, flat, and convex areas (Ostwald- 
Freundlich formula, assuming the surface tension 
to be constant). 

The surface tension is a function of direction 
for a crystal, as is clear from the frames shown 
in Fig. 1 for a closed drop of solution (see [2, 3] 
for description of the method of making these closed 
drops). The migration of material continues until 
20;/1; becomes the same for all faces present in 
the equilibrium form of the crystal [4] (oj is sur- 
face tension, / is central distance). Here the vari- 
ation of oj between faces is important, as well as 
the shape and size of the interfaces. The final 
frame shows the crystal as a body of nearly con- 
stant curvature, the deviations being due to the 
anisotropy in the surface tension. The solubility of 
this crystal is the same at all points. 

The transformation of dendrites can be ex- 
plained on this basis. It is clear from the film that 
material is transported from convex parts to con- 
cave ones; but these parts are often widely sepa- 
rated, so there is an additional tendency for ma- 
terial to migrate from adjacent plane areas to con- 
cave ones. 

Figure 2 illustrates the transformation of flat 
dendrites of thickness only a few times the wave- 
length of light. The minor branches appear as_ 
though withdrawn into the trunk. Rather thicker 
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Fig. 1. Transformation of dendritic ammonium chloride crys- 
tals in a closed drop of solution, x 525. Time intervals; a-b, 
10 min; a-c, 1h; a-d, 2h 20 min; a-e, 4 h 30 min; a-f, 36h. 


dendrites in addition tend to break up, a process 
that would, at first sight, seem to contradict the 
Gibbs-Curie-Wulff law. However, the migration 
from plane and convex faces to concave ones ex- 
plains the effect, although one might initially sup- 
pose that material from the base of a branch would 
be deposited elsewhere, because the largest nega- 
tive curvature occurs at the base of a branch (Fig. 
3), so this ought to be the site of deposition. The 
negative curvature seen in the photographs is ac- 
companied by a large positive curvature in a sec- 
tion parallel to the base of the dendrite, and the 
latter causes the base to dissolve. A closely 
spaced array of secondary branches causes the 
process to accelerate, on account of transfer of 
material from adjacent parts to the concave area 
at the base. This transport continues until the cur- 
vature at that site is largely eliminated. Even be- 
fore this has been completed, the base has some- 
times been so extensively attacked from both sides 
that the secondary branch is detachedfrom the main 
one. This detachment occurs the sooner the thinner 
the secondary branch, as Chernov has shown [5] in 
terms of a relation of detachment time to thickness. 
This time can be measured quite closely (to 
about 20%) with the microscope (error of measure- 
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Fig. 2, Conversion of thin ammoni- 
um chloride dendrites, x 675. Time 
intervals: a-b, 4h; a-c, 30h. 


ment 2 sec). For this purpose we used 30 second- 
aries on a single dendrite, which gave the relation 
of time to thickness. This agreed reasonably well 
with Chernov's theoretical curve. 

This detachment leaves a sharp spike at the 
base (Fig. 3f), which has a large positive curva- 
ture, so it is rapidly lost. The detached branch al- 
so becomes more nearly isometric; material is de- 
posited on it, although it everywhere has positive 
curvature, especially at the ends, where it dis- 
solves. Some of the dissolved material is deposited 
on the sides. A detached branch separated from the 
base by only a few diameters gives rise to appre- 
ciable deposition even on the flat parts of the base, 
even though the detached part has not yet become 
isometric. This transfer is clear from the reduc- 
tion in area, although the part may increase some- 
what in width. However, even an isometric dendrite 
is unstable in the presence of a nearby base; its 
material is gradually transferred to the latter. 

These transfer processes continue until all the 
crystals have the same effective curvature. The 
rate of conversion to equilibrium (isometric) form 
is dependent on the separation and disposition of 
the surfaces that differ in curvature. Convection 
currents here play no part, and the processes are 
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Fig. 3. Transformation and decomposition of dendritic NH,Cl crystals, 
x 975. Time intervals: a-b, 6 sec; a-c, 15 sec; a-d, 40 sec; a-e, 50 
sec; a-f, 1 min 25 sec; a-g, 4 min; a-h, 5h. 


mainly ones of diffusion; surface diffusion may 
play some part with isolated crystals. 

The conversion rate is also dependent on the 
temperature; a stage can be frozen by appropriate 
cooling, as in many metallurgical techniques. Dif- 
fusion still occurs in these frozen structures, but 
only at extremely low rates. 

All the secondary branches (whether detached 
or not) are lost by the time the conversion is fairly 
well advanced; the specimen consists of interwoven 
main trunks from different dendrites. The process 
does not stop at this stage, for the trunks eventual- 
ly develop necks and break up, the fragments be- 
coming rounded grains. The larger grains then 
grow at the expense of the small ones. The speci- 


men is then a set of randomly oriented grains in no 
way resembling the original dendritic structure. 

Metal castings cannot be observed in this way. 
Very often the processes do not go to completion, a 
frozen structure remaining. Extremely varied poly- 
granular structures may be produced, on account 
of the anisotropy in crystal structure, surface ten- 
sion, impurity distribution, etc. 

We are indebted to Professor G.G. Lemmlein 
for direction in this work. 
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Surface-tension effects are considered in relation to the breakup of elongated crystals 
and inclusions; the conversion time is evaluated as a function of several factors and is 
compared with the experimental evidence. 


Introduction 


In 1951-1956, Lemmlein and Kliya reported 
[1-3] experiments on the effects of surface tension 
on interfaces between saturated aqueous solutions 
and crystals of sodium nitrate, ammonium chloride, 
and silver nitrate. An effect dependent on minimum 
surface energy was observed: liquid inclusions and 
branches of dendrites (both effectively long cylin- 
ders) spontaneously broke up to products that then 
took up equilibrium (nearly spherical) shapes. Fig- 
ure 1 of [8] gives time-lapse photomicrographs that 
illustrate this. 

Here we discuss these effects in the isotropic 
approximation, which gives an estimate of the de- 
composition time as a function of surface tension, 
temperature, concentration, etc. 

A dendrite is essentially equivalent to an inclu- 
sion, so we consider the latter first. 

An elongated cylinder is not the equilibrium 
shape, so a system containing such an inclusion 
will tend to alter the shape of the inclusion at con- 
stant volume. A sphere has the least surface ener- 
gy in the isotropic approximation, but the forma- 
tion of two or more spherical inclusions will also 
reduce the surface energy (the number actually 
formed is dependent on the ratio of length to diam- 
eter). This in no way implies that all of these cases, 
or even the most favorable, will actually occur; the 
detailed mechanism of the transformation is also 
decisive. The only mechanism available here is dif- 
fusion transport within the inclusion (the crystal 
has infinite viscosity and does not allow of hydro- 


dynamic transformation); this arises from differ- 
ences in the concentrations of saturated solutions 
above surfaces differing in curvature. 

Consider the inclusion of Fig. 1 (ratio of 
length 1 to diameter 2a small). The concentration 
at A is less than that at B, and the diffusing par- 
ticles travel as shown by the arrows, The process 
tends to shorten the inclusion while widening it at 
the center, so the shape tends to come to equilib- 
rium. Now let 1 >> 2a. It can be shown that the 
concentration gradients are appreciable only over 
distances of about 2a near the ends, so that the latter 
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have no effect on the center for 1 >>2a. This is 
why a long inclusion splits up into several spherical 
inclusions. 

This may be considered more closely. Let us 
suppose that section OO' (Fig. 2) is slightly nar- 
rower than adjacent parts (e.g., because it was 
formed that way initially), the radius at that point 
being r. If this change of shape tends to reduce the 
surface area, then it will give rise to a flux of ma- 
terial that tends to reduce r further. This will con- 
tinue until the inclusion breaks into two, or, alter- 
natively, further reduction in r does not reduce the 
energy, e.g., because the ends are approached, 
with a resulting increase to radius a. This rather 
complex problem will not be discussed in detail; a 
simpler approximation will be used. 


The mean curvature near OO' is 


24 aaa 


by] = 


4 1 4 
pe Vas erg (1) 


Thomson's formula gives us that the concentration 
of the saturated solution near OO' will be less than 
that in adjacent areas; diffusion then continues un- 
til the inclusion splits into two, provided that (1) 
always applies. But the narrowing may be such that 
R >> r while the ends are far away, in which case 
(1) will apply not only for all r less than a certain 
critical value r* (governed by a and by the posi- 
tion) but will actually become stronger as r de- 
creases. 

Precisely the same may occur for a dendrite 
(rod) surrounded by saturated solution; the only dif- 
ference is that the signs of the curvatures are re- 
versed, and hence also are the directions of the dif- 
fusion fluxes. This causes the dendrite to break up. 

The breakup time can now be evaluated. 


Conversion of an Inclusion 


Figure 2 indicates the scale of the region with- 
in which diffusion is effective; 1 is several times a. 
The part between QQ' and PP' is of length about 2a, 
and here the concentration changes from 


~ Co (1 — sre 


— 4 om 
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where c,, is the concentration above a flat surface, 


o is the surface tension, m is the mass of a mole- 
cule of the substance, p is the density of the solid, 
k is Boltzmann's constant, and T is absolute tem- 
perature. The density j of the diffusion flux is 
: De,,om / 4 4 
I~ DapkT. (Ge ar ) 
where D is the diffusion coefficient. 
We assume that all the material entering 
through PP' and MM' is deposited between these 
points; then r as a function of time t is given by 


dr De,,om / } 4 
“dt pk Tal (= = aM 
The limits of integration for r are ry (radius of 
initial section) and zero; neglecting the change in a, 
the time Tt needed is found as 


ekT atl a 
De,om” a rn (2) 


This result is discussed below. 


Conversion of a Dendrite 


Breakage of a Large Dendrite. Fig- 
ure 3 shows the shape of the weak point and the dif- 
fusion flows near it; here the diffusion occurs on 
the outside and so is governed mainly by the con- 
centration at infinity, which we assume to be de- 
termined by remote dendrites of radius about a, 
which is 


om 
ae cn(1 ae “ara 


The attack is most rapid where r is small 
(Fig.3), so we consider only a spherically sym- 
metrical solution to the steady-state diffusion equa- 
tion. Then 


The material balance gives us, as before, that 


=faser(t_ty 
fr) ’ 


dr 
pera r 


dt 
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oie Fon in rs. a Zea i (3) 
Detachment of a Branch. This is of 
especial interest as being the most common; it may 


be discussed, as above, by reference to Fig. 4, but 
here rp << a, so (3) simplifies to 


so 


ork Te 
a 3De,om (4) 
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Fig. 4 


Flat Dendrites. This case has no direct 
practical significance, but it reproduces the main 
features of the effects seen in Kliya's experimental- 
ly simpler system [3], in which ammonium chlo- 
ride dendrites grew in the space between two glass 
slides. The detachment of branches was recorded 
by cinemicrography; T = T (rp) was measured. This 
case is discussed as above. 

The concentration over the dissolving area is 


a (1 fn ar ) 


The material is deposited on surfaces whose mean 
radius is about d, because radii larger than this 
are impossible for plates a distance 2d apart. These 
surfaces may lie at a distance R > d; diffusion is 
confined essentially to two dimensions, so 


ip War d 
k a. 
ekTr In R 
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ee dr 
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But R >> d, r)<d, and the logarithm is aslow- 
ly varying function, so we have roughly 


c~ ion [= ln (1 — “1 ) fy 5 (2) In =. (5) 


Discussion 
Equations (2)-(5) all contain! 


e2?kT 
De,,om (6) 


This shows that the conversion is accelerated at 
high concentrations and by large diffusion coeffi- 
cients. 

The quantities in (6) are functions of tempera- 
ture, so T may be found as a function of tempera- 
ture. For example, c.. may increase rapidly, so 
the process accelerates toward high temperatures, 
especially if p decreases, as is usually the case. 
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Fig. 5 


‘These expressions cannot be deduced by the method of dimen- 


sions, because the above result for r is not the only combination 
of the quantities that has the dimensions of a time; but only Dcon- 
tains time, so we may put T ~1/D directly. 
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Comparison of (2), (3), and (5) shows that the 
conversion time has a logarithmic singularity as 
ry) tends to a or d, so the conversion is initially 
very slow if the initial deviation is small or if rjis 
comparable with the distance between the slides. 
Any narrow spot arising by accident will be leveled 
off if rp >> d; it will not become narrower, and no 
breakage will occur. 

Small ry (pronounced initial distortion) cause T 
to vary in different ways with the initial radius: 
slightly for inclusions, in proportion to r for bulky 
dendrites, and in proportion to r3 In rp for flat den- 
drites. 

Typical numerical values are D ~ 107 em? per 
sec;p ~ 3 g/cm’; o ~ 3 - 104yn/em?; m ~ 10 g; 
Co ~ 0.1 g/em*; T ~ 300°K; ry ~ 5° 107 cm; they 
give T ~ 4 min. 

Figure 5 compares (5) (full line) with Kliya's 
[3] results (points), the coefficient of proportion- 
ality being chosen to given R = 10d.? 

There are many instances where small crys- 
tals or cavities are altered in shape by surface ten- 
sion in a calculable way; experiment here can pro- 


vide reliable values for the surface tension as a 
function of crystallographic direction for such con- 
tacts. 

We are indebted to G.G. Lemmlein for advice 
and assistance in this work. 
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A technique devised in 1947-1950 is described, which now forms the basis for the 


commercial production of these crystals. 


These crystals (denoted in the industry by KT 
and those of ethylene-diamine tartrate (EDT) are 
used as piezoelectric plates of low temperature co- 
efficient of frequency. KT is better than EDT in 
that its elastic constants vary less with tempera- 
ture, and also because its electromechanical coup- 
ling coefficient is higher. This means that crystals 
of larger resonance range may be made. KT and 
EDT sometimes form good substitutes for quartz, 
e.g., in the filters of our 12-channel high-frequen- 
cy telephones. 

Mason [2] has described the properties of KT 
crystals in detail. 

These crystals have not found commercial use 
in other countries, because of insuperable difficul- 
ties in routine production; EDT is used, in spite of 
the advantages of KT, because the crystals are 
much more readily grown [38]. 

There appears to be no published literature on 
the growth of KT crystals. 

This salt (K,C,H,O,!4H,O) is made by neutraliz- 
ing the acid or potassium hydrogen tartrate with 
caustic potash or potassium carbonate. The crys- 
tals are colorless and readily soluble in water, 
density 1.98 g/cm’. The water of crystallization 
is very firmly bound; no effluorescence occurs 
even on prolonged exposure to temperatures up to 
100°C, but the water is lost completely around 165 


(it) 


°C. The water is not lost under vacuum, so the ma- 
terial may be used in evacuated devices. The crys- 
tals have little tendency to take up water; they are 
unaffected by air of normal humidity,but deliquesce 
at or above 80% relative humidity. 

The crystals are monoclinic (dihedral axial 
class); Fig. 1 illustrates one, the principal forms 
being the {100} and {001} pinakoids and the {111} 
and {111} dihedra, with the {010} monohedron also 
commonly present. The {020} monohedron is very 
rare; weak faces of the {101} and {101} pinakoids- 
are also seen, with occasionally a few dihedral 
faces with high indices, such as {113} : 

The polarity along the diad axis (the b axis) is 
seen as a difference in growth rates along that axis. 
The angle of the monoclinic cell is 89°10'. There 
is prominent cleavage on {100} and {001} faces. 
Figure 2 shows a KT crystal of the usual habit 
grown from a point seed. 

*Read at the First All-Union Conference on Piezoelectricity, 
November 24, 1952 [1]. 


Fig, 2 
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Aqueous KT is neutral; KT does not dissolve 
in alcohol or other organic solvents. The solubility 
in water is very nearly a linear function of tem- 
perature: 100 g of water at 25°C contains about 92 
g, as against 108 g at 50°C. This high solubility 
means that a slight change in water content (e.g., 
by evaporation) can give rise to appreciable super- 
saturation. The solution is of low viscosity. High 
supersaturations do not occur, because the com- 
pound crystallizes spontaneously at supersatura- 
tions substantially lower than those found for EDT. 
This adds to the difficulty of growing KT crystals. 

Slightly supersaturated solutions give crystals 
elongated along the b axis, the growth in the posi- 
tive direction being the more rapid. The (020) face 
grows the most rapidly and is usually lost, followed 
by the {111} dihedron. The {111} faces on the 
negative direction of the b axis grow at only half 
the rate of the {111} faces, while the rates for the 
{100} and {oo1} faces are lower by factors of 
about three. These values are average ones, be- 
cause temperature, supersaturation, etc., have 
some slight effect. Crystals grown from neutral 
solutions have very nearly the same habit as those 
grown from alkaline ones. 

The temperature affects the habit by acceler- 
ating the growth of the i100} and {001} faces, so 
the crystals produced by hot solutions are more 
nearly isometric, while those from cold solutions 
are highly elongated along the b axis. 

We began tests on the growth of KT crystals in 
1947, but for a long time obtained unsatisfactory 
results with pure solutions; the pinakoids hecame 
coated with a cloudy layer when the crystals were 
20-30 mm in size (Sometimes earlier), and this re- 
tarded the growth. This was followed by a network 
of small cracks, which then extended. It seemed 
that traces of impurities were damaging the crys- 
tals, since good results were obtained only with 
certain batches of salt, and then only during the 
first and second cycles (on account of impurity 
buildup). Some tests were done on constant-tem- 
perature growth, but no consistent results were ob- 
tained. 

It was clear that the crystals were very sensi- 
tive to traces of impurities, so various materials 
were added to the solution; alkali had been found to 
be of value with Rochelle salt and ammonium phos- 
phate, as it appears to suppress the effects of some 
impurities. A solution that had failed to produce 
any satisfactory crystals was treated with 50 g of 
KOH per liter; the very first run gave perfectcrys- 
tals, and further tests with these solutions contain- 
ing much KOH showed that the cloudiness and crack- 
ing in the pinakoids were completely suppressed. 
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Ammonia had the same effect as KOH, but KOH 
was preferred, because this had less tendency to 
alter in concentration with lapse of time. Poor re- 
sults were obtained with solutions containing NaOH, 
which reduced the growth rates of many faces, es- 
pecially {111} , and gave rise to cloudy films. 

The alkaline solutions were more stable than 
the neutral ones, and they allowed higher supersatu- 
rations, and hence higher growth rates; the maxi- 
mum supercooling without spontaneous crystalliza- 
tion or production of defective crystals was about 
two degrees. Growth rates of 5 mm/ day on the b 
axis are possible with these KOH solutions. 

Long-term production of these crystals has 
shown that the best KOH is 15-25 g per liter of so-- 
lution, because the KOH tends to reduce the solu- 
bility of the tartrate: 100 g of water at 25°C dis- 
solves about 85 g when the KOH content is 15 g per 
liter, and about 100 g at 60°C. 

Our first tests were done with point seeds of 
various orientations, but the best results were ob- 
tained with seeds elongated along the b axis with 
the fast-growing end facing outward; as would be 
expected. Unfortunately, the shape of the resulting 
crystals was unsatisfactory for cutting, on account 
of the excessive elongation and small area of the 
{100} faces. Filter crystals should either be of 
large cross section or have broad {100} faces. The 
first type may be produced satisfactorily from 
broad seed plates cut parallel to (111) or (020) 
faces, but the method has the disadvantage that the 
large area hinders regeneration. It is better to use 
rod seeds of much smaller area, and these allow 
one to produce crystals of large cross section in 
the two most useful styles. The first style has the 
length of the rod parallel to the edge between the 
(111) and (111) faces, which is analogous to the use 
of rods with ammonium phosphate, which are cut 
parallel to the intersection of adjacent faces of the 
dipyramid [4]. Square KT crystals of large cross 
section are produced. The second style is a rod 
parallel to the intersection of the (111) and (100) 
faces, which is held with the face of the (111) di- 
hedron upwards. This gives large broad crystals 
whose cross section is a parallelogram elongated 
along the c axis. These crystals are the best, be- 
cause their {100} faces are very prominent [5]. 
Figure 3 shows the best orientation of the rod seed 
and a crystal grown from such a seed. Figure 4 
shows a homogeneous crystal of weight about 900 g, 
which was grown from a seed 3-4 mm high, 2 mm 
wide, and 50-60 mm long. Stepped intergrowths 
are usually produced if the length of the seed devi- 
ates by more than 2° from the intersection between 
(100) and (111) faces. This tendency increases with 
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Fig. 3 : 


the length of the seed, so long seeds must be care- 
fully checked for orientation. The upper and lower 
faces must also be very smooth; grooves or re- 
cesses can give rise to intergrowths. The seeds 
must also be free from internal defects, because 
these may give rise to cracks. 

Vigorous stirring or movement of the crystal 
is essential to the production of homogeneous crys- 
tals; continuous circular motion of the crystal gives 
satisfactory results, but the growth rate should then 
be kept low, because the solution tends to stagnate 
on some faces. Reciprocating motion is better,e.g., 
rotation through 60-90° at 20-30 cycles per minute. 

Routine production is based on cooling solu- 
tions saturated at 45-50°. The temperature coeffi- 
cient of the solubility is not large, so the method 
requires a large ratio of solution volume to weight 
of finished crystal; a liter of solution saturated at 
50°C yields only 100-125 g of crystal on cooling to 
room temperature. Crystals grown by cooling from 
higher temperatures have considerable internal 
stress, so they often crack. 

The initial supercooling should not exceed 1°C; 
0.6-0.8°C is most usual. The crystals are sensi- 
tive to the supersaturation, so particular attention 
must be given to the sealing of the apparatus, other- 
wise loss of evaporation may produce excessive 
supersaturation and hence damage to the crystals. 

A supersaturation sufficient to produce a 
growth rate of 5 mm/day has a marked tendency to 
cause spontaneous crystallization; the most econom- 
ic rate for routine production is. 3-4 mm/ day. 

The crystallizers we have developed have been 
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described [6,7]. A single crystallizer will produce 
9-10 crystals each of weight about 500 g in a 20- 
day cycle. Crystals of weight over 1 kg are readily 
grown merely by using less seeds in the crystal- 
lizer. 
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The 24 groups of three-dimensional Patterson 
syntheses were first derived in graduate work by 
N.I. Shulepova at Gor'kii University in 1949, but 
the work was not published. Buerger [1] published 
a derivation in 1950, but with the omission of one 
group, which he corrected in the same year [2]. In 
1950 also, Kitaigorodskii [3] repeated Buerger's 
error and listed only 23 groups. 

In 1958, Lipson and Cochran [4] published the 
astonishing statement that the two-dimensional Pat- 
terson groups are derived from the 17 ordinary 
ones by the addition of a center of symmetry. 

All possible misunderstandings are removed 
if we define the Patterson groups as the product of 
the 11 Laue classes by the 14 Bravais lattices sub- 
ject to the condition (and this is where Buerger and 
Kitaigorodskii erred) that Laue class 3m may be 
oriented in two ways with respect to the hexagonal 
(primitive) lattice: the mirror planes may be per- 
pendicular to the shortest edges of the cell or par- 
allel to them. The complete list of the Patterson 
groups is then as follows: 


ee Patterson groups 
classes 
teed ‘Ny fee a 
EE Age 2. P2/m 3. C2/m 
Ill. mmm 4. Pmmm 5. Cmmm 6. Fmmm 7. Immm 
IV. 4/m 8. P4lm 9. T4/m 
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V. 4/mmm 10. P4/mmm 11. [4/mmm 
Vio3 NP, IED ANB, TRS as 
VII. 3m 44. P3m1 15. P31m 16. R3m 
VIII. 6/m 17. P6/m 

1X. 6/mmm 18. P6/mmm 

NILES 19. Pm3 20. Fm3 21. Im3 


. Pm3m 23. Fm3m 24. Im3m 


The same result is obtained if we say that the 
Patterson groups are the symmorphic space groups 
corresponding to the 11 Laue classes. 

The 7 two-dimensional Patterson groups are 
found in the same way: 


Me 2? dp, 

II. mm2 2. pmm2 3. cmm2 
III. 4 4. p4 

IV. Amm 5. pamm 

VEG 6. p6 

VI. 6mm 7. p6mm 
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The ionic refractions are given here for group 
Villa elements in their various states of ionization. 

The refractions have been calculated by 
Kordes' method from the ionic radii. 

It has been shown [2] that the ionic radius may 
be calculated from the electronegativity E of the 
atom, which Sanderson [8] gives as E=Z/ (47 Ej-r°), 
where Z is the number of electrons, Ej is the elec- 
tron density of an inert-gas atom with the same 
number of electrons, and ris the covalent radius. 

The covalent radius is dependent on the valen- 
cy for group VIII metals, so the atomic E of group 
VII metals cannot be based on the usual covalent 


radii. This difficulty can be avoided if we calcu- 
late the covalent radii from the atomic radii [4], be- 
cause the metallic and covalent bonds are essential- 
ly similar in nature, the differences being due to 
change in the coordination number (12 for metallic 
bonds, 6 for covalent bonds in these elements). 

Table 1 gives the covalent radii (from the 
atomic ones, with correction for the coordination 
number), the electron densities for isoelectronic 
inert-gas atoms, and the atomic and ionic E. The 
last are given by Sanderson's formula AE = 2.08n - 
(Eat)”“, where AE is the change in E when the atom 
loses n electrons. 


TABLE 1 
E for ion 
Element 7: E Eat 
+2 +3 4 
Fe 1.24 | 1.02 3.43 11.13 14.98 — 
Co 1.20 1.04 3.99 411.47 15.414 = 
Ni 1.19 1.06 3.74 11.79 15.80 19.87 
Ru 1.28 1.33 3.77 11.84 15.87 19.94 
Rh 1.28 1.34 3.82 11.95 16.01 20,08 
Pd 1.34 1.35 3.61 41.51 — 19.41 
Os 1.29 1.63 5.18 14.64 19.37 24.10 
Ir here 1.64 5.22 14.72 19.47 24 22 
Pt 1.32 1.65 4.90 14.11 _— 23,32 
TABLE 2 


—__—_—_——— eee 


Degree of ionization 
Element poke nas sa 

Ej r Ej r Ej Tr 
Fe 0.97 0.81 0.94 ORE: — = 
Co 0.99 0.80; 0.97 0.725 — = 
Ni 1.02 0.80 0.99 0.725 0.97 0.66; 
Ru gee? 0.86 AES 0.78 1.30 0.715 
Rh iLeSy 0.865 slay 0.78 kgeal 0.72 
Pd (| as 0.88 = — s Taps y 0.73 
Os 1.64 0.94 1.60 0.825 1.59 0.765 
Ir 1.62 0.94 1.64 0.82; - 1.60 0.765 
Pt 1.63 0.92, = — 1.64 0.78 
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TABLE 3 
emo SE ee 


Degree of ionization 


Element 

+2 | +3 | +4 
Fe 1.14 6 Be — 
Co 444 1,09 — 
Ni 1.08 1.06 1.03 
Ru 4237 1233 4220 
Rb 1,41 433 128 
Pd 1.52 — 434 
Os 4-57 1.41 1,28 
Ir 4.57 1,41 se) 
Pt 1.69 — 1.39 


The ionic E then give the ionic radii in group 
Villa as 


iZ } 


4 : 
ra TE Ej * Fion 


Table 2 gives Ej for Me**, Me**, and Me**, to- 
gether with the radii. 

These ionic radii have been compared with 
those given by Bokii [4], Ahrens [5], and Stockar 
[9]. 


The calculated radii lie within the limits of 


variation of the radii given by these various sources. 


Table 3 gives the ionic refractions calculated 
from Kordes' formula R = 1.354(r - g2/(n-1))4.5 , 
where r is ionic radius, Z is valency, and n is the 
Born repulsion factor. 
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These results agree closely with the empirical 
values Rett 1.15 cm?, Ratt 0.7 cm®, and 


3 


Ryyte 0.6 cm’, the mean deviation not exceeding 


0.3 cm’, which does not exceed the spread of the 
experimental results for the ionic refractions of 
the first triad of metals in group VII. 
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1. It has been shown [1] that the thermoelas- 
tic stresses in anisotropic bodies are very import- 
ant to the theory of annealing. 

Here we consider this for a very simple par- 
ticular case in two dimensions, with allowance for 
temperature terms. Here, the components u,v,w 
of the displacement vector are dependent only on 
the variables x and y, withu = u(x,y), v = v(x, y), 
w = 0, which gives us the components of the de- 
formation tensor as 


Ou | Ov 


to) Ov 
eg “u = ay? = og dy2 fy = Tyz = Yr = 9, 


(1) 
Hooke's law, as extended via the thermal-ex- 
pansion tensor @j;, is put as 


Sx = My Fx + AyySy + Ay39z + Aygtay + HT’, 
Sy = AyoFx + Age Sy + Ay3Fz + Agetxy + Ml’, 
Ez = Ayg9x + Ag3Sy + A335 + Agetay + Xgl", 


Tou = 68x + A2gSy + Agg%z + Aggtay + HoT. 
Then (see [2]): 


Ox = Byox + Bioty + Bretxy + Yul, 
fy = ByeFx + Booty + Posty + Yo2Z’, 
Yay = PreSx + BogSy + Booty + Y127"; 


4 
o, = — ae (43% x ++ AggFy -L- Azer xy + 331’). 


Here, Bij =a 4j3033/ a 33 Gi, j=1,2,6) 


oo 


413%ag A23%33 


A36%33 
438 


Ta = Sy 1 (by = 1 ify == a 
433 433 


(4) 


The equations of equilibrium and for the boun- 
dary conditions are as in the theory of elasticity 
without the temperature-dependent terms (we as- 
sume that external forces are absent): 


ie SH) 
ex a oy 0, 
do Ot 
aay 2} =.) 
Oy a a 2% 


(5) 
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3x COS (NZ) + Tz cos (ny) = 0, 
Ttxy COS (nx) + 3, cos (ny) = 0. 


(6) 
Compatibility of the deformations requires that 
Oe 


2 
lad ey 


RG 
aut T Ga8 = 


da Oy 


> 


and this, with (2), becomes 


3 
e) nay 
26 Ox2 


8 


n) 


°o,. 0c O27 e2T ac 
y xY x 
Bir Gye + Bra Gya + Bis Ger paler ae Bia 33 + 
eT do, as, Cie eT 
+ T22 G3 Bis ray Bes Say P88 Bz ay 112 9g dy 


(7) 


* * 
The substitutions 0x = 0x, —c4T, Oy = Oy —CooT, 
T xy = Txy Cy2T, in which cy;, Cy9, andcy, satisfy 


Yar = Birra + Brece2 + Bielie, 
Yaa = BreCir + BogCos + BesCra» 
Yoo = BieC11 + BooCos + BogCro, 


enable us to reduce the problem to the ordinary 
one of the stress distribution in a body subject to 
bulk and surface forces whose components are 

oT oT 
Sie Gi ie By? 


Lot OF aT 
713 Bie OPE ay’ 


Xx 


Y 
Xn = C4,T cos (nz) + cy9T cos (ny), 
Yn = Cy2T cos (nx) + Cy9T cos (ny). 


2. Consider the stress in a long cylinder of 
radius R cut from a cubic crystal along a cube edge. 
The result is to be applied in annealing theory, so 
we assume that the sides of the cylinder are heated 
or cooled at a constant rate th. The temperature 
T within the cylinder is given by [3] 


(8) 


BRIEF COMMUNICATIONS 


in which t is time and » is thermal conductivity. 

The coordinate system has its Z axis along the 
axis of the cylinder, the X and Y axes being in the 
plane of the cross section; then the tensor for the 
elastic constants simplifies to give [4] 


Qy = Qyy = Agg = Ag3, Ag = Ayg = Ajg = Ugg, 
Aig = Ag = Age = 0, Ogg = Ay. 


We seek a solution to (5) subject to (6) and (7) in 
the form 


3, = A(x? + 3y? — R?), 
oy = A(3x* + y? — R3), 
Txy = — 2Azy. (9) 
It can be shown that this choice is such as to satis- 
fy the equation of equilibrium and the boundary con- 
ditions, A being given by (7) as 


_ ah a,— Ag 


ax 6 (a? — a3) +2 (a4, — a3) + a3 


The solution implies stresses 0, given by (3) as 


2Aa ah 
6, = — — (27? — R*) + Ga (R? — r?), 


Teme oe 


{n addition, we must have oz = 0 atz=+H/2, in 
which H is the length of the cylinder, if the result 
is to be applicable to a finite cylinder with free 
‘ends; but oz is not dependent on z, so this condi- 
tion cannot be satisfied. The principle [5] allowsus 
to proceed as follows, however. We replace the 
stresses by the equivalent force and apply to the 
base of the cylinder a force equal in magnitude but 
opposite in direction to the above force. The state 
of stress in simple tension or compression is not 
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dependent on the elastic properties [6], so ox, Oy; 
and T xy are unaltered, while for 0, we have 


3, = (2r? — R) {=a ae an. 


ay 8xa, 


The radial and transverse components of the stress 
are 


3, = A(r? — R?), 
dg = A(3r? — R?), 
Tre = 0. 


The standard result is obtained automatically for 
an isotropic cylinder. 

We are indebted to V.R. Regel' for proposing 
the topic and for direction of the work. 
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These stresses often arise in technical applica- 
tions of crystals, and they have a major influence 
on the choice of growth and annealing conditions, as 
well as on the thermal stability in pyroelectric and 
piezoelectric devices. This makes it of interest to 
deduce the stresses produced in such plates by non- 
uniform temperature distributions.! 

The plate is taken to be cut in any way from 
the crystal; with it we associate a Cartesian coor- 
dinate system whose origin lies in the central plane 
and whose Z axis is normal to the surface. The 
temperature distribution is @ = 0(z), but otherwise 
it is entirely arbitrary. 

In terms of the components of the stress ten- 
sor, the equation for equilibrium is 


06, rs 


Bee 


(1) 
with 0), = 0 at the surface (no external forces); 
this gives 


5,2 =) 


ip sh, hy 
(2) 
for all points. 

The stresses and bending moments averaged 
over the thickness are zero at the edges: 


34; = 0, 203; => 0, @. ji => il. 2y, (3) 
The bar in (3) denotes averaging over the thickness 
in the usual way: 


+a 
\ 94 (z) dz, etc., 


—a 


1 
FO ae 


in which a is half the thickness. 
The solution is a function of z alone, so (3) is 
satisfied throughout the plate. 
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Hooke's law in generalized form? gives 
Sig = Signet Ter + %i79. (4) 


and with (8) we have 


ij = 00, Bey = 420, f =1, 2). (5) 


Only three of the compatibility equations are 
not satisfied identically: 


Oe; 5 ee 
at =0 (i, j=1, 2). (6) 


From (5), we have the solution to (6) as 


\ 


which gives us a system of linear equations for the 
stresses: 


RET Re A Ply 
Si = ai;( 520-2 +8) (i, 7=1, 2), (7) 


Sijkl Sxl = %ij ler. 20-2 +5—6) ei eee (8) 


The solution to (8) may be put as 


on = Bia (5 20-2 i 6—6) (Z, = il. Ze (9) 


in which the thermoelastic coefficients 6;, are as 
follows, if we use the ordinary rules [2] for convert- 
ing the subscripts: 


(10) 


1No thermoelastic stresses are produced by uniform distributions, 
or ones linearly dependent on the rectilinear coordinates, as for an 
isotropic body [1]. 
Here and henceforth we assume summation with h 

with respect to the 
subscripts as appropriate, 
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Here, 


Siz Sin S16 
N=]S12 S22 Sag 


Sig Soe See 


and Amn is the algebraic complement of element 

Smn in A. We have z@ = 0 for the special case of a 

distribution symmetrical about the mean plane, and 
six = Bir (U — 9). (11) 

Table 1 gives formulas for the fj, for plates 
cut normal to the physics axes in terms of the usu- 
al settings® of the crystals. 

For instance, for ethylenediamine tartrate ha a 
monoclinic crystal widely used in piezoelectric de- 
vices, we have for plates cut perpendicular to the 
X, Y, and Z axes (in dynes/ cm” - deg),° 


Byy = 105-107, See = 0.63-107, Bry = 0.05- 107, 
B22 = 0.99-10%. 3.2 = 0.7 +107. Byy = 0.56- 107. 
By: = 9. Box = 0.16-107. Bxy = 0. 


The stress ellipse degenerates to a circle for 
a plate cut normal to an axis of order three or 
higher, as would be expected. 

This case is closely similar to that for an iso- 
tropic plate having forces parallel to the surface 
and equal [5] to 


7= (2-24 5-0), (12) 


vA 


where 


ae aly 


l[—u 


E being Young's modulus and uw Poisson's ratio. 
Formula (12) is the common basis of calcula- 
tions on heating and cooling [6], annealing [5], 
quenching [7], and thermal stability [8] for isotrop- 
ic plates, so (9)-(11) enable us to calculate these 
effects for any crystalline plate. Of course, allow- 


ance must be made for the anisotropy in the therm- 
al conductivity [9]; the heat transfer will be gov- 
erned by the conductivity normal to the surface. 

All the above formulas are applicable to any 
plate with rectilinear [10] anisotropy, not merely 
to crystal plates. 
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Natural quartz crystals are usually twinned, 
the dauphiné and Brazilian laws being the most 
common, one or another of which accounts for a 
large part of almost any natural crystal. 

A Brazilian (optical) twin consists of right- 
and left-handed parts in antiparallel positions, the 
positive rhombohedron of the right-handed part be- 
ing intergrown with the same of the left; the same 
applies to the negative rhombohedra. 

The symmetry of such a twin, if ideal, is 6-m, 
which is higher than that of quartz itself. An ideal 
Brazilian twin shows no piezoelectric response, 
because all the piezoelectric moduli for this class 
are zero [1]; plates cut from such crystals thus 
find no piezoelectric applications. 

A feature commonly overlooked is that the 
elastic-constants tensor for such a twin, {Sil aa ; 
is exactly analogous with the tensor for a single 
crystal of quartz, which has 3:2 symmetry. Thus, 
the mechanical modes of a plate cut from a twin 
will be as for such a crystal; in particular, the 
resonant frequencies and temperature coefficients 
of these will be the same [2]. 

Brazilian twins commonly occur as polysynth- 
etic intergrowths (alternating thin layers of the 
two forms); the two components frequently differ 
very much in size, so such an intergrowth contains 
areas in which the twinned parts are small, al- 
though the size of such areas may be inadequate to 
allow the cutting of truly single-crystal plates of 
the required size. 

It will be clear from the above that such a 
twinned crystal may be cut into resonant plates 
whose properties are precisely as for single-crys- 
tal ones; this is possible if the twins lie in a part 
of the plate not covered by the metal electrodes. 
This is often so, because most quartz plates are 
nowadays only partly metallized; in particular, the 
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electrodes of high-frequency plates of AT and BT 
cuts are usually placed as shown in Fig. 1. 

To test these ideas, we prepared a set of 38 
plates containing varying proportions of Brazilian 
twins, and also a set of 25 plates cut from single 
crystals, all cuts being AT (35°30' to the Z axis) 
and nominally of frequency about 5 Mc. These 
plates could also be excited at their third harmonic 
(size 20 X 15 x 0.3 mm). The electrodes were as 
shown in Fig. 1b (spot electrodes 8 mm in diam- 
eter). 

The twins usually lay toward the edges of the 
plate and were not covered by the metal, but in 
some cases the thin twinned parts passed under the 
electrodes. These plates were mounted in K-19 
holders, which were evacuated to 107! to 1072 mm 
Hg. 

The plates were tested as follows: 

1. The activity was measured at the funda- 
mental and at the third harmonic with a Shembel 
oscillator. 

2. The temperature coefficient of the frequen- 
cy was measured between +20 and +70°C with the 
plate excited at its third harmonic. 

3. Parameters » (circuit magnification) and 
Cy (static capacitance of crystal and holder) were 
measured [38]. 
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Fig. 2, Examples of disposition of twins on resonant 
plates whose response does not differ from that of a 
single crystal. 


4. The spectrum was recorded. 
Rg (equivalent loss resistance) was also meas- 
ured. 


The twins produced no change in the frequency 
coefficients. These plates had been ground in the 
same mount with the single-crystal ones to the 
same size and had the same frequency. Also, all 
twinned plates (except one) had temperature coeffi- 
cients equal to or less than those of the single- 
crystal plates. The one poor plate was checked 
with an x-ray goniometer and was found to have 
been wrongly cut, so we may say that Brazilian 
twins produce no appreciable distortion in the mode 
of oscillation. Moreover, the twinned plates in- 
cluded none vibrating at several frequencies. 

The activity (as expected) was found to be di- 
rectly related to the amount of twinned material 
under the electrode; it was as for single-crystal 
plates when there was no twinned material under 


Fig. 3, Examples of twins in plates 
whose response was well below nor- 
mal. 


the metal, or even when the twinned part under the 
metal was only a thin band (Fig.2). The same ap- 
plies to x and Rg, which were also dependent on 
the particular part of the twin found under the 
metal (Fig. 3). 

Quartz crystals containing Brazilian twins, 
which are usually rejected for piezoelectric use, 
may thus often be employed to make filter plates 
without showing performance inferior to that of 
single-crystal ones. 

We are indebted to I.S. Zheludev for valuable 
advice and for criticisms of the manuscript. 
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This demonstration is prepared as follows. A 
finely ground and slightly warmed glass slide is 
touched with a small brush previously dipped in 
molten salol (melting point 42.5°C) to produce two 
fairly thin lines meeting ata small angle. Thespe- 
cimen is then heated again and is allowed to cool to 
room temperature. Salol readily becomes super- 
cooled, so it does not usually crystallize spontane- 
ously. 

The demonstration is started by initiating crys- 
tallization at one of the four ends with a needle 
bearing solid salol; the preparation is at once trans- 
ferred to the microscope (used with ordinary light 
at a magnification of about twenty). The crystals 
grow rhythmically and (an interesting point), the 
growth front bends at the junction of the streaks to 


end up perpendicular to all four branches (Fig. 1). 
This can be demonstrated to a large audience by 
using a projection microscope, but this requires 
care to insure that the preparation does not become 
heated by the illumination system, because this can 
halt the rhythmic growth. 

We have previously given a qualitative explana- 
tion for the rhythmic growth of salol crystals [1,2] 
which has not previously been described in this 
form. 
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Fig. 1. The crystallization was started in the top left branch and propagated in the directions shown by the arrows. 
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Spiral growth layers have recently been re- 
ported for many organic and inorganic substances, 
but none of these is a ferroelectric. Recently, we 
observed spirals on barium titanate, which were 
produced during thermal analyses concerned with 
the crystallization of BaTiO3. 

The layers occur on the most abundant faces of 
the {100} form (Fig. 1); they have straight edges. 
Figure 1 is a crystal seen in slightly oblique re- 
flection, which reveals simultaneously the spirals 
and the domain pattern. The parts of the spiral are 
not parallel to the edges of the crystal or to the do- 
main boundaries, which run parallel to those edges; 
the deviations are 5-6°, The reason for this maybe 
as ae the spirals on salol crystals [1]. The angu- 
lar deviation @ is governed by the ratio of the velo- 
city v of the tangential displacement to the velocity 
V of the elongation of the periphery of the spiral. 

Oblique illumination (Fig. 2) was used to re- 
veal the domains directly under the center of the 
spiral. The crystal evidently contains several re- 
gions differing in orientation of the domain boun- 
daries directed along face diagonals. Directly under 
the center of the spiral lies an opaque region resem- 
bling a spherical inclusion; near the spot and along 
the diagonals lie regions in which the domains with 
mutually perpendicular directions overlap. The 
four sectors around the center have the directions 
of the domain boundaries changing in a way such as 
to make the boundaries in adjacent sectors mutual- 
ly perpendicular, 

The final domain structure is dependent on the 
initial stresses and on any change in these pro- 
duced by the domain structure [2], so we may sup- 
pose that the structure seen here is due to the 
stresses arising near the screw defect, whichacts 
as a growth center. 
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Chernysheva [3] has described related effects 
in the twinning of Rochelle salt; the distribution of 
the twins (domains) is governed by the structural 
imperfections, and this becomes more complex in 
parts rich in inclusions (microscopic foreign par- 
ticles, inclusions of solution). Polarized light re- 
veals a regular array of light and dark regions 
around the inclusions, but the photographs show 
that the boundaries near the twins (light or dark) 
all have the same direction. This domain structure 
in Rochelle salt differs from that in our barium ti- 
tanate crystal. 

The BaTiO; crystals were platelets with faces 
of the {100} form, which were produced from a 
melt of barium chloride, barium carbonate, and 
titanium dioxide (2:1, 4:1), which was used in 
corundum crucibles of 20 em® capacity at 1400- 
1480°C. The temperature was measured with a 
differential platinum — platinum — rhodium thermo- 
couple. 


Fig. 1. Face of a barium titanate crystal 
Carrying a spiral growth layer, x 150. 
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Fig. 2. Internal structure of the same 
crystal, x 150, 


The temperature and time spent in the molten 
state were tested for their effects on the heat of 
crystallization and on the morphology of the crys- 
tals. The differential cooling curve was smooth 
and without halts if the melt was heated to tem- 
peratures up to 1450°C, but not maintained at the 
maximum temperature, because only a small quan- 
tity of material was crystallizing, so the tempera- 
ture difference from the Al,O3 standard was inap- 
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preciable. An effect became apparent only when 
the melt was kept at 1450°C for 3 h, presumably 
because sufficient titanate is formed in solution to 
give a fair crop of crystals on cooling, these re- 
leasing heat and so delaying the fall in temperature, 
which is seen on the differential curve as an exo- 
thermic peak. These were the conditions that gave 
crystals with spirals on their faces. An effect is 
detectable for shorter times (about 1 h) but only if 
the temperature is raised to 1480°C. 
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These structures on quartz crystals have been 
reported several times, but without detailed illus- 
tration or description. 

Willis, in 1952 [1], merely stated that he had 
observed these spirals on the R faces of natural 
quartz, while Weill [2] gave a photomicrograph of 
a spiral on a (0001) surface of a quartz inclusion in 
a molybdenite block from Queensland (Australia). 
This spiral is very difficult to distinguish in the 
journal, partly because the state of the surface 
generally was poor. More recently, it has been 
reported [3] that the accessory growth on the R 
faces of artificial quartz consisted of concentric 
layers and also spiral structures; but the plates 
show only concentric layers and include no spiral 
structure. 

Here we give several photomicrographs of ac- 
cessory growth spiral structures on the (0001) sur- 
faces of an artificial quartz crystal. 

This crystal was grown in this institute from 
a seed plate in the form of a disk 2.5 mm thick and 
35 mm in diameter cut perpendicular to the L; axis 
of a single crystal. The initial process was one of 
regeneration; the circular seed plate rapidly be- 
came hexagonal, and the area of the slow-growing 
R faces gradually increased, whereas that of the 
fast-growing basal surface rapidly diminished. The 
growth rates of the faces of the small r rhombo- 
hedron were substantially higher than those of the 
R faces under these conditions {they were roughly 
the same as for the basal surface), so the r faces 
did not arise during growth. At the end of regen- 
eration the crystal had faces corresponding to three 
simple forms: {1010} hexagonal prism, {1011} 
major rhombohedron, and {8191} trapezohedron. 
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The basal surface grew by nucleation and ex- 
tension of numerous conical accessories, which 
often became trigonal pyramids with rounded faces 
and edges, the orientation being such as to set the 
rounded edges toward the emergence of the nega- 
tive ends of the Ly axes (Fig. 1). 

It appears that amorphous silica may be pre- 
cipitated from the soda solution while the auto- 
clave is cooling during the final stage of growth; 
this forms a thin film over the surfaces that tends 
to mask details. Figure 2c shows parts of sucha 
film. This may be adequately removed from the 
prism and positive rhombohedron faces by boiling 
in concentrated NaOH, but this treatment was in- 
effective for the basal surface, except in the case 
of crystals at the very top of the autoclave, which 
remain above the solution during the cooling. Very 
occasionally such crystals have basal surfaces free 


di 
/ 
m™ / 


bs Aged 
m(i100) 
| a 
m7 7 


a 


Eigen 


BRIEF COMMUNICATIONS 481 


Fig. 2. Accessory pyramids with obvious spiral structure on a (0001) surface of artificial quartz, x 100. 


from this film, and one such crystal showed spiral 
growth on the basal surface; three spirals are 
shown in Fig. 2, magnification 100. 

The positive print reveals the trigonal struc- 
ture of the spirals and shows that the (0001) layers 2. 
are anisotropic in growth rate. The spirals of Figs. 
Za and 2b are adjacent on a single crystal (right- 
handed); their senses are opposite, which is evi- 
dence that the sense of the spiral is unrelated to 
the optical sign of the crystal. 

We are indebted to N.V. Gliki for advice and 
assistance in photographing these spirals. 
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The Crystallography Symposium in Madrid 


(April 2-6, 1956) ended with excursions to Southern 
Spain, including centers of the former Moorish cul- 


ture. Some of the relics of Seville (Alcazar), 


Granada (Alhambra), and Cordoba include walls 

covered with simple geometric patterns, many of 
which form striking illustrations of the 17 (single- 
color) and 46 (two-color) infinite two-dimensional 
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Fig. 5 Fig. 6 


Symmetry groups which are the most rich in sym- 
metry but which, for some reason, seldom appear 
in tapestries and other patterns. 

The striking hexagonal pattern at Alhambra 
(12th to 13th century) consists entirely of current 
international symbols for hexad, triad, and diad 
axes (Fig. 1); itillustrates group p6. 

Figure 2 shows a tetragonal pattern correspond- 
ing to the p4'g'm two-color group and consisting of 
a single, very simple figure repeated. The same 
pattern is seen in the border of Fig. 1. 

Figure 3 would illustrate group p6' if the al- 
ternately black and white symbols for the triad 
axes were identical, i.e., if there were only stars 
of two colors or only ditrigons. The difference be- 
tween the figures reduces the symmetry to p3. 

Figure 4 represents interwoven bands of sym- 
metry P622, i.e., one of the 230 (80) groups. 

Figure 5 has the p6mm symmetry reduced to 
orthorhombic by the use of two halftones. The 
border at the top is common in tapestries and has 
symmetry pa'. 

The tetragonal patterns of Figs. 6-8 are quite 
common in parquet flooring and elsewhere. 

The stone floors of the Moorish buildings most 
frequently have the herringbone pattern character- 
istic of our wood-block floors; they have the barely 
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Fig. 10 


perceptible orthorhombic symmetry pgg (Fig. 9). 
All the older buildings have decorative works 
often in the form of regular reentrant polyhedra 
such as those illustrated in volume 27 of the Great 
Soviet Encyclopedia. However, these are more 
commonly based not on the regular dodecahedron 
but on combination of the cube with the octahedron 
(the coordination polyhedron giving closest packing) 
to give alternating square and triangular faces, 
these bearing the corresponding pyramids. 
Finally, Fig. 10 illustrates the "Symmetry" 
sculpture in the Prado Museum, Madrid. 
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In my paper of the above title I stated that this 
error occurs in a variety of forms in textbooks on 
crystallography; it is essentially the impermis- 
sible replacement of an asymmetric figure or point 
by the most symmetrical figure, i.e., a sphere or 
a point by implication assigned the symmetry of a 
sphere. 

Ansheles criticizes not my entirely clear for- 
mulation, but some idea entirely of his own, which 
is concerned not with the replacement of an asym- 
metric figure by a symmetrical point, but with the 
replacement of a symmetrical figure by some un- 
specified point. The same applies to the comments 


by the editor, in which he speaks of some crystals 
as though shown in my illustrations, whereas the 
latter do not show any crystals. Ansheles's com- 
ments, and those by the editors of the Zapiski, are 
irrelevant to my paper, so they require no addi- 
tional reply from me. It remains only to note that 
such comments could appear in a journal published 
by the Academy, whereas the objections to them 
were not published. 


*Zap. Vsesoyuzn. Min. Obshch., 84, 239 (1955). 
{ Zap. Vsesoyuzn, Min. Obshch., 83, 408 (1954). 
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